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(57) Abstract 

The present invention provides a method for detecting and quantifying mRNA in a sample. The mRNA that can be detect- 
ed has a unique sequence. The method includes immobilizing a First polynucleotide to an insoluble support. The first polynucleo- 
tide has a first sequence that hybridizes to the unique sequence on the mRNA. After immobilization of the First polynucleotide, 
the sample is applied to the insoluble support under conditions that allow the unique sequence on the mRNA to hybridize with 
the first polynucleotide. Thereafter, a labeled probe can be applied to the insoluble support. This labeled probe has a sequence 
thereon that hybridizes to a portion of the mRNA other than the unique sequence. The application of the labeled probe is per- 
formed under conditions that allow the labeled probe to hybridize with mRNA immobilized on said support, if present. Finally, 
the amount of the labeled probe immobilized on the support is measured to provide an indication of the amount of mRNA pres- 
ent in the sample. 
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METHOD FOR MEASURING MESSENGER RNA 

Background of the Invention 

Field of the Invention 

This lnv^nHn" 1 c r^lar^H rn o m<»tk^^ f ™ ^ _ r>v * ' — t-> v ' * * - 

-i — ..it-.j.ju.. i- . *..~* w u u i. ■ >- ». * «^\-" *wi iiivujuitii^ * i — j t^ci ixi>^\ ^ 1 1 1 K l \ a\ | dnu 

5 polynucleotide-immobilized supports useful in this method. Specifically, this invention is 

related to measuring methods and reagents in which various types of mRNA can be 
measured at the same time. In one embodiment of this invention, the mRNA need not be 
purified from the cell, and the mRNA can be rapidly quantified with high sensitivity. 
Description of Previous Techniques 

10 Previously well-known techniques for measuring mRNA levels include the Northern 

Blot Method (See, e.g., J. Sambrook, et al.. Molecular Cloning, 2nd ed„ Cold spring Harbor 
Laboratory Press, 1989 hereafter referred to as "Molecular Cloning"), the Dot Blot Method 
(See, e.g., Molecular Cloning), Ribonuclease Protection Assay (See, e.g.. Molecular Cloning), 
and the Reverse PCR Method (H.A. Erlich (ed.). PCR Technology-Principles and 

15 Applications for DNA Amplification, Stockton Press, New York, NY, 1989). The disclosures 

of all of these previous references are hereby incorporated by reference. 

The Northern Blot Method detects mRNA which is separated by molecular weight 
using electrophoresis. The separated molecules are immobilized on a membrane and 
mRNA levels are determined by hybridization of marker DNA probes. With this method, 

20 the molecular weight of the subject mRNA must generally be known. By using this method 

a complex between a labeled probe. mRNA and an insoluble membrane is formed. The 
amount of label in the complex at the appropriate site on the membrane is then measured 
coiorimetrically or by radioactive emissions. 

4 

The Dot Blot Method detects mRNA by hybridization with a labeled probe after 
25 purified mRNA is immobilized as a spot on a membrane. This method can be used to 

rapidly examine many samples at one time. However, both the Dot Blot and the Northern 
Blot methods require the use of a specific labeled probe for each mRNA which is desired 
to be detected. 

The Ribonuclease Protection Assay method is used to detect mRNA after it is 
30 initially purified. The purified mRNA is hybridized with a labeled RNA probe, and then 

treated with a ribonuclease that specifically digests single stranded RNA. The labeled RNA 
probe will form a duplex in the solution and be protected from digestion. This is then 
immobilized on a membrane after separation of the double-stranded RNA/probe complex. 
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The reverse PCR method is used where cDNA is synthesized using reverse 
transcriptase and purified mRNA. This cDNA can be immobilized to a membrane. The 
synthesized cDNA is then amplified using standard PCR (Polymerase Chain Reaction) 
methods. After separating cDNA according to its molecular weight by electrophoresis, it 
5 becomes a single band. The DNA is then immobilized on a membrane, and detected by 

hybridization with a labeled DNA probe. However, after PCR, no quantitative information 
remains due to the amplification of the original message. 

All of the foregoing methods take extended periods of time to perform since the 
mRNA used must be in at least a partially purified form. Also, electrophoresis of the 
10 samples and immobilization to a membrane is necessary 1 . Both of these steps reduce the 

reliability of quantification of the mRNA. 

Thus, there remains a need to resolve the above problems, and to provide a method 
of measuring for rapidly determining messenger RNA levels and a measuring reagent useful 
in such a method. 
15 Summary of the Invention 

In the mRNA assay methods of the present invention, a first and a second 
polynucleotide probe can be used. The inventors have discovered that immobilization of 
probes (a first polynucleotide) rather than samples, on a water-insoluble support can provide 
a basis for the rapid quantification of multiple unpurified simples containing mRNA. 
20 Applicants have provided an efficient method of identifying appropriate specific probes for 

this purpose. In preferred forms of the present invention, the first polynucleotide probe has 
a nucleotide sequence complementary to a region which is unique in the sample to the 
mRNA and does npt bind other portions in the target mRNA. The inventors also 
discovered that a single kind of second polynucleotide probe can be used for simultaneous 
25 quantification of various mRNAs. However, any of a variety of probes that bind to a 

different region of the mRNA than the. first probe can be used for this purpose. In this 
regard, the inventors have found that a second polynucleotide probe containing sequences 
complementary to a polyadenylic acid tail of mRNA can be used in such a simultaneous 
assay method. Based on these and other discoveries, the inventors established the present 
30 invention. 

One embodiment of the present invention is a highly sensitive, quantitative and rapid 
method for detecting and quantifying mRNA in a sample without the need to purify mRNA 
from cells. In this embodiment of the invention, the method comprises the following steps: 
(a) identifying a polynucleotide sequence that is unique to the mRNA: 
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(b) immobilizing a first polynucleotide to an insoluble support, the first 
polynucleotide having a first sequence complementary to the sequence unique to the 
mRNA; 

. (c) incubating the sample with the insoiubie support under conditions 
wherein the unique sequence will hybridize with the first polynucleotide, thereby 
immobilizing mRNA present in the sample to the insoluble support: 

(d) washing non-immobilized components of the sample from the 
insoluble support; 

(e) labeling mRNA on the support in a manner that label is incorporated 
onto the support related to the amount of mRNA on the support: and 

(f) measuring the amount of label immobilized on the support. 

In a preferred embodiment, the above step of labeling mRNA on the support involves 
incubating a labeled probe, most preferably poly-d(T). with the insoluble support under 
conditions wherein the labeled probe will hybridize with mRNA on the support, and the 

15 labeled probe bears a label and is complementary to a part of the mRNA other than a part 

of the mRNA complementary to the first probe, thereby immobilizing to the insoluble 
support the label on the labeled probe which has hybridized with the mRNA; and thereafter 
washing any non-immobilized labeled probe from the insoluble support. 

In addition, the above method of labeling nucleic acids in step (e) can preferably 

20 involve labeling the oligonucleotides with a nucleic acid stain, most preferably from the 
group consisting of ethidium bromide, yoyo-1 and toto-1. 

In a further preferred embodiment, the above method can utilize a label, wherein 
the label comprises a radionuclide, or alternatively biotin. In the method using biotin, the 
incubation can preferably be with enzyme-conjugated streptavidin: followed by measurement 

25 of the amount of enzyme bound to the insoluble support. In this embodiment step (f) can 

additionally comprise adding streptavidin bound to alkaline phosphatase to the insoluble 
support. 

In still another embodiment of the preferred method, the mRNA can code for a p 
receptor, and wherein the first polynucleotide comprises the sequence of SEQ ID NO:6. 
An alternative embodiment of the above method includes mRNA coding for the a subunit 
of Gs protein, wherein the first polynucleotide comprises a sequence complementary to a 
sequence selected from the group consisting of SEQ ID NO:l. 145-146. 209-212. 268-287, 
705 and 708. 



30 
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Yet another alternative mRNA sequence that can be used with the present invention 
method codes for the Gi-2 protein, wherein the first polynucieotide comprises a sequence 
complementary to a sequence selected from the croup consisting of SEQ ID NO:703, 806, 
141 143. 187-197 and 238*253. 
5 The present method can also preferably involve the use mRNA coding for substance P 

receptor, wherein the first polynucleotide comprises the sequence of SEQ ID NO:5. 

mRNA encoding Gi-3 protein can also be used in the present method, wherein the 
first polynucleotide comprises a sequence complementary to a sequence selected from the 
group consisting of SEQ ID NO:704, 707, 144, 198-208 and 254-267. 
10 Still another embodiment of the above method can involve the mRNA coding for 

the GM protein, wherein the first polynucleotide comprises a sequence complementary to 
a sequence selected from the group consisting of SEQ ID NO:702. 805. 141, and 224-237. 
An even further preferred embodiment includes mRNA coding for G protein, wherein the 
first polynucleotide comprises a sequence complementary to a sequence selected from the 
15 group consisting of SEQ ID NO:700, 701, 139, 140. 148-159, 162-172. 176-186 and 3-4. 

Most preferably, the first sequence hybridizes to a unique sequence of the mRNA 
and not to other polynucleotides present in the sample, and even more preferably, the 
sample comprises a cell lysate. 

In an alternative preferred embodiment, the label in the above method can be 
20 measured by light emitted therefrom, and step (f) of the method comprises measuring the 

amount of light emitted by the label. The method of measuring the light can preferably 
involve; recording the amount of light on film: and measuring the exposure of this film using 
a densitometer. In addition, the light measurement can alternatively comprise adding 
ATTOPHOS and measuring fluorescence emitted using a fluorimeter. 
25 The present invention method in step (a) above can advantageously include a 

computer program, with the computer program preferably involving use of an H-site model. 
The H-site model of step (a) is preferably used to identify the first polynucleotide probe of 
the method by the steps of: 

specifying a minimum melting temperature for the first nucleotide 
30 probe and the nucleotide sequence specific to the organism: 

specifying a nucleation threshold that places a minimum value on the 
number of base pairs at any nucleation site; 
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determining the melting temperatures (T m ) of the first nucleotide 
probe and the sequence specific to the organism at every possible 
hybridization point; and 

selecting the nucleotide probe having the highest T m value. 
5 The computerized method of identifying the first polynucleotide probe can preferably 

have the melting temperature determined by the formula: 

T m = 81.5-16.6(log[Na])-0.63%(formamide) + 0.41(%(G + C)).600/R wherein Log[Na] is the 
log function of the sodium concentration. 0.063%(formamide) is the concentration of 
formamide, %(G + C) is the percentage of matched GC base pairs, and N is the probe 
10 length. 

Additionally the method of step (a) can be used to identify the labeled probe 
referred to above using the H-site computerized model by the steps of: 

specifying a minimum melting temperature for the labeled probe and 
the nucleotide sequence specific to the organism; 
15 specifying a nucleation threshold that places a minimum value on the 

number of base pairs at any nucleation site: 

determining the melting temperatures (T m ) of the first nucleotide 
probe and the sequence specific to the organism at ever)' possible 
hybridization point; and 
20 selecting the nucleotide probe of the proper length having the lowest 

T m value. 

A method of identifying the second nucleotide probe can involve determining the 
melting temperature by the formula: 

T m = 81.5-16.6(log[Na])-0.63%(formamide) + 0.41(%(G + C))-600/N. wherein Log[Na] is the 
25 log function of the sodium concentration, 0.063%(formamide) is the concentration of 

formamide, %(G + C) is the percentage of matched GC base pairs, and N is the probe 
length. 

A different embodiment of the present invention is a method of inhibiting RNase 
activity in a sample containing RNA in a buffer system containing SDS and EDTA, 
30 comprising adding vanadyl ribonucleoside complex to the sample. 

Preferably, there is at least a five-fold molar excess of EDTA over RNA in the sample and 
most preferably, the vanadyl ribonucleoside complex is added in combination with proteinase 
K. In an alternative preferred embodiment, the vanadyl ribonucleoside complex is added 
in combination with RNasin. 
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Yet another different embodiment of the present invention is a polynucleotide- 
immobilized support that is useful in detecting and quantifying the amount of (J receptor 
mRNA in a sample, comprising: 

an insoluble suppoi i, 

5 a first polynucleotide bound to the support, the first polynucleotide 

comprising a sequence that is complementary to a sequence unique to p receptor 
mRNA. 

More preferably, the sequence unique to p receptor mRNA comprises a sequence selected 
from the group consisting of SEQ ID NO:6. 
10 Still another preferred embodiment of the present invention is a polynucleotide- 

immobilized support that is useful in detecting and quantifying the amount of jun oncogene 
mRNA in a sample, comprising: 

an insoluble support; 

a first polynucleotide bound to the support, the first polynucleotide 
15 comprising a sequence that is complementary to a sequence unique to jun oncogene 

mRNA. 

Preferably, the sequence unique to the jun oncogene mRNA comprises a sequence selected 
from the group consisting of SEQ ID NO:ll-35. 30-53. 62-138. 310-342. 600-606, 614-620, 
and 628-634. 

20 Another embodiment of the present invention is a polynucleotide-immobilized 

support that is useful in detecting and quantifying the amount of G protein mRNA in a 
sample, comprising: 

an insoluble support; 

a first polynucleotide bound to the support, the first polynucleotide 
25 comprising a sequence that is complementary to a sequence unique to G protein 

mRNA. 

Preferably, the sequence unique to G protein mRNA comprises a sequence selected from 
the group consisting of 1, 2, 700-709, 805. 806, 139-159. 162-172 and 176-309. 

Another preferred embodiment of the present invention is a polynucleotide- 
30 immobilized support that is useful in detecting and quantifying the amount of Substance P 

receptor mRNA in a sample, comprising: 
an insoluble support: 
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a first polynucleotide bound to the support, the first polynucleotide 
comprising a sequence that is complementary to a sequence unique to Substance P 
receptor mRNA. 

Preferably the sequence unique to Substance P receptor mRNA comprises the sequence of 
5 SEQ ID NO:5. 

Brief Description of the Figures 
Figure 1 illustrates an immobilization procedure of a first polynucleotide probe onto 
an insoluble support by the maleimide method. 

Figure 2 illustrates an immobilization procedure of a first nucleotide probe onto an 
10 insoluble support by the carbodiimide method. 

Figure 3 schematically shows the principle of an mRNA quantification method in 
accordance with the present invention. 

Figure 4 (I) shows a photographic image of an X-ray film exposed to 
chemiluminescent light derived from the human p 2 adrenergic receptor cDNA probe on a 
15 Northern membrane, in which various concentrations of the human f* 2 adrenergic receptor 

mRNA were immobilized. (II A) shows a photographic image of an X-ray film exposed to 
chemiluminescent light derived from the oligo-(dT) probe in accordance with the present 
invention in which various concentrations of the human P 2 adrenergic receptor mRNA were 
trapped onto the first oligonucleotide-immobiJized plastic plate. .(Il.B) shows a densitogram 
20 of the image shown in (II.A.). (II.C) graphically depicts the relationship between relative 

intensity of chemiluminescent signals and concentrations of mRNA of the image shown in 
(II.A). 

Figure 5 graphically depicts the linear relationship between relative intensity of 
chemiluminescent signals and concentrations of mRNA in a separate experiment from that 
25 shown in Figure 4.11. 

Figure 6 graphically depicts a calibration line for the rat Gs protein mRNA of the 
invented method. 

Figure 7 shows a photographic image of a polaroid film exposed to chemiluminescent 
light derived from oligo-(dT) probe by the invented method, in which Gi-2 and GS protein 
30 mRNA in various cells were trapped onto the first polynucleotide-immobilized plastic plate. 

Figure 8 shows a photographic image of an X-ray film exposed to chemiluminescent 
light derived from the Gi-2 and Gs protein cDNA probes on the Southern membrane, in 
which PCR-amplified G protein specific DNA from 4 different cell lines was immobilized. 
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Figure 9 graphically depicts a calibration line for the human substance P receptor 
mRNA of the invented method. 

Figure 10 shows a gel showing the effect of various RNase inhibitors on mRNA 
preparations containing SDS and EDTA. 
5 Figure 11 shows a gel of five different G protein oligonucleotides amplified with the 

G*> and G 4 PCR primers and also provides Southern blots using each of the five G protein 
sequences as a probe. 

Figure 12 shows a gel of samples containing varying amounts of each of five different 
G protein oligonucleotides (as indicated by number of " + " symbols) amplified with the G 2 
10 and G 4 PCR primers and also provides Southern blots using each of the five G protein 

sequences as a probe. 

Figure 13 shows a gel of XZAP cDNA libraries from rate pituitary (P), kidney (K) 
and intestinatal (I) amplified with each of five different G protein primers, as indicated. 

Figure 14 shows a gel of 500 bp DNA from cDNAs of human IM9 and Jurkat cells 
15 amplified with G 2 and G 4 PCR primers. 

Figure 15 is a graph showing the relationship between YOYO-1 fluorescense and 
amount of immobilized oligonucleotide. 

Figure 16 is a graph showing the time course of YOYO-1 analysis. 

Figure 17 is a graph showing the dose response of YOY-O-1 concentration. 
20 Figure 18 is a bar graph showing the reproducibility of YOYO- 1 staining based on 

quantitation of immobilized oligonucleotides. 

Figure 19 exhibits a Polaroid picture of agarose gel electrophoresis of the PCR 
amplified ds-cDNA of jun oncogene from human peripheral blood leukocytes.' 

Figure 20 is a simplified block diagram of a computer system illustrating the overall 
25 design of this invention; 

Figures 21 A-21C show display screen representations of the main oligoprobe design 
station dialog windows of this invention; 

Figure 22 is a flow chart of the overall invention illustrating the program and the 
invention s sequence and structure: 
30 Figure 23 is a display screen representation of the Mitsuhashi probe selection 

diagram; 

Figure 24 is a display screen representation of the probeinfo and matchinfo window; 
Figure 25 is a display screen representation of the probesedit window; 
Figure 26A is a printout of the probesedit output file; 
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Figure 27 is a flow chart of the overall k_diff program of the Mismatch Model of this 
invention, including its sequence and structure; 

Figure 28 is a flow chart of the k_diff module of this invention; 

Figure 29 is a flow char: of the hashing module of this invention; 
5 Figure 30 is a flow chart of the tran module of this invention; 

Figure 31 is a flow chart of the let_dig module of this invention; 

Figure 32 is a flow chart of the update module of this invention; 

Figure 33 is a flow chart of the assembly module of this invention; 

Figure 34 is a flow chart of the seqload module of this invention; 
10 Figure 35 is a flow chart of the readl module of this invention; 

Figure 36 is a flow chart of the digjet module oi this invention; 

Figure 37 is a flow chart of the qjrolour module of this invention; 

Figure 38 is a flow chart of the hit_ext module of this invention; 

Figure 39 is a flow chart of the colour module of this invention; 
15 Fi&xre 40 is the first page of a printout of a sample file containing the output of the 

Mismatch Model program of this invention; 

Figure 41 is a flow chart of the H-Site Model, stage I, covering the creation of a 
preprocessed preparation file of this invention: 

Figure 42 is a flow chart of the H-Site Model, stage II, covering the preparation of 
20 the target sequence(s); 

Figure 43 is a flow chart of the H-Site Model, stage III. covering the calculation of 
MPSD data; 

Figure 44A is the first page of a printout ot a sample file containing output of 
Mismatch Model program; 

25 Figure 44B is the first pace of a printout of a sample file containing output of H-Site 

Model program; 

Figure 45 is a flow chart of the processing used to create the Mitsuhashi probe 
selection diagram (MPSD); 

Figure 46 is a flow chart of processing used to create the matchinfo window; 
30 Figure 47 is a printout of a sample target species file; 

Figure 48 is a printout of a sample preparation file. 

Detailed Description of the Invention 

Target mRNA 
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Previously, in order to genericaUy diagnose an infectious disease, procedures utilizing 
DNA were generaUy performed. Although, with DNA diagnosis, the existence of an external 
gene can be confirmed, it cannot be determined if the known gene is in an unexpressed 
carrier state or expressed in the diseased state. 
5 In living organisms, once a gene is activated, the DNA coding for that gene is 

transcribed into mRNA, and the mRNA is subsequently translated into protein. Thus, the 
activation of the external gene can be determined by measuring the levels of mRNA. 

Currently, to genetically diagnose a protein deficiency state, most experimenters have 
limited themselves to analysis at the DNA level. However, this method cannot be used to 
10 determine abnormalities of transcription or translation. If mRNA levels could be measured 

easily, the diseased state might be more precisely recognized. In a diseased state where 
mRNA cannot be detected it can be said that there is a defect in transcription, while if 
normal levels of mRNA are detected, there is probably an abnormality in translation. 

In addition, the amount of mRNA produced in the cell reflects the amount of 
15 expression of a gene, so it would be useful in the course of an abnormal condition, such as 

an illness or drug reaction, to examine mRNA levels. As discussed in greater detail below, 
measuring any of a large number of mRNAs can provide clinically useful information. 
Particular examples include mRNA for the p receptor, substance P, G protein and 
oncogenes, such as jun. 

20 For the aforementioned reasons, the ability to measure mRNA would be useful for 

diagnosing and recognizing the pathophysiology of various disease states. This method 
would be especially useful for diagnosis of hereditary diseases, cancer, and infectious 
diseases. 

A preferred target of this invention is mRNA of eukaryote cells which contain 
25 polyadenylic acid tail on their S'-terminal end. Once the specific nucleotide sequence of the 

mRNA is known, virtually any mRNA can be measured using the methods of the present 
invention. Examples of such mRNAs are human G protein, human P ? receptor, substance 
P, lymphocyte surface antigen, immunoglobulin, collagen, and adrenaline receptor. 
Tggt Samples 

30 Test samples containing target mRNA can be obtained from a variety of sources. 

Preferred sources include living samples or cell solutions which have been treated to 
inactivate ribonuclease (RNase) activity. Purification of mRNA for standard samples can 
be conducted according to the methods usually associated with Northern Blot procedures 
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or Dot Blot procedures method. Kits for the purification of mRNA are commercially 
available. 

A test sample can be prepared from a living sample by treatment with 10 mM of 
ethyJenediaminetetraacetic acid (EDTA) (pH 8.0). 0.2 M NaCl. 0.5% of sodium dodecyl 
5 sulfate (SDS), 500 Unit/ml of RNase inhibitor. K)mM of Vanadyl Ribonucleoside Complex 

and 200 ug/1 of Proteinase K (hereafter called Lysis buffer) to lyse the cells and inhibit 
RNase activity. After lysis of the cells, NaCl concentration is adjusted to 0.5M. 

In the methods of the present invention that involve probing samples containing 
RNA. or in any procedure in which it is desired to prevent the degradation of RNA, it is 
10 advantageous to inhibit the activity of any ribonucleases (RNases) which may be present. 

One Rnase inhibitor which can be used is Vanadyl Ribonucleoside Complex (VRQ 
(Bethesda Research Laboratories. Gaithersbure. MD). VRC has been reported to be useful 
during cell fractionation and in the preparation of RNA. and has been shown not to 
interfere with the phenol extraction or ethanol precipitation of RNA. In addition, VRC 
15 does not affect other cytoplasmic components of cells. Therefore. VRC is an ideal inhibitor 
of RNase in many experimental procedures. 

However, prior art procedures for inhibiting RNases with VRC taught that VRC 
should not be used in buffer systems containing EDTA or SDS. which are commonly used 
in the field of molecular biology. The reason for this prohibition was that it was believed 
that the Complex would dissociate in the presence of these buffers, leading to a loss of 
RNase inhibiting activity In fact, the BRL insert accompanying the VRC product 
recommends that a five- to ten-fold molar excess of EDTA be added to an RNA solution 
containing VRC in order "destroy" the VRC prior to ethanol extraction of the RNA from 
the solution. The apparent inability to use VRC together with common buffers that include 
25 EDTA and SDS thus presented a major impediment to the exploitation of VRC as an 

RNase inhibitor. 

We have discovered, however, that VRC is in fact an effective RNase inhibitor even 
in the presence of SDS and/or EDTA. Thus. VRC can be used in assays which use buffers 
including SDS or EDTA, where heretofore it was believed that VRC would not be effective 
in such systems. VRC is an effective RNase inhibitor at the concentrations of EDTA and 
SDS that are normally used when manipulating RNA or when performing a variety of other 
molecular biology techniques. For example, we have found that VRC effectively inhibits 
RNase in a buffer solution comprising approximate!) 1 mM EDTA. We also found VRC 
to be effective in .5% SDS solutions, and is believed to be effective in solutions ranging up 



20 



30 
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to 2% SDS, more preferably up to 1% SDS. VRC can. of course, also be used in other 
solutions including EDTA and SDS. 

We have found that VRC is a particularly potent inhibitor of RNases when used in 
combination with Proteinase K. Proteinase K is also available from BRL. As shown in the 
5 gel in Figure 10. mRNA prepared from U937 cells (human macrophage cell line) was 

protected from RNase degradation by a combination of VRC and Proteinase K. Lane 1 of 
the gel shows the mRNA from a cell preparation which included VRC, Proteinase K. and 
RNasin, while lane 2 represents the mRNA from a cell preparation that included VRC and 
Proteinase K. RNasin is available from Promega of Madison. WL The distinct band 10 in 

10 lanes 1 and 2 matches the band seen in lane 7, which contains pure clonal cDNA from an 

RNase-free preparation of U937 cells, thus showing that the mRNA in the preparations of 
lanes 1 and 2 did not experience substantial mRNA degradation. 

A comparison of lanes 1 and 2 with lanes 3-6 shows that VRC in combination with 
Proteinase K inhibits RNase activity in the above-mentioned mRNA preparation from U937 

15 cells to a far greater extent than either Protein K alone (lane 4), Proteinase K in 

combination with RNasin (lane 3), or a commercial RNase-inhibiting preparation sold under 
the name TastTrack" (available from In Vitrogen of San Diego. CA) (lane 5). None of 
lanes 3-5 exhibit the distinct band 10 representing undegraded mRNA. On the contrary, 
lane 4 (Proteinase K alone) and lane 5 (FastTrack) have the' same smear of degraded 

20 mRNA as band 20 in lane 6 (no inhibitors). The gel shown in Figure 10 also shows the 

effectiveness of VRC in a buffer solution of I mM EDTA and .5% SDS. which is the buffer 
used in the U937 cell preparations tested, since lane 2 (VRC and Proteinase K) shows less 
mRNA degradation than lane 4 (Proteinase K alone) or lane 3 (Proteinase K and RNasin). 
First Polynucleotide Probe 

25 In the practice of the present invention, a first polynucleotide probe is bound to an 

insoluble support. This first polynucleotide is used to trap the target mRNA to a solid 
support. Thus, the first polynucleotide probe comprises a nucleotide sequences which can 
hybridize with a sequence of the mRNA that does not occur in other RNA in the sample. 
Such a sequence is referred to herein as a unique region of the mRNA. Preferably, the first 

30 polynucleotide probe does not bind to portions of the target mRNA other than the unique 

region. 

A unique region can be identified by any of several methods. These include: (1) 
comparing the nucleotide sequence between several potential mRNAs and their target 
mRNAs; (2) getting a candidate sequence which can be compared with other nucleotide 
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sequences; and (3) computer searching for nucleotide sequences which do no: pair with the 
target mRNA. Information on various known mRNA sequences can be obtained from 
several genetic information databases. 

The first polynucleotide probe is preferably an oligodeoxyribonucleotide (DNA) 

5 rather than an oligoribonucleotide (RNA) since DNA is more stable than RNA. The 

number of nucleotides is not restricted; however, if an oligodeoxyonucleotide is used as the 
first probe, a preferred length is 17-40 nucleotides. If a cDNA of target RNA is available, 
a part or the entire cDNA can be used as the first nucleotide probe. The first nucleotide 
probe can also be easily manufactured by using a DNA synthesizer or by a commercial 

0 manufacturer. 

In this invention, the first nucleotide probe used is immobilized to an insoluble 
support. It may be chemically modified for immobilization, such as by adding an amino 
group onto the f'-terminal or 3*-terminal of the polynucleotide. 
Insoluble Support 

5 In the practice of the present invention, the first nucleotide probe is immobilized to 

an insoluble support. The insoluble support can be any of a variety of insoluble materials, 
such as a plastic plate, membrane filter, microtiter plate or any other insoluble material to 
which polynucleotides can be attached. The insoluble support is preferably made of a 
material which can immobilize the first nucleotide probe by covalent bonds. For example, 

0 a plastic plate, such as polystyrene, which has a carboxyi group or amino group on the 

surface thereof can be used. Such plates are available commercially. However, procedures 
for introducing these residues are also well known in the art. One commercially available 
plastic plate which has carboxyi group, on its surface is Sumilon microtiter plate MS-3796F 
made by Sumitomo Bakelite. A plastic plate having an amino group on the surface is the 

5 Sumilon microtiter plate MS-3696F. 

Immobilization of First Polynucleotide 

In the practice of the present invention, the first polynucleotide probe is immobilized 
onto the insoluble support. Various methods of immobilizing polynucleotides to the 
insoluble support are available, including covalent binding, ionic binding, and the physical 

0 absorbance method. However, the covalent binding method is preferred. Thus, in certain 

embodiments of the present invention, the polynucleotides are immobilized to microtiter 
plates which exhibit functional groups, such as carboxyi residues, amine residues, or hydroxyl 
residues on the surface thereof. 
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In a preferred procedure for immobilization of the polynucleotide to an insoluble 
support exhibiting a functional group, the 5*-terminal end of the polynucleotide is covalently 
linked to the functional group. Any of a variety of methods for covalent binding of 
polynucleotides to these functional groups can be used. Examples of preferred weii-known 
5 methods include the maleimide method and carbodiimide method. 

The maleimide method, illustrated in Figure I. involves the reaction between a 
substance containing a maleimide group and another material containing a sulfhydryl residue 
(SH). In order to attach the 5* end of a polynucleotide to an immobilized support using the 
maleimide method, the 5' end of the polynucleotide is reacted with a maleimide compound. 

10 A suitable maleimide compound is sulfosuccinimidyl-4-(N-maleimidomethyl) cyclohexane-1- 
carboxylate (sulfo-SMCC). 

The SH residue is provided on the support by a reaction between a support having 
an amino group and succinimidyl-S-acetylthioacetate (SATA), followed by deacetylation 
using hydroxylamine (NH ; OH). (Sulfo-SMCC and SATA are readily available from a 

15 variety of commercial sources, including the Pierce Company.) The resulting SH group on 

the support can then react with the maleimide group on the 5* end of the polynucleotide to 
form a polynucleotide-immobilized support. One problem we have experienced in 

the use of the maleimide method is that the SH group on the plate can react not only with 
an amino group at the 5' end of the polynucleotide, but can also react with primary amino 

20 groups on the purine bases, adenine and guanine. In order to assure that the 
polynucleotides are immobilized at their 5' ends, the amino groups on the purine bases can 
be protected by pairing the polynucleotide to a complementary polynucleotide prior to 
immobilization. After immobilization, the complementary polynucleotide can be removed 
through denaturation. such as through heating. 

25 The carbodiimide method is illustrated in Figure 2. This method involves a reaction 

between an amino group and a material containing a carboxyl residue using a carbodiimide 
compound. An example of a carbodiimide compound is l-ethyI-3-(3-dimethylaminopropyl) 
carbodiimide hydrochloride (hereafter called EDC). This reaction can be enhanced with 
N-hydroxysulfosuccinimide (hereafter called Sulfo-NHS). Both EDC and Sulfo-NHS are 

30 available from well known commercial sources, including the Pierce Companv. 

In the practice of a preferred carbodiimide method for attaching polynucleotides to 
a support, a support having a carboxyl residue attached is used. EDC is activated by 
reaction with Sulfo-NHS. This activated EDC is reacted with support containing surface- 



WO 93/15221 



PCT/US93/00977 



-15- 

bound carboxyl residues. This can then be reacted with polynucleotides having an amino 
group at their S'-terminal ends, resulting in a polynucleotide-immobilized support. 

In order to assure that the polynucleotides are immobilized at their 5* end, primary 
amino groups on adenyl, guanyl and cytosyl croups can be protected by hybridizing the 

5 nucleotide to a complementary polynucleotide prior to immobilization After 

immobilization, the complementary polynucleotide can be removed through denaturation, 
such as through heating 

We have found that non-specific binding of activated amino or carboxyl residues on 
insoluble supports can be effectively reduced or eliminated by treating the plates to which 

10 polynucleotides have been immobilized with a primary amine compound, preferably glycine. 

Second Polynucleotide Probe 

The second nucleotide probe contains sequences which can also hybridize to target 
mRNA that is trapped by the first bound nucleotide sequence. Preferably, the second 
polynucleotide probe hybridizes to a different portion of the target mRNA than does the 

15 first polynucleotide probe. In an especially preferred embodiment, the second nucleotide 

probe consists of nucleotides which are complementary to the polyadenylic acid tail of the 
target mRNA. The second nucleotide probe is also preferably DNA rather than RNA due 
to the higher stability of DNA. The number of nucleotides is not limited; however, a 
polynucleotide of 15-40 nucleotides is preferable. Thus a poly (dT) sequence of 15-40 

20 nucleotides is a preferred second polynucleotide probe. Such a second polynucleotide probe 

would serve as a universal second probe for virtually any target mRNA from eukaryotic 
cells. However, if a cDNA to the target RNA is available, a part or all of this cDNA can 
also function as a second polynucleotide probe. The second nucleotide probe can also be 
easily manufactured by a DNA synthesizer or commercial manufacturer. 

25 The second polynucleotide prot>e is preferably labeled in order to easily detect its 

presence. A variety of chemical substances are available which can label the target 
substance. For example, a variety of radionuclides can be used, such as the radioisotopes 
32 P, 35 S, 3 H. and 125 L Enzymes can also be used to label the second polynucleotide probes. 
Suitable enzymes include alkaline phosphatase, lucifenise. and peroxidase. Other labels 

30 might be chemical compounds such as biotin. avidin. streptavidin. and digoxigenin. Labels 

which provide a colorimetric indication or a radionuclide are preferred. In a preferred 
embodiment, biotin is attached to the nucleotide probe, followed by an avidin. such as 
streptavidin. conjugated to an enzyme such as alkaline phosphatase. The presence of the 
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enzyme can be detected by various substrates, such as ATTOPHOS which provides a 
flourescent marker that can be detected through fluorimetry. 

In certain embodiments, part or all of the labeling compound is substituted with part 
or ail of the second nucleotide probe. Techniques tor substituting portions of 
5 polynucleotides are well known for a variety of labels and kits are commercially available 

for many of these. For example, one can use a DNA 3'-end Labeling Kit or oligo- 
deoxythymidylic acid to label the probe. 

Nucleic acid can also be "labeled" by staining with a nucleic acid stain. Thus, where 
a relatively large amount of nucleic acid is present, ethidium bromid (EtBr) can be used to 
10 identify the presence of nucleic acid. This is useful within the context of the present 

invention, for example, if mRNA has hybridized to a short immobilized anti-sense probe 
then substantially more nucleic acid will be immobilized than if no mRNa had hybridized. 
However, more sensitive stains are more preferable. Such more sensitive stains include the 
various cyanine nucleic acid stains, such as POPO, BOBO. YOYO and TOTO available 
15 from Moecular Probes (California). These stains are described, e.g., in Science, 257:885 

(1992). Particularly preferred stains for use in the context of the present invention include 
the shorter wavelength forms, TOTO-1 and YOYO-L still more preferably YOYO-L As 
little as four picograms of stained DNA can be detected by visible fluorescence upon 
stimulation with a transilluminator or hand-held UV lamp. Thus, these stains provide a 
20 particularly easy and sensitive method of identifying the presence of nucleic acids. 

We tested the ability of YOYO-1 to stain oligonucleotides immobilized to wells as 
discussed herein. We first immobilized various known amounts of oligonucleotides to wells 
and washed to remove non-immobilized oligonucleotides. We then added 1:1000 dilution 
of YOYO-1 in water. We then used a fluorimeter directly without washing. The relation 
25 between pmoles of oligonucleotide is shown in Figure 15 with open circles. We also 

incubated the TOTO-1 stained immobilized oligonucleotides for ten minutes, washed and 
added water, followed by use of the fluorimeter. The washed results are shown in dark 
circles in Ficure 15. It can be seen that washing does not sienicantly effect the amount of 
staining and that the amount of staining is related to the amount of oligonucleotide. 
30 We also tested the time course of YOYO-1 staining over the course of one hour. 

We included a control in which no oligonucleotide was immobilized to the plate. The 
oligonucleotide-immobilized plate is shown in open circles in Figure 16 and the control in 
dark circles in that Figure. It can be seen from Figure 16 that after an initial spike, 
relatively constant staining is found. 
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We further tested the dose response of constant amounts of oligonucleotides with 
oligonucleotides ..mobilized in wells (open circles in Figure 17) and control wells with no 
oligonucleotide (closed circles in Figure 17). The difference between the immobdized and 
non-immobilized is shown in Figure 17 as triangles. It can be seen that a sharp .ncre.se in 
5 fluorescence occurs between 10" 4 and lO' 3 dilution. 

We also tested the reproducability of staining of both oligonucleotide-immobilized 
wells and non-immobilized wells. We repeated the experiment five times and graphically 
depicted the data for both oligonucleotide immobilized ( + ) and non-immobilized (-) wells 
in Figure 18 It can be seen that the data was substantially similar for each experiment. 
10 Thus, the data depicted in Figures 15-18 shows that use of YOYO-1 as a staining 

agent provides a reliable indication of the amount of oligonucleotide present. 

/V ppHratinn nf Sample 

In the preferred methods of the present invention, a sample containing mRNA is 
applied to the insoluble support to which the firs, polynucleotide probe has been bound. 
15 Any tarcet mRNA found in the sample is allowed to hybridize with the first polynucleotide 
probe This can be accomplished by incubation at temperature dependent on a variety of 
factors, as is well known to those with ordinary skill in the art. These factors include the 
length of complementary nucleotide sequences, the ratio of Guanine or Cytosine (GC 
contents) within the entire base content in the complementary nucleotide sequences, NaCl 
20 concentration in buffer, the number of bases which mismatch in the complementary 
nucleotide sequence, and the type of nucleotide. In a preferred form of this invention, the 
following equation can be used to calculate the preferred incubation temperature: 
T inc , 16.6 x log (M) + 0.41 (GC) + 81.5 - 675/n - 15 (©C). 
In the equation shown above. M is the NaCl concentration (M) in solution. GC 
25 represents GC contents (%). n represents the length of nucleotide sequences (the number 

of nucleotides). 

The incubation time can also be determined according to methods described in the 

Molecular Cloning manual. 

The time "for incubation is preferably from 1 hour to overnight, and the sample 
30 should preferably be uently swung during incubation. Incubation is preferably performed 

in an appropriate buffer solution. The same buffer used to hybridize RNA and DNA in the 
Northern Blot or the Dot Blot Methods can be used. The buffer is preferably prepared in 
a way so as no. to contaminate it with RNase. If any RNase contamination is present, the 
activation thereof should be controlled to be as low ; ,s possible RNase-free buffers are 
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commercially available, such as within mRNA purification kits or lysis buffer as shown 
before. 

In order to eliminate RNase activity from water used in the methods of this 
invention, the water is preferably treated with Diethylpyrocarbonate (DEPC). The preferred 
5 DEPC treatment involves addition of 0.1% DEPC to the water, followed by storage 

overnight at 37°C and sterilization in an autoclave. 
Washing Procedures 

Preferably, after incubation, unbound components of the sample are washed from 
the insoluble support. Appropriate solutions for washing include the lysis buffer as shown 
10 before or the buffer included in the mRNA purification kit. 

A pplication of Second Polynucleotide Probe 

The second polynucleotide probe is applied to the insoluble support after any 
washing of unbound components of the sample, under conditions as described above in 
connection with the hybridization of mRNA in the sample with the first nucleotide probe. 
15 This is preferably followed by washing, as described above. 

Quantitation 

The amount of mRNA in the sample is quantitated by measuring the amount of 
second polynucleotide probe bound to the insoluble support after application thereto in 
accordance with the present invention. In this regard, a physical or chemical quantity or 

20 activity of the label on the second polynucleotide can be measured by several techniques, 

including optical density, emitted-light intensity, or radiation. The label itself can provide 
this indication or can require other compounds which bind thereto or catalyze the label. 
Other mechanisms for detecting label include the use of compounds that can chemically 
react against the label, detection of colored label, detection of light emission, detection of 

25 radiation or of catalytic ability. 

In a preferred measurement technique, the label on the second polynucleotide is 
biotin. The presence of this label can be detected by reaction with the enzymes Peroxidase 
or Alkaline phosphatase. These enzymes can be specifically directed to biotin by 
conjugation with avidin or streptavidin. The presence of the enzymes is then detected by 

30 addition of an appropriate substrate to provide a detectable color-developing or light- 

emitting reaction. The alkaline phosphatase labeled streptavidin can be readily obtained 
from the commercial market. A preferred substrate for alkaline phosphatase is adamantyl- 
1.2-dioxetane phosphate (AMPPD). Upon reaction with the alkaline phosphatase. AMPPD 
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will emit light at h wavelength of 447 nm. This light can be detected in accordance with 
techniques known in the art. 

In the reaction between alkaline phosphatase and AMPPD, an enhancer such as 5-N- 
tetradecanoyl-amino-fluorescein can be added. 5-N-tetradecanoyl-amino-fluorescein has the 
5 ability to convert light of 477 nm wavelength to 530 nm wavelength which is more readilv 

detectable. 

Other preferred labels include an antigen, such as digoxigenin or an antibody. An 
antigen can be detected by its ability to bind to an antibody directed thereto. Such 
antibodies, or an antibody directly used to label the second polynucleotide probe, can be 
0 detected by their ability to bind protein A. The antibody itself can be labeled directly with 

a radionuclide, such as 12:> 1, or can be labeled by binding thereto of protein A labeled with 
the radionuclide. The radionuclide can then be detected in accordance with techniques well 
known in the art. 

Well-known techniques for the detection of label include the detection of label which 
5 appears by a color-developing reaction by a spectrophotometer. Other such techniques 

include the detection of a light emitting reaction using X-ray film or an instant camera, and 

the detection of raditation using a radiation counter. 

The emission reactions are recorded by X-ray film or instant camera film in the dark 

room. The X-ray film which is exposed by emission reactions are recorded as a blot, so this 
D shading of the blot can be measured by a densitometer. If one uses an instant camera such 

as a Polaroid, the picture is read by a scanner to decide the location of the blot on the 

computer, measuring the shading of the blot using graphic analysis software. 

To measure the amount of target mRNA in an unknown sample solution, an 

arithmetic plot for showing the relationship between a standard solution of target mRNA 
5 and relative concentration in a blot can be used. The standard solution for target mRNA 

need not be purified target mRNA. Thus, the standards can also be a mixture in which the 

amount of target mRNA is known. 

Schematic Example 

A schematic representation of an exemplary method in accordance with the present 
0 invention is shown in Figure 3. The target mRNA I is shown having a polyadenylic acid tail 

2. The first nucleotide probe 3 is immobilized to an insoluble support 4, then mixed the 
target mRNA 1. The target mRNA 1 hybridizes with the first nucleotide probe 3. thereby 
immobilizing the target mRNA to the insoluble support. Then, a second nucleotide probe 
5 which bears a label 6 hybridizes with a portion of the target mRNA other than the portion 
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hybridized with the first nucleotide probe. As shown, the second nucleotide probe 5 
comprises poly (dT) and hybridizes with the polyadenylic acid tail 2 of the target mRNA 1. 
The resulting complex 7 can then be identified by detecting the presence of the label 6. 
Time-Savings of Invention 
5 Advantageously, the present invention provides a method of measuring mRNA which 

can be performed in far less time than previously known methods. Table 1 hereinbelow 
compares the time needed to measure DNA using the method of this invention and a 
Northern blot method. 
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10 



15 



20 



25 



30 



35 



40 



TABLE 1 

Comparison For Measurement Of mRNA 



Present Invention 



Time 
(hours) 



Northern Blot Method 



Time 

(hours) 



- Lysis of Cells 
-Hybridization of the 

first nucleotide probe 

to mRNA 



1 



» hybridization of the 
second nucleotide probe 
to mRNA 

•washing and blocking 

* detection 



1 

0.2 



* Purification of mRNA 
►Separation of mRNA 

by electrophoresis 

with Agarose Gel 
►denaturing of mRNA 

moving the mRNA to 

membrane 

"immobilization of mRNA 
to membrane 

► hybridization of the mRNA 
and probe 

" washing and blocking 

> detection 



48 



1 
2 

12 

4 
1 



Total 



about 5 hours 



about 70 hours 



As shown in Table L the methods of this invention generally take 5 hours or less to 
complete, while the Northern blot method generally requires about 70 hours. These times 
do not include the labeling time (about 12 hours) of the nucleotide probes nor the time for 
immobilization of the first nucleotide probe because these constructs can be prepared prior 
to the measurement. 

Particular examples of target mRNA that can be measured to obtain clinically useful 
information are provided below. The examples included herein are intended to illustrate 
certain preferred embodiments of the present invention. These examples are not intended 
to limit the invention in any manner. 
B Receptor 

A potential target for using mRNA measurement to determine a diseased state 
utilizes the mRNA for the p receptor. The p receptor is a protein located in human nerve 
tissue. In particular, abnormalities in the (J : receptor has been found to be closely 
correlated with asthma. Thus, measuring the mRNA for p 2 receptor can be used to 
determine the pathophysiology of asthma patients, and could also be used to assess the 
effectiveness of anti-asthma acents. 
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Measurement of Human p 2 Receptor mRNA 
(Preparation of Calibration Line) 

(1) Preparation ot the First Nucleotide Probe 

5 An oligodeoxyribonucleotide probe containing the nucleotide sequence shown below 

with an amino group attached to its 5' end was synthesized by Genosys, Inc. (TX). In this 
sequence. A, C, G and T indicate adenine, cytosine, guanine, and thymine, respectively. 
Based on computer analysis (DNASIS, Hitachi Software Engineering America, Ltd., CA), 
this sequence was determined to be homologous to the human P 2 adrenergic receptor 
10 mRNA. 

5'-NH 2 -ATG CTG GCC GTG ACG CAC AGC A-3" (Seq ID NO: 6) 

(2) Immobilization of the First Nucleotide Probe Onto Insoluble Supports 

Both EDC and sulfo-NHS (Pierce, IL) were dissolved in DEPC-treated water at 
concentrations of 20 mM and lOmM. respectively. EDC/Sulfo-NHS solution was then 
15 prepared by mixing equal volumes of both EDC and sulfo-NHS. The first nucleotide probe 

was dissolved in DEPC-treated water at a concentration of ljig/iil and then mixed with the 
EDC/Sulfo-NHS solution in the ration 1:25 (VohVol). 

50 \il of this probe solution was added to each well of a microtiter plate (MS-3796F, 
Sumitomo Bakelite, JAPAN), which were known to express carboxyl groups on their surface. 
20 After incubation at room temperature overnight, the reaction solution was removed by 

aspiration. 

(3) Preparation of mRNA 

DNA plasmid pUCl8 containing the human adrenergic receptor cDNA (pTF\ 
American Type Culture Collection) was digested with Pstl and Eco RL and mixed with the 

25 pBC phagemid (Stratagene. CA) which had been predigested with the same restriction 
enzymes. The resultant plasmid was linearized by PstL extracted once with phenol- 
chloroform, precipitated in ethanol. and used to transcribe the human p n receptor mRNA 
using the Riboprobe system and T7 RNA polymerase (Promega, WI). The transcribed 
human (J 2 receptor mRNA was treated with RNase-free DNase to digest the remaining 

30 template cDNA, then purified with one cycle of phenol extraction and ethanol precipitation. 

The resultant mRNA showed up as a single band on agarose (FMC byproducts. ME) gel 
electrophoresis. Messenger RNA concentration was measured by a spectrophotometer at 
260 and 280 nm against standard mRNA solutions prepared at concentrations of 25 pg. 250 
pg, 2.5 ng. 25 ng and 250 ng. respectively. 
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(4) Hybridization of the First Nucleotide Probe to mRNA 

In order to remove RNase from the well, the firsi nucleotide probe-immobilized 
wells were treated with 250 nl of lysis buffer containing 0.5 M NaCl at 45oC for one hour. 
The buffer was removed from individual weiis by aspiration, and 5u \l[ ol lysis buffer 
5 containing varying concentrations of standard human p 2 receptor mRNA was added into 

each well. These solutions were incubated at 59° C for one hour to allow hybridization. 

(5) Preparation of the Second Nucleotide Probe 

The second nucleotide probe was a 15-mer oligo(dT) labeled with biotin at the 5' 

end. 

0 (6) Hybridization of the Labeled Second Nucleotide Probe to the Polv(A) tail of mRNA 

After hybridization of the first nucleotide probe and mRNA described in step (4), 
the hybridization solution was removed by aspiration and the individual wells were washed 
with 250 }il of lysis buffer. Fifty ni of lysis buffer containing 0.5 M NaCl and 1 >il of the 
biotinylated second probe solution (35 pmol/pl) was added into each well, and incubated 

5 at room temperature for an additional one hour. 

(7) Measurem ent of Chemical Activities of the Labeled Second Nucleotide Probe 
Following the hybridization described in step (6), the hybridization solution was 

removed by aspiration and wells were washed once with 250 \il of lysis buffer. A blocking 
buffer consisting of (0.05%(w/v) Tween 20, 500 mM NaCl. 100 mM Tris-HCl, pH 7.5) was 

0 added into each well and incubated at room temperature for five minutes to reduced 

nonspecific binding. Reaction solutions were then removed by aspiration, 50 \il of blocking 
solution containing 1/5000 vol. of alkaline phosphatase-conjugated streptavidin solution 
(Clontech, CA) was added into each well, and incubated at room temperature for 30 
minutes. Each well was then washed with 250 ^1 of washing buffer containing 500 mM 

5 NaCl, 100 mM Tris-HCl. pH 7.5. and 50 |il of Lumiphos 530 solution (Boehringer- 

Mannheim. IN) whole plates were placed on X-ray film (XAR5. Kodak, NY) in specific 
holders (MYC. Kodak, NY). X-ray film was exposed to a generated chemiluminescent light 
at room temperature for one hour, and then developed using a film processor (QX-400, 
Konica. Japan). Relative intensities of individual wells on X-ray films were quantified by 

0 densitometry (model 620. Bio-Rad, CA). 

(8) Results 

The results of Example 1 are shown in Figure 8. The inset (A) shows results of an 
X-ray film exposed to chemiluminescent light during the detection of mRNAs by the 
invented method; (B) shows a densitographic representation of the same imace as in (A): 
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(C) exemplifies a calibration line of the mRNA in which the X-axis represents 
concentrations of mRNA and the Y-axis represents the relative intensity of 
chemiluminescence in an individual well . Wells 1. 2. 3. 4 and 5 of Figure 8 correspond to 
the amounts of tested mRNA at 25 pg, 250 pg, 2.5 ng, 25 ng and 250 ng. respectively. As 
5 shown in (A), mRNAs were detected in amounts over 25 pg by the invented method. The 

calibration line shown in (C) shows the high linearity of signal when increasing 
concentrations from 25 pg to 250 ng of mRNA. 

COMPARISON 1 
Quantification of mRNA hv Northern Blot 

10 (1) Agarose Gel Electrophoresis of mRNA 

Ten mM of sodium phosphate buffer, pH 7.0 containing 0.1 % DEPC was heated 
to 37oC overnight and autoclaved. A \% agarose solution was prepared using this 10 mM 
DEPC-treated phosphate buffer. The solution was chilled to ca. 50°C. and iodine acetate 
was added to 0.2% (W/V). A gel slab was formed (6.5 cm x 10 cm) with combs and placed 

15 at room temperature for one hour to allow hardening. Five \lI of DEPC-treated water 

containing 250ng, 25 ng, 2.5 ng, 250 pg or 25 pg of mRNA described in example 1.(3), were 
mixed with 7 \il of loading buffer (125 \x\ DMSO. 42 ^1 40% glyoxal, 3 \xl 1 M sodium 
phosphate buffer, pH 7.0), and incubated at 50oC for one hour. Samples were applied to 
the gel, and electrophoresed in a solution of DEPC-treated 10 mM phosphate buffer at 

20 constant 6V/cm for 15 minutes without circulation. An additional one hour of 

electrophoresis was performed with circulation. 

Following electrophoresis, gels were immersed in 50mM sodium hydroxide containing 
0.5 ng/ml ethidium bromide (Sigma, MO) for 20 minutes, and 0.1 M Tris-HCl buffer, pH 
7.5 for 40 minutes. Gels were then transferred into fresh 0.1 M Tris-HCl buffer and 

25 photographed under UV light. 
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(2) Northern Transfer of mRNA Onto Membranes 

Filter paper (3MM Chr, Whatman, NJ) and a nylon membrane (Magnagraph, 
Micron Separations, MA) were immersed in sterile water for 20 minutes and then placed 
in lOx SSPE (15 M Nad. 115 mM sodium phosphate, 11.5 mM EDTA in water, pH 7.4 
5 adjusted with 10 N NaOH). A blotting sponge (Stratagen, CA) was also immersed in lOx 

SSPE. The 3MM filter, nylon membrane, agarose gel, and sponge were assembled from the 
bottom in a Posi-blot Pressure blotter (Staratagen, CA), to cause the mRNA to transferred 
onto the nylon membrane by positive pressure (835 mmHg for 45 minutes). The nylon 
membrane was dried, and the mRNA was permanently immobilized onto the nylon 
10 membrane by UV light cross linkage at 120 mjoules (Stratalinker 1800. Stratagen. CA). 

(3) Preparation of Biotinvlated cDNA Probe 

DNA probes were prepared from p 2 receptor cDNA by using a DNA labeling kit 
(Boehringer-Mannheim. IN) and biotin-labeled deoxyuridylic acid (Clontech. CA). In brief, 
a (J 2 receptor cDNA -containing plasmid, pTF was digested with Pstl and £a>RI. The 

15 resultant plasmid and cDNA were separated by agarose gel electrophoresis. A band of p 2 

receptor cDNA was cut out and removed from the gel using the Gene clean kit (Bio 101, 
CA). 0.1 fig of above cDNA dissolved in 5 \il of water was incubated at 95©C for 10 
minutes, and then incubated with 2 \il of hexanucleotide mixture, 2 \l\ of dNTP and 1 \xl of 
Klenow enzyme (DNA labeling kit, Boehringer-Mannheim, IN) at 37oC for one hour. 

20 Following incubation, the reaction was terminated by addition of 2 \i\ of 0.2 M EDTA, pH 

8.0, 2.5 \il of 4.0 M LiCl and 75 \il of 100% ethanol, and cooled at -70oC for 3.0 minutes to 
precipitate cDNA. After centrifugation, the pellet whs washed once with 70% ethanol and 
then dissolved in 50 \il of water. 

(4) Hybridization of Biotinvlated cDNA Probe to mRNA 

25 The mRNA containing membranes were incubated in hybridization buffer containing 

0.5 M NaCl (ECL, Amersham. IL) at 45°C for one hour. Biotinylated cDNA probe 
described above were first heat denatured at 95©C for 10 minutes, immediately cooled on 
ice. and then added to the hybridization buffer. Hybridization was continued at 40oC 
overnight. 

30 (5) Chemical Activity Measurements of Labeled cDNA Probe 

After hybridization, membranes were washed twice in 0.5x SSPE solution containing 
36% (W/V) urea. 0.4% (W/V) SDS and 0.5 M NaCI at 65<>C for 15 minutes each, followed 
by washing twice in 2x SSPE with 0.5 M NaCl at room temperature for five minutes each. 
Membranes were then incubated with 100 mM Tris-HCl buffer. pH 7.0. containins 2% 
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blocking reagent (Boehringer-Mannheim. IN) and 150 mM NaCl at room temperature for 
three hours. 1 /5000 vol. of alkaline phosphatase-conjugated streptavidin (Clontech. CA) was 
added to the above solution and incubated at room temperature 30 minutes. Membranes 
were washed twice with 100 mM "Iris-HQ, pH 7.0, containing 150 mM NaCl at room 
5 temperature for 30 minutes each followed by one washing with alkaline phosphate buffer 

(100 mM NaCl. 5 mM MgCU in Tris-HQ buffer, pH 9.5) at room temperature for 2 
minutes. Membranes were then soaked in Lumi-Phos 530 (Boehringer-Mannheim, IN), and 
quickly placed between two sheets of transparent plastic film. Chemiluminescent membrane 
signals were exposed to X-ray film (XAR5, Kodak, NY) in specific holders (MYC, Kodak, 
10 NY) at room temperature for one hour. 

(6) Results 

The Northern blot results of Comparison 1 are shown in Figure 8(1), in which the 
chemiluminescent signals of mRNA were exposed to an X-ray film. Lanes 1, 2, 3. 4. and 5 
of Figure 8(1) correspond to 25 pg. 250 pg, 2.5 ng. 25 ng and 250ng of applied mRNA, 

15 respectively. Because this cDNA probe is approximately 100 fold larger than the 

oligonucleotide probe used in Example 1, the amount of incorporated biotin is also 100 fold 
more than the oligonucleotide. This result indicates that the cDNA probes used in 
Comparison I are approximately 100 fold more sensitive than the oligonucleotide probe. 
However, as seen in Figure 8(1), a conventional Northern blot can only detect samples 

20 containing greater than 25 ng of mRNA. This is far less sensitive than the invented method 

as exemplified by Figure 8(1). Also, the invention as described by Figure 8(1) used less 
sensitive oligonucleotide probes when compared to the cDNA probes used on the Northern 
Blot. 

EXAMPLE 2 

25 Detection of Human Receptor mRNA From Unknown 

Samples And Evaluation of Reproducibility 
(1) Preparation of the First Nucleotide Probe and Its Immobilization Onto the 
Insoluble Support 

The first nucleotide probe was prepared and immobilized onto an insoluble support 
30 as described in Example 1. 
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(2) Preparation of Standard mRNA 

mRNA was prepared as described in Example 1 

(3) Preparation of Cell Lvsates 

Human HL60 promyelocyte leukemia ecll:> were collected by centrifugation. wasned 
with PBS (8 g NaCl. 0.2 g KCL 1.44 g Na 2 HP 4 and 0.24 g KH : P0 4 in 1 liter of water. pH 
7.4). Cells were then resuspended in lysis buffer at a concentration of 1 x 10 7 cells/ml. The 
solution was repeatedly passed through a steriJe plastic syringe fitted with a 21 cause needle 
and incubated in a slow-shaking waterbath at 45°C for one hour. 

(4) Hybridization and Chemiluminescent Detection 

As described in Example 1, hybridization of the first nucleotide probe, preparation 
of the second nucleotide probe, hybridization of the labeled second nucleotide probe to the 
poly(A) tail of mRNA. and chemical activity measurement of the labeled second nucleotide 
probe was carried out. 

(5) Preparation of Calibration Line 

As described in Example 1, a calibration line of known p 2 receptor mRNA 
concentrations was created (Figure 9). The calibration line expresses a linear co- 
relationship between mRNA concentrations mapped to a logarithmic scale and the relative 
chemiluminescent intensity read from X-ray film. 

(6) Evaluation of Reproducibility 

Concentrations of p ; receptor specific mRNA from HL60 cell lysates were 
determined three times. In Table 2, each concentration of p : receptor specific mRNA was 
determined from the calibration line described in section (5) above. The mean and standard 
error (S.E.) of data were calculated. Cell numbers were counted by a phase hemacytometer 
(American Scientific Products. IL). As shown in Table 2, three separate assays exhibited 
similar mRNA values at a range of 0.96 to 1.12 fg/ceU. This suggests that the invented 
method provides reproducible method of quantifying specific mRNA levels. 
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TABLE 2 

Concentrations of Receptor-Specific mRNA in HL60 Ceils 
No. of assay mRNA Concentration [tg/cell] 

5 



1 i.Ol 

2 0.96 

3 1.12 

10 

mean ± S.E. 1.03 ± 0.04 

G Proteins 

Cell surface receptors for hormones and neurotransmitters are known to be coupled 
to intracellular heterotrimeric GTP-binding proteins (G proteins) composed of a, p and t 

15 subunits. Once receptors are activated by specific ligands. receptor-coupled G proteins 

transduce signals to intracellular secondary effector systems, such as adenylyl cyclase, 
phospholipase C and ion channels. 

Recently, the molecular mechanisms of the interaction between receptors and G 
proteins have become less clear, because extensive molecular qloning has identified more 

20 subclasses of G protein a subunits than there are known cellular functions of G proteins. 

Furthermore, it has been found that a single receptor is also capable of coupling to multiple 
G proteins. For example, histamine H 2 receptors in human HL-60 cells were coupled to two 
different G proteins, which activated adenylyl cyclase and phospholipase C, respectively. It 
has also been reported that GABA-B receptors were capable of coupling to G 0 . G 0 * and G^ 

25 j proteins, but not G w proteins. Moreover. CHO cells have been shown to express multiple 

G proteins, which activate the same phospholipase C pathways but are coupled selectively 
to different receptors. Some receptors may alter their G-protein subclass association 
depending on the state of differentiation of the cell or on the particular phase of the cell 
cycle.One example would be the mRNA coding for G protein. G protein binds to and 

30 controls the function of various receptors in the cell membrane. It catalyzes the intracellular 

reactions of cAMP and calcium ion free reaction systems. These molecules are used in 
various cellular transmitter systems. 

G proteins are believed to be involved in causing various disease states. For 
example, a genetic deficiency of G s proteins is the molecular basis of hereditary 
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osteodystrophy. Pituitary tumors in acromegalic patients rv.ve been shown to contain 
mutant Gs proteins. G proteins are also involved in invasive and metastatic melanoma cells. 
Rat models of streptozotocin-induced experimental diabetes suggest that the levels of 
rnRNA for various subclasses of Ga proteins aie significantly altered from normai control 
5 rats. Furthermore, cellular functions of pertussis toxin-sensitive G proteins were shown to 

be significantly impaired in atherosclerotic porcine coronary arteries, while G protein 
function in leukocytes of patients with mania was hyperfunctional. Thus, those diseases 
which are known to be related to a malfunction of the G protein include atherosclerosis, 
bipolar disorder, diabetes mellitus, melanoma, and pituitary tumors. 

10 Currently available immunological detection methods (Western blots) and mRNA 

detection methods (Northern blots) are not sensitive and require a lot of cellular material, 
making it difficult to study the role of G proteins in disease. It is not easy to ascertain the 
quantity of expressed G protein. If the mRNA for the G protein was measured, it could be 
used to recognize the pathophysiology of patients having one of the foregoing ailments. The 

15 effects of agents designed to arrest or control these conditions could also be measured using 
this technique. G protein consists of a, (J and y subunits. The a subunit is known to have 
several allelic variations. Thus, by measuring each variation of G protein mRNA, the 
amount of expressed protein for each variation can be estimated. 

EXAMPLE 3 

20 Quantification of Human G Protein a 

Suhunit-Specif'ic mRNA 

(1) First Nucleotide Probe Preparation 

The first nucleotide probes used against Gs protein-specific mRNA and Gi-2 protein- 
specific mRNA were oligodeoxyribonucleotides. These nucleotide sequences, shown below, 
25 contain an amino group on their 5* end and were synthesized by Genosys (Texas). 

The probe sequence used for quantifying Gs mRNA is shown below (Seq ID No. 1): 

5'-NH r TTC ATC CTC CCA CAG AGC CTT G-3' 

30 The sequence of probe used for quantifying G i-2 mRNA is shown below (Seq ID No. 

2): 



5"-NHn-ATG GTC AGC CCA GAG CCT CCG G-3" 
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(2) Immobilization of the Second Nucleotide Probe Onto Insoluble Supports 

The first nucleotide probes were dissolved in DEPC-treated water at l^g/^L Probes 
were then mixed with the EDC/Sulfo-NHS solution described in Example 1 Paragraph (2) 
in a ratio of 1:25 (Vol:Vol). A 50 \i\ sample was added to each well of a Sumilon microtiter 
5 plate (MS-3796F, Sumitomo BAkelite, JAPAN) and incubated overnight at room 

temperature. Following incubation, the reaction solution was removed with an aspirator. 

(3) Preparation of Standard mRNA 

A plasmid vector pGEM2 containing rat Gs protein cDNA (provided by Dr. R. R. 
Reed, Johns Hopkins University, U.S A.) was linearized by digestion with Nhe] (Promega, 
10 WI). Gs protein-specific mRNA was prepared by the method described in Example 1 

Paragraph (3) using T7 RNA polymerase (Promega. WI). The resultant mRNA 
concentrations were determined by spectrophotometry and subsequently diluted with DEPC- 
treated water as described in Example 1 Paragraph (3). 

(4) Preparation of Cell Lvsates 

15 Human B lymphocyte-derived IM9 cells were centrifuged. washed with PBS, and 

resuspended in lysis buffer at 1 x 10 7 cells/ml as described in Example 2. Cells were then 
lysed by repeated passage through a sterile plastic syringe fitted with a 21 gauge needle, and 
then incubated in a slow-shaking waterbath at 45°C for one hour. 

(5) Hybridization of the First Nucleotide Probe to mRNA .* 

20 In order to remove RNase contamination in the microtiter well, 250 |xi of lysis 

buffer was added to each well at 45oC for one hour as described in Example 1. The buffer 
was aspirated, and 50 ^1 of lysis buffer containing 0.5 M NaCl and various concentrations 
of human Gs protein mRNA was added into each well. Each sample was incubated at 51°C 
for one hour to allow hybridization as described in Example 1. 



25 



WO 93/15221 



PCT/US93/00977 



-31- 

< 6 ) Hybridization of the L abeled Second Nucleotide Probe to a Polv(A) Tail of mRNA 1 
and Measurement of Chem ical Activities of the Labeled Second Nucleotide Pmhe 

After hybridization of the first nucleotide probe, the hybridization solution was 
aspirated off and individual wells were washed once with 250 \il of lysis buffer. Fifty \il of 
5 lysis buffer containing 0.5 M NaCl and 1 \il of the biotinylated second nucleotide probe 

solution (35 pmol/^1) was added into each well, and incubated at room temperature for one 
hour to allow hybridization. 

Chemical activities of the labeled second nucleotide probe was measured as 
described in Example 1. 

10 (7) Preparation of Calibration Line 

The calibration line of standard Gs protein-specific mRNA is shown in Fieure 10. 
The calibration line expresses a relationship between the common logarithm of mRNA 
concentrations and the relative intensity of chemiluminescent signals exposed to X-ray film. 
(8) Evaluation of Reproducibility 

15 Concentrations of Gs protein and Gi-2 protein specific mRNAs from IM9 cell iysates 

were measured six times to test reproducibility of the present invention. The methods used 
for Gi-2 were similar to that described above for Gs, with the incubation temperature used 
for hybridization of Gi-2 mRNA being 59°C Table 3 lists the mRNA concentrations as 
obtained from the calibration line, and, the mean, and the standard error. Cell numbers 

20 from each sample were counted with a phase hemacytometer (American Scientific Products, 

IL). As listed in Table 3. all six assays reported similar mRNA concentrations for both Gs 
and Gi-2 mRNA. This suggests that the current invention provides a reliable, reproducible 
means of measuring specific mRNA concentrations. 
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TABLE 3 



15 



No. of Assay 


mRNA Concentration [fg/ceil] 




Gs 


Gi-2 




1 


24 8 


4 6 


2 


10.6 


8.8 


3 


18.8 


7.2 


4 


35.5 


6.8 


5 


13.6 


11.5 


6 


33.0 


13.7 


mean ± S.E. 


22.7 ± 3.8 


8.77 ± 1.24 



EXAMPLE 4 

Simultaneous Quantification of Various mRNAs 
Derived From Cell Lines hv the Invented Method 
20 in Conjunction with Reverse PCR Method 

(1) Simultaneous Detection of Gs and Gi-2 Protein-Specific mRNAs in Various Cell 
Lvsates bv the Invented Method 

Using the method described in Example 3. specific mRNAs for both Gs and Gi-2 
protein from various cell lvsates, such as Jurkat human T-cells. IM9 human B cell. U937 
25 human monocytes and HL60 human granulocytic cells were analyzed. Chemical activities 

of the labeled second probe were exposed to Polaroid film (667, Polaroid. MA) in a specific 
holder (model 901. Tropix. CA) at room temperature for 5 seconds, instead of exposing to 
X-ray films as described in Example 3. 

(2) Simultaneous Detection of Gs and Gi-2 Protein-Specific mRNAs bv the 
30 Conventional Reverse PCR Method 

(a) cDNA Preparation and Amplification hv PCR 

mRNA was purified using the FastTrack kit (Invitrogen. CA) directly from various 
cell lvsates as described in Example 3 Paragraph (4). Samples were dissolved in 10 \i\ 
DEPC-treated water, and cDNA was synthesized using the Superscript kit (Gibco/BRL Life 
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Technologies, MD). Resultant cDNA was extracted once with phenol, precipitated in 
ethanol. and d issoived in 10 |il of DEPC-treated water. 

One ill of solution containing cDNA was mixed with J \i\ each of 10 mM dATP, 
dCTP, dGTP and dTTP (Pharmacia LKB Biotechnology, NJ), 1 ^1 of 25 mM MgCU. 1 \il 
5 each of G protein-specific PCR primers. 0.5 (il of Taq polymarese (Promega. WI), and 5 \l\ 
of PCR buffer (Promega) to a final volume of 50 (il. Sense and antisense G protein-specific 
PCR primers are exibited below 

AG CA CCATTGTG AAG CA G A TG A (sense) (Seq ID No. 3) 

CTCTG G CCTCCCA CATC AAA CA (5' to 3' antisense) (Seq ID No. 4) 
10 These sequences are known to be highly conserved among various G proteins, with 

this pair of primers being known to amplify Gs. Gi-1. Gi-2, Gi-3 and Go protein specific 
cDNA from data gathered in preliminary experiments. Followed by pre-denaturation at 
95oC for 10 minutes, PCR was carried out in a DNA thermal cycler (Perkin-Elmer Cetus, 
CT) with 30 cycles of annealing at 55©C for 1.5 minutes, extension at 72©C for 4 minutes, 
15 and denaturation at 95oC for 1.5 minutes, respectively. 

(b) Aenrose gel electrophoresis of amplified cDNA 

A 1.2% agarose gel containing 0.00005% ethidium bromide, 0.001 M EDTA and 
Tris-acetate. pH 8.0 was prepared. Ten \i\ of each test solution was mixed with 2 \x\ of gel- 
loading buffer (0.25% bromophenol blue. 0.25% xylene cyanol and 15% Ficoll (Type 400, 
20 Pharmacia LKB Biotechnology, NJ; in water) and heated at 65©C for five minutes prior to 

loading. Electrophoresis was performed at 6 V/cm for one hour. After electrophoresis, the 
cDNA bands were examined under UV light illumination. 

The gel was then immersed in 0.25 M HC1 for 30 minutes (depurination), 0.5 M 
NaOH containing 1.5 M NaCl for one hour (denaturation). followed by 1 M Tris-HCl, pH 
25 7.5 containing 1.5 M NaCl for one hour (neutralization). 

(c) Transfer of Amplified cDNAs onto Memhranes (Southern Blon 
Amplified cDNAs were transferred onto nylon membranes and immobilized by UV 

light as described in Comparison 1(2). 

(d) Preparation of Biotin-Laheled cDNA Probes 

30 Ten ng each of Gs and Gi-2 cDNA clones were amplified by PCR as described in 

(2)(a) above, with the modification that 20% of dTTP was replaced with biotin-labeled 
dUTP (Clontech). In preliminary Southern blot analysis, these biotin-labeled PCR products 
have been shown to specifically hybridize to Gs and Gi-2. respectively. 

( e ) Hybridization of Labe led cDNA Probes to Southern Membranes 
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As shown in Comparison Example 1 Par;c aph 4. biotin-labeied PCR products were 
heat-denatured at 95°C for 10 minutes, then hyiv.dized to Southern membranes, on which 
PCR-ampIified cDNAs from cell lysates were immobilized. 

(f) • Chemical Activity Measurement oi Labeled cDNA Probes 
5 Chemical activity of labeled cDNA probes was measured as described in Comparison 

1 Paragraph (5). 
(3) Results 

The results of Example 4(1) are shown in Figure 11. Polaroid film was exposed to 
chemiluminescent light derived from an oligo-(dT) probe by the invented method, in which 
10 Gi-2 and Gs protein mRNA from various cells were trapped onto the first oligonucleotide- 

immobilized plastic plate. Phosphate-buffered saline (PBS) was used as a control. As 
shown in Figure 1 1. mRNA was detected from all cell lines. 

The results of Example 4(2) are shown in Figure L2, in which Gi-2 and Gs cDNA 
probes were used on a Southern Blot. PCR-ampiified G protein specific DNA from four 
15 different cell lines were immobilized onto a membrane. Figure 12 is a copy of the X-ray 

film detailing the chemiluminescent results of this probe. As shown in Figure 12, rnRNAs 
are detected in all cell lines. Therefore, the invented method provides a high enough 
reliability to be used for quantification of mRNA. 

EXAMPLE 5 

20 Methods for Amnlifvinc G Proteins With PCR Primers 

(1) Materials 

The cDNAs of rat G protein a subunits (G^. G } _ y G kV G s and G 0 ) were provided 
by Dr. R.R. Reed (Johns Hopkins Univ.. MD). XZAP libraries of rat pituitary and intestine 
were provided by Dr. D.G. Payan (Univ. Calif. San Francisco). Kirsten murine sarcoma 

25 virus transformed rat kidney cells (KNRK), human IM9 B-lymphocytes and human Jurkat 
T-lymphocytes were obtained from American Type Tissue Culture Collection. Rockville, 
MD). Cell culture media. Superscript (Gibco/BRL, Gaithersburg, MD), reagents for PCR 
(Prornega. Madison, WI), ECL (Amersham, Arlington Height. IL), Genius, Lumi-Phos 530 
(Boehringer-Mannheim, Indianapolis. IN), FastTrack (Invitrogen San Diego, CA), X gtlO 

30 library of human HL-60 cells. biotin-dUTP. alkaline phosphatase-conjugated streptavidine 

(Clontech, Palo Alto, CA). dNTP (Pharmacia. Piscataway. NJ) were obtained from the 
designated suppliers. Other chemicals were purchased from Sigma (St. Louis, MO). 

(2) C?ll cmUmh? 
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KRNK cells were grown in Dulbecco's modified Eagles medium containing 10% fetal 
calf serum, 100 U/ml penicillin and 100 pg/ml streptomycin at 37°C in 5% C0 2 /95% air. 
Cells were fed every other day and passaged at 70-90% confluency with 0.1% trypsin in 
CA 2 + -Mg 2 + -free saline containing 0.02% EDTA. IM9 and Jurkat cells was grown in RPMI 
5 1640 containing 10% fetal calf serum. 100 U/ml penicillin and 100 Kig/ml streptomycin at 

37°C in 5% CO^/95%. Cell viability was more than 90% as assessed by the exclusion of 
trypan blue. 

(3) Primer design 

Rat clones of Ga proteins (G i4 (RATBPGTPB), G^ (RATBPGTPA), 
10 (RATBPGTP), G s (RATBPGTPD), and G 0 (RATBPGTPC) were retrieved from GenBank 

release 65.0 (HIBIIO. Hitachi America. Brisbane, CA). The nucleotide sequence similarity 
among these clones were then analyzed by the multiple alignment program (DNASIS. 
Hitachi). We have initially identified 7 highly conserved areas among them. These 
conserved nucleotide sequences were then analyzed against all mammalian sequences in 
15 GenBank in order to identify other similar sequences. The designed oligonucleotides (Seq. 

ID No. *) were synthesized by Genosys Biotechnologies (Woodlands, TX), and suspended 
in water at 100 pg/ml. 

(4) PCR 

One ill of the template DNA was mixed with 1 mM each of dATP, dGTP. dCTP and 
20 dTTP. 1 til of each PCR primers, 1 ^1 of 25 mM MgCU. 5 |il PCR buffer, and 0.5 p\ of Taq 

polymerase (18). PCR was then carried out in a DNA thermal cycler (model 480, Perkin- 
Elmer Cetus, Norwalk, CT) with 30 cycles of annealing temperature at ranging from 37° C 
to 65° C for 1.5 min, 72° C extension for 4 min followed by 95° C denaturization for 1.5 nun. 
In separate experiments, 35% of dTTP was replaced with biotin-dUTP in order to prepare 
25 biotin-labeled probes. 

(5) Southern blot 

PCR products were separated by electrophoresis in 1.2% agarose, and stained with 
ethidium bromide (19). Gels were then depurinated in 0.25 N HC1 for 30 minutes and 
denatured in 0.5 N NaOH containing 1.5 M NaCl for 30 minutes. The gels were then 
30 neutralized with 1.0 iM Tris, pH 7.6 containing 1.5 M NaCl for 30 minutes. Gels were then 

placed onto nylon membranes (MagnaGraph. MSI, Westboro, MA) prewetted in 10XSSPE 
for 10 min, and DNA was transferred onto membranes by positive pressure at 75 mmHg for 
60 minutes (Posiblot. Strata gene, La Jolla. CA). The DNA from the gel was then cross- 
linked to the membranes with ultraviolet light at 120 mjoules (Stratalinker. Stratagene). and 
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the membranes were incubated with hybridization buffer (ECL) containing 5% blocking 
reagent (ECL) and 0.5 M NaCl at 40°C for more than 1 hour. Heat denatured biotin- 
labeled PCR probes were then added, and hybridization was continued overnight. The 
membranes were washed four times for 15 minutes each lime with primary wash buffer (OJx 

5 SSPE. 36 w/v% urea, 0.4w/v% SDS) at 45-65° C then washed twice for 5 minutes with 

secondary wash buffer (2x SSPE) at room temperature, and were incubated with the 
blocking buffer (Genius) for at least 3 hours at room temperature. Alkaline phosphatase- 
conjucated streptavidine (1:5,000 dilution) was then added, and incubation was continued 
for an additional 30-60 minutes at room temperature. The membranes were washed four 

10 times for 15 minutes with buffer A (100 mM Tris, pH 7.5, 150 mM NaCl) at room 

temperature, were washed for 2 minutes once with buffer C ( 100 mM Tris, pH 9.5, 100 mM 
NaCl. 50 mM MgCl 2 ), and soaked in Lumi-Phos 530 for approximately 1-2 minutes. The 
membranes were then wrapped with transparency films, and chemiluminescent signals were 
allowed to expose X-ray films (XAR-5, Kodak, Rochester. NY) for between 10 minutes and 

15 1 hour 

(6) mRNA preparation and cDNA synthesis 

The cells were washed with phosphate buffered saline three times, homogenized in 
lysis buffer (FastTrack), and then incubated at 45° C for 1 hour to eliminate any RNase 
activity. NaCl concentrations were adjusted at 0.5 M. and an oligo (dT) cellulose tablet was 

20 added to lysis buffer. Incubation was then continued at room temperature for an additional 

40 minutes. After oligo (dT) cellulose was washed with binding buffer (FastTrack) four 
times, bound mRNA was eluted with DEPC-treated water. Concentrations of mRNA were 
determined in a spectrophotometer (Hitachi. U-2000. Irvine, CA) at OD^q. The first strand 
cDNA was synthesized from a template mRNA in the presence of 50 mM Tris. pH 8.3, 75 

25 mM KC1, 3 mM MgCU, 10 mM DTT, 0.5 mM each of dATP, dCTP, dGTP, and dTTP, poly 

(dT) as a primer, and reverse transcriptase (Superscript) at 37° C for 1 hour. Second strand 
cDNA was then synthesized in the same tube, containing 25 mM Tris, pH 7.5, 100 mM KC1, 
5 mM MgCU. 10 mM (NH 4 ) 2 S0 4 , 0.15 mM P-NAD*. 250 \lM each of dATP, dGTP, dCTP 
and dTTP. 1.2 mM DTT, 65 U/ml DNA ligase. 250 U/ml DNA polymerase, and 13 U/ml 

30 RNase H (Superscript) for 2 hours at 16° C. Synthesized cDNAs were then extracted once 

with an equal volume of phenol:chloroform:isoamyl alcohol (25:24:1), precipitated with 
ethanol, and resuspended in FUO. 

(7) Graphic presentation 
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Data on Polaroid films and X-ray films was scanned by Stratascan (Stratacene) with 
optimization of signal-to-noise ratio, then edited with desk top publishing software 
(PageMaker. Aldus. Seattle, WA). 

EXAMPLE 6 

5 Characterization of Cellular Components of G Protein 

a Subunits Using PCR 

(l) introduction 

Using the methods of Example 5. we have now identified and characterized two 
highly conserved oligonucleotide sequences among five different a subunits of G proteins, 
10 SEQ ID NO:7 and SEO ID NO:8. which can be used as PCR primers. We have found that 

the sequences of all the subclasses of G proteins can be amplified under the same PCR 
conditions using these two oligonucleotides as PCR primers, including both G protein 
sequences obtained from a mixture of rat Get protein clones and cDNAs derived from 
various human tissues. 

15 Although G proteins have been analyzed extensively from a biochemical and 

immunological point of view using various antibodies, antibody production without any cross- 
reactivity among various subclasses or with high species specificity has been quite difficult 
to obtain. Therefore, recent experiments have focused on Northern blot analyses to identify 
G protein-specific mRNA from various tissues or cells in different species. However, 

20 Northern blots require experienced handling and protection from RNase contamination, in 

addition to a large amount of starting cellular materials. 

In contrast to these conventional methods, PCR technology is more convenient and 
practically useful, because it requires less material than a Northern blot analysis and has 
great sensitivity. However, it is difficult for PCR to quantify the amount of DNA or mRNA 

25 in starting materials. Interestingly, the final PCR products obtained using the novel G- 

protein PCR primers of the present invention reflect the relative composition of each of the 
subclasses of Get proteins present in the starting materials. This is probably because the five 
different Ga proteins cDNAs are amplified at a similar rate with a single set of PCR 
primers under the same PCR conditions. If known mixtures of each of the subclasses of Ga 

30 protein clones are assayed together with unknown test samples, as shown in Fieure 13, the 

relative composition of Ga proteins can be determined fairly precisely. Therefore, the 
present method is ideal tor the characterization of Ga proteins in various tissues and cells. 

Performing PCR with the primers of the present invention is also useful in clinical 
and diagnostic assays in the detection of disease. Since G protein abnormalities have been 
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associated with hereditary diseases, cancer, forms of diabetes, and other diseases, the 
present PCR primers for detecting and quantifying G proteins can be used to detect these 
diseases and assess their severity. 

In the present study, we have identified G^, Gj_ 2 , G ; _ 3 , G s and G 0 using the PCR 
5 primers of the present invention. Although recent cloning has identified more subclasses 

of G proteins, all of these newly identified G protein cDNAs showed high degree of 
homology to other known G proteins. Therefore, it is expected that the primers of the 
present invention will amplify these subclasses as well, and that the present PCR the present 
PCR technique can also be applied to these new G proteins. Moreover, this PCR method 
10 can be utilized to clone unique G protein genes as well. 

(2) Primers 

We have designed two 22-mer oligonucleotides (G 2 and G 4 -SEQ ID NOS:700 and 
701) as PCR primers. As shown in Table 4 below, these oligonucleotides contain sequences 
which are highly conserved among five different Gee protein cDNAs. having only 0-4 base 

15 mismatches per sequence. No mismatch was found in the 4 bases at the 3* end of the G 2 - 

sense and G 4 -antisense sequences. Furthermore, G>* and G 4 have no self-complementary 
sequences more than 3 base pairs in a row (data not shown). In order to analyze whether 
G 2 and G 4 are common to all the Get proteins, but not to other unrelated sequences, a 
homology search (DNASIS) of G 2 and G 4 sequences was carried<out against all mammalian 

20 sequences in GenBank. As a result, G-> and G 4 were found to be common to all the types 

of Ga proteins and rhodopsins of various species, but less homologous to other unrelated 
sequences (data not shown). 
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TABLE 4 

Two consensus oligonucleotides (G2 and G4) 
among five different cDNAs of G protein a subunits 



Consensus sequence (# of mismatch) [SEO ID NO:] 



C2 G4 

ACCACCATTGTCAAGCAGATGA [700] Length (bp) TGTTTGATGTGGGAGGCCAGAG (701) 

10 AGCACaATTGTGAAGCAGATGA (1) [702] 476 TGTTTGACGTGGGAGGCCAGAG (1) £8053 

Gi *2 AGCACCATcGTcAAGCAGATGA (2) [703] 479 TGTTTGATGTGCGtGGtCACcC (3) [806] 

Gi-3 AGtACtATTGTGAAaCAGATGA (3) [704] 476 TGTTTGATGTaGGtGGCCAaAG (3) [707] 

Cs AGCACCATTGTGAAGCAGATGA (0) C705] 524 TGTTcGATGTGCGcGGCCAGcG (3) [708] 

GO AGCACCATTGTGAAGCAGATGA <0> [7063 47V TGTTTGAcGT tGGtjGGCCAGcC (4) [709] 



TABLE 5 

Nucleotide sequence similarity of PCR products 
among 5 different Go protein clones. 



20 Gi-1 Gi-2 Gi-3 Gs GO 

Gi -1 * 73.4 76.8 50.2 64.0 

Gi * 2 73.4 - 67.8 52.3 65.9 

Gi-3 76.8 67.8 - 47.6 62.6 

25 Gs 50.2 52.3 47.6 . 55.6 

Go 64.0 65.9 62.6 55.6 

(3) PCR nf rat cloned Ga protein cDNAs 

The cloned rat G^, G L2 , G^, G s and G 0 cDNAs were diluted at 10 ng/jiL and 1 liI 
of DNA was mixed with 1 mM each of dATP, dGTP. dCTP and dTTP, 1 liI of each PCR 

30 primers, 1 nl of 25 mM MgCU 5 Li] PCR buffer, and 0.5 llI of Taq polymerase as described 

in the Materials and Methods. PCR was then carried out with 30 cycles of annealinc 
temperature at 45° C for 1.5 min, 72° extension for 4 min followed by 95° C denaturation 
for 1.5 min, 72° C extension for 4 min followed by ( J5°C denaturation for 1.5 min. PCR 
products were then separated in 1.2% agarose gel, and stained with ethidium bromide 

35 (upper panel). Marker indicates a 0.6 Kb fragment of Hindlll digested k DNA. DNAs 

were then transferred onto nylon membranes, hybridized with subclass specific biotin-labeled 
probes as described in the Materials and Methods. The membranes were washed with 
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mismatches between G 0 and G 4 is higher than others between G 4 and the sequences. 
If 1 or 2 of the 5 Get protein cDNAs were present in smaller quantities than the others, the 
amounts of amplified cDNA were relatively correlated with the starting concentrations of 
cDNAs (Figure 13, lane I, 2, 3, S, 9). Furthermore, if 1 of 5 of the Ga proteins' cDNA is 
5 more abundant than that of the others, this G protein gene was amplified more than others 

(Figure 13, lane 6, 7). 

(5) PCR of rat cDNAs from various tissues 

Rat XZAP cDNA libraries of pituitary (P) and intestine (1), and cDNA from the 
mRNA of rat kidney cells (KNRK) were used as PCR templates. Then. PCR and Southern 
10 blot were carried out as described in part (3) of this Example. 

Using this PCR method, Ga protein genes were amplified not only from cloned 
cDNAs, but also from various rat cDNAs (Figure 15). In XZAP cDNA libraries from rat 
pituitary glands and cDNA from rat kidney KNRK cells. G 0 was more abundant than G s , 
Gj_ 2 and G i0 , and G M was undetectable (Figure 15. lane 1, 2). XZAP cDNA library of rat 
15 intestine contained more G h2 . G,. 3 , and G s and G 0 (Figure 15. lane 3). 

(6) PCR of human cDNAs 

cDNAs from the mRNA of human IM9 B-Iymphocytes and Jurkat T-lymphocytes 
were used as PCR templates. Then, PCR was carried out with annealing temperature at 
45° C (1,3) or 55° C (2. 4). Southern blot were then performed <as described in part (3) of 
20 this Example. 

According to sequence analyses, the PCR products of rat G M , G^, G,_ 3 , G s and G 0 
sequence amplification exhibited a high degree of homology to human G protein cDNAs, 
see Table 6 below). Furthermore, as shown in Figure 16. PCR with a pair of G. and G 4 
primers could amplify 500 bp DNA from cDNAs of human IM9 and Jurkat cells. Unlike 
25 rat cDNAs (Figure 15), both 1M9 and Jurkat cells contained all the subclasses of Ga 

proteins (Figure 16). However, G i0 is relatively more abundant in IM9 cells, while G s and 
G 0 were more in Jurkat cells than IM9 cells (Figure 16). 
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Nucleoride sequence similarity of PCR products between 
rat and human G proteins 



5 Length (bfjyo. of 

Rat Human mismatch % Similarity 



Gi-l 476 476 62 87.0% 

Gi-2 479 479 44 90.8% 

10 Gi-3 476 476 40 91.6% 

Gs 524 482 52. 89.2% 

Go 479 479 39 91.9% 



It will be evident to one having ordinary skill in the art that a variety of sequences 
could serve as sense or antisense primers for PCR methods or as probes for the detection 

15 of mRNA as described herein. A method for identification of such sequences that are 

either common to a variety of G proteins or specific to a particular species is provided 
hereinbelow. In the preferred embodiment of this method of identification, a computer 
program is used to identify the sequences. Through use of such a program, we have 
identified a large number of both common and specific primers and probes. Provided as 

20 Tables 7 through 13 are various sense sequences identified through the use of such a 

program that are useful as G protein probes and primers. 

All of the sequences listed in these tables are useful within the context of the PCR 
methods of the present invention. The complementary antisense sequences are also useful 
in certain aspects of the invention. As will be known having ordinary skill in the art, for 

25 common probes that are similar, but not identical to target sequences, stringency conditions 

can be varied (e.g. by changes in temperature and salinity) so that such probes will hybridize 
or fail to hybridize with a particular target sequence. Thus, also included within the present 
invention are sequences thar are capable of hybridizing with the same sequences as either 
the sense sequences listed or their anti-sense counterparts. 

30 Additional probes for G protein also include the following: 

Common G protein prohes 

S'-CTCTGGCCTCCCACATCAAACAo (SEQ ID NO: 139) 
5 '-TCATCTG CTTC A CA ATGGTGCT- 3 * (SEQ ID NO: 140) 
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G:'-l 5 -GTTTTCACTCTAGTTCTGAGAACATCo" (SEO ID NO: 141] 

Gi-2 5'-CAAAGTCGATCTGCAGGTTGC-3' (SEQ ID NO: 142) 

5'-ATGGTCAGCCCAGAGCCTCCGG-.V (SEQ ID NO: 143) 

Gs 5-GCCTTGGCATGCTCATAGAATT-3" (SEQ ID NO: 145) 

S'-TTCATCCTCCCACAGAGCCTTG-.V (SEQ ID NO: 146) 
Go 5'-CGCATCATGGCAGAAAGCAG-3' (SEQ ID NO: 147) 
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Tafcle 7. G protein common primer (G2-S) 



G2-S 

AGCACCATTGTGAAGCAGATGA 

5 



Human SEQ ID NO: 

Gi-1 HUMGNBPAI 148 ::T 

Gi-2 HUMGIAA 149 ::: 

10 Gi-3 HUMGIAB 150 : : : 

Gs HUMGNPAS 151 : : : 

Go HUMGOAQ01 152 : : : 

Rat 

15 Gi-1 RATBPGTPB 153 : : : 

Gi-2 RATBPGTPA 154 : : : 

Gi-3 RATBPGTP 155 ::T 

Gs RATBPGTPD 156 : : : 

Go RATBPGTPC 157 : : : 

20 Gx RATGXA 158 : : : 



Highest matched sequences in GenBank 

HUMADECYC adenyl cyclase 159 :::::::::::::::::::::: 
RATACOAlacyl-coA oxidase 160 GC: : : : : :G: : : : : :A: : : : : : 
25 HUMTGASEtransglutaminase 161 CT: ::::::::::::: :CC:AC 
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Table 8. G protein common primer (G4-AS) 



G4-AS 

SEQUENCE ID NO: TGTTTGATGTGGGAGGCCAGAG 



Human 

Gi-1 HUMGNBPA 162 :::::::::::::::: T :::: : 

Gi-2 HUMGIAA 163 ::::::::::::: T :: T ::: C : 

Gi-3 HUMGIAB 164 :::::::::: A: :T: : : : : A: : 

Gs HUMGNPAS 165 ::::::: C : : : : : T: : : : : : C: 

Go HUMGOAQ01 166 ::::::: C: : C: :::::::: C: 

Rat 

Gi-1 RATBPGTPB 167 ::::::: C ::::::::::::: : 

Gi-2 RATBPGTPA 168 ::::::::::::: T :: T ::: C : 

Gi-3 RATBPGTP 169 :::::::::: A: :T: :: : :A: : 

Gs RATBPGTPD 170 :::: C :::::::: C: ::::: C : 

Go RATBPGTPC 171 :::::::C::T::G::::::C: 

Gx RATGXA 172 :: G : G ::::::::::: G :::: : 

Highest matched sequences in GenBank 

HUMLDLRRL LDL-receptor 17 3 CTGA: :::::::::::: : TACT 

MUSHEPGFA hepatocyte growth 174 GAG :G:::::::T::A:::::: 

HUMKEREP epidermal keratin 17 5 AG :::::: AT ::: GG :::::: : 
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Table 9. Gi-1 protein specific, human-rat common 
primer (Gi-1) 



10 



15 



20 



SEQUENCE ID NO: 



Human 








Gi- 


1 


HUMGNBPAI 


176 


Gi- 


2 


HUMGIAA 


177 


Gi- 


3 


HUMGIAB 


178 


Gs 




HUMGNPAS 


179 


Go 




HUMGOAQ01 


180 


Rat 








Gi- 


•1 


RATBFGTPB 


181 


Gi- 


■2 


RATBPGTPA 


182 


Gi- 


-3 


RATBPGTP 


183 


Gs 




RATBPGTPD 


184 


Go 




RATBPGTPC 


185 


Gx 




RATGXA 


186 



Gl-1 

GATGTTCTCAGAACTAGAGTGAAAAC 



TT: :::::::: CT: CCCCTGTCCCCT 
• ••••••• • f* • • rz • * • • • • f • • 

AC: : :G:C: : :T: : :TCCTG: :C: :G 



:G: :GC:G: :CC:T 

• • lt*b» • Vj • • • 

:GCACAATTA:TTA: : : 
:CA:C: : :C: : : :C: :G 
TCA: : : :GAG: :C: : :A:CC 



:G 
:G 



CG 
CA 
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Table 9A. Human- rodent conmon Gi-1 specific probes. 



Name Gi 1-735 (5'-3') Gi1-1l31 (5'-3') SEQ ID NOS: 

in GenBank GATGTTCTCAGAACTAGAGTGAAAAC TGTCAGTTTGAAGACCTCAATAAAAr, (respectively) 



10 Rax 

Gs RATBPGTPO A-GCACAATTA- TTA CG AACAGAAA - A AA - A - - - G - - - T 

Gi-1 RATBPGTPB 

Gi-2 RATBPGTPA G- -GC-G- -CC-T G-* A-CA G G C- 

Gi-3 RATBPGTP TC-G--G G-- --C T--G--CCG--- 

15 GO RATBPGTPC --CA-C- --C-- • -C-C--C ACA AG-AAA- -CCGCTC 

Human 

Gi-1 HUMGNBPA 

GenBank** 

HUMCPH192 C---C---G--C---G-—A G 

20 MUSGBPA A-C---CC--T--C 



22%, 


225 


226, 


227 


228, 


229 


230, 


231 


232, 


233 


234, 


235 


236 




237 





25 



Size of PGR products using Gil -735 and Gil-1131 as primers. 
The highest matched sequences in GenSank (release 68.0) 
next to the C protein DMAs. 
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Table 10. Gi-2 protein specific, human-rat common 
primer (Gi-2) 



Gi-2 

5 SEQUENCE ID NO: GCAACCTGCAGATCGACTTTG 



20 



Human 
















Gi-1 


HUMGNBPAI 


187 


: G : GGT: : A: 










Gi-2 
Gi-3 


HUMGIAA 
HUMGIAB 


188 
189 












: ACGG: :AA: 


: : :T: 








Gs 


HUMGNPAS 


190 


AGT: : :A: :T 










Go 


HUMGOAQ01 


191 


CTTCT: : : : : 


: : G : TG : : 


> : :C 


Rat 
















Gi-1 


RATBPGTPB 


192 


: G : GAT : : A: 


A: : : 


• • 






Gi-2 
Gi-3 


RATBPGTPA 
RATBPGTP 


193 
194 












AA •••G*«*« 


:T:T 


;TT 






Gs 


RATBPGTPD 


195 












Go 


RATBPGTPC 


196 


ATTCT : : : : : 


: : A : TG : 


: : :C 


Gx 


RATGXA 


197 


• T • • • A • C a a 


:T:T 


:T: 


: :C 
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Table 10A. Hunan- rooent common Gi-2 specific prooes. 



10 



15 



20 



Name 
in GenBank 



Gi2-742 C5'-3»> 
AAGATGTTTCATCTGGGTGGTC 



Gi2-1102 CS 1 ^') 
AAGGAGATCTACACCCACTTCA 



Rat 






Gs 


RATBPGTPO 


G-C-A-A ATCAAG- 


Gi-1 


RATBPGTPB 




Gi-2 


RATBPGTPA 




Ci-3 


RATBPGTP 


--A A C- 


Go 


RATBPGTPC 


•G-C C--T--G--C- 


Mouse 






Gi-2 


HUSGI 




H Linen 






Gi-2 


HUHGIR 




GenBank 


*• 






HUMCMPF 


C TGC 




MUSHSP84A 





■-C--G CG-GC-ACGC 

A- • T C 



CTA- 



-T 1- 

-C-CTCA-C 



StO ID NOS: 
(respectively) 



-A-G- 



238, 


239 


240, 


241 


242, 


243 


244, 


245 


246, 


247 


248, 


249 


250, 


251 


252 




253 





* : Size of PCR products using Gi2-742 and G i 2 - 1 1 02 as primers. 
25 **: The highest matched sequences in GenBank {release 68.0) 

next to the G protein DMAs. 
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Table 11. Gi-3 protein specific, human-rat common 
primer (Gi-3) 



SEQUENCE ID NO: 



Gi-3 

TCTTCGGACGAGAGTGAAGAC 



10 



15 



20 



Hutna n 

1 1 UiUU 1 1 








Gi- 


1 


HUMGNBPAI 


198 


Gi- 


2 


HUMGIAA 


199 


Gi- 


3 


HUMGIAB 


200 


GS 




HUMGNPAS 


201 


Go 




HUMGOAQ01 


202 


Rat 








Gi- 


1 


RATBPGTPB 


203 


Gi- 


2 


RATBPGTPA 


204 


Gi- 


3 


RATBPGTP 


205 


Gs 




RATBPGTPD 


206 


Go 




RATBPGTPC 


2 07 


Gx 




RATGXA 


208 



• • • f* 21 * 21 • • r P* • • « « 

• • • >■■» *»• • X* 

G ■ • f\ • • • • • CC • C • i 


: : :A: : 
A: : : : : 






: : :C: : : :GA: : : CGTG; A: : 


* • • O 21 • 21 • • T 1 • • • • 

• • • w A • i\ 


: : : :A: : 






CGCAT: : :G: :C: C 


: : :T:TG 
: : : : CCA 



AAC: : : 



G: A: 



: CACCAT : 
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Table 11A. Human-rodent common Gi-3 specific probes. 



10 



15 



20 



Name 
in GenBank 


Rat 




Gs 


RATBPGTPD 


Gi-1 


RATBPGTPB 


Gi-2 


RATBPGTPA 


Gi-3 


RATSPGTP 


Go 


RATBPGTPC 


Hunan 




Gi-3 


HUMGTPBP 


GenBank 


** 




HUMAPOB 




HUM 1 HSR 02 



Gi3-407 C5 ' -3 • > 
TTGTTTTAGCTCGCAGTGCTGA 



Gi3-730 f5'-3') SFO ID NOS: 

GAGGGAGTGACAGCAATTATCT ( respect i vely ) 



CC-CAAGT - GATC* 



-TC A-T-AT T--C--C- 



A GA- 


GGC-CTA-T 








CCTACAC- 


--A--AT C 




-T--C--G 




CG 


-T - - AGTC-TTAC- 


TC - - 


AT-* 


* TCAACGAC 


G C- 


- A T C 



















254, 


255 


256, 


257 


258, 


259 


260, 


261 


262, 


263 


264, 


265 


266 




267 





25 



Size of PCR products using Gi3-407 and Gi3-730 as primers. 
The highest matched sequences in GenBank (release 68.0) 
next to the G protein ONAs. 
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Table 12. Gs Protein specific, human-rat common 
primer (Gs) 



SEQUENCE ID Nu: 



Gs 

AAT"TCTATGAGCATGCCAAGGC 



10 



Human 
Gi-1 
Gi-2 
Gi-3 
Gs 



HUMGNBPAI 2 09 

HUMGIAA 210 

HUMGIAB 211 

HUMGNPAS 212 



GGGCGG: : : : T: : : C : : : : : CT 
: : : ATG : : : : : GCA: : ; : GCTA 
C : G : : : : ACTA: : : T : : : : CTC 



15 
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Table 12A. Human- rodent common Gs specific probes. 



in GenBank 



Gs-246 (S'^') 
CCCAACAAAAAGATCGAGAAGC 



Gs-a2t (5'-3') 
GGACAAAGTCAACTTCCACATG 



SEQ ID NOS: 
C respectively) 



10 



15 



20 



25 



Rat 

Gs 

Gi-1 
Gi-2 
Gi-3 
Go 

Mouse 
Gs-1 
Gs-2 

Human 
GS-1 
Gs-2 

GenBank 1 



RATBPGTPO 

RATBPGTPB 

RATBPGTPA 

RATBPGTP 

RATBPGTPC 

HUSGS 
MUSGTPAMU 

HUHGSA1R 
HUMGSA2R 

HUMMBP2 
RATAFPM 



CG-G---G-T CCGCA 

CAA-GG CA-A 

TGAGGA T-AA T 

TT7GG-G-G 7A T 



C-CGGC-G-GG-GCGC 

A C-C--C-G-T---T 

-TCAGC — A-TA CC 

CCGAGG* -CGGCGA-G 



•A--C- — 



-A-- 



268, 


269 


270, 


271 


272, 


273 


274. 


275 


276, 


277 


278, 


279 


260, 


281 


282, 


263 




285 


286 




287 





Size of PCR products using Gs-246 and Gs-824 as primers. 
The highest matched sequences in GenBank (release 66-0) 
next to the Gs protein DNAs. 
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Taole 13. Go Protein specific, human-rat common 
primer (Go) 



Go 

5 SEQUENCE ID NO: CTGCTTTCTGCCATGATGCG 



Human 



Gi- 


1 


HUMGNBPAI 


213 


GAC : : : : : 


: TAA: : : TGACA 


Gi- 


2 


HUMGIAA 


214 


: : :TGA: : 


» C • .A • G • C a • • • 


Gi- 


3 


HUMGIAB 


215 




: AGA : CTTC ACA 


Gs 




HUMGNPAS 


216 


TCAACGA : 


• • • • G • • • C^lTC 


Go 




HUMGOAQ01 


217 







Rat 



Gi-l 


RATBPGTPB 


218 


: : :A: 


: : :CCA:G:T: : :TT 


Gi-2 


RATBPGTPA 


219 


T: : : : 


GA:CC: : :GCG: : : : 


Gi-3 


RATBPGTP 


220 






Gs 
Go 
Gx 


RATBPGTPD 
RATBPGTPC 
RATGXA 


221 
222 
223 




. :G: : : AG : A : G : C : T 


:A: : 


: C : T : : : : C : : : CTG : 
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TAbie 13A Human-rodent common Gc specific proba*. 



in GenBani 



Go-1224 ( V-3' ) 
AGCACATCCTCCGAACCAG 



Go-1397 (5' -3' ) 
TATGACCAGGTGCTCCACG 



SEC ID NOS 

l respectively ) 



10 



15 



20 



25 



Gs 

Ci-1 

Ci-2 

Gi-3 

Gc 
Mouse 

Go -A 

Go-B 
Hamster 

Gc-1 

Gc-2 
Human 

Go 

GtnBanlc* 



RATBPGTPD AAAAAC-A-CT 

RATBPGTPB G--G-CG GG 

RATBPGTPA CCAT --TTCCT--A 

RATBPGTP G-CCA-AGA GGA 

RATBPGTPC - 

MUSGOASA 

MUSGOASB 

HAMHITAOl 

HAMHITA02 

HUMGOAH2) 



MUSPC3AA 
RATACACA 



G 



-T--TC--- 



AG AC--G-TTC 

7 T--TGCT- 

-TCTGTT GG-GA 

A TA-T-AATTT 



2B8. 


289 


290. 


291 


292, 


293 


294, 


295 


296, 


297 


296. 


299 


30C . 


301 


302. 


303 


30*. 


305 


306, 


307 


308 




309 





Size of PCR products usins Go- 122* and Go-1397 as primers. 
The highest matched sequences in GenBank (release 66.0) 
next to the G protein DMAs . 
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Substance P 

Another example of clinically useful measurement of mRNA would be for the 
mRNA of substance P. Substance P is a neurotransmitter. Because the nerves expressing 
substance P are involved in pain receptor pathways, the ability to measure the mRNA of 
5 substance P would be useful in the development of analgesics. 

EXAMPLE 6 

Detection of Human Substance P Receptor mRNA From 
Substance P Receptor cDNA-Transf erred Cells. Vector- 
Transfected Cells, and Evaluation of Their Reliability 

10 (1) Preparation of the First Nucleotide Probe 

The first nucleotide probe used for quantification of substance P receptor-specific 
mRNA was an oligodeoxyribonucleotides containing a 5' amino group. This probe was 
manufactured using a DNA synthesizer (model 380B, Applied Biosystems, CA) and is shown 
below: 

15 5'-NH 2 -GGA CTT ATG AGA AAG CGT ACC A-3' (Seq ID NO:5). 

(2) Immobilization of the Second Nucleotide Probe Onto Insoluble Supports 

The first nucleotide probe were dissolved in DEPC-treated water at a concentration 
of 1 ng/jil. then mixed with the EDC/Sulfo-NHS solution described in Example 1 Paragraph 

(2) in a ratio of 1:25 (Vol:Vol). 50 \i\ of the above mixed solution was added to each well 
20 of a Sumilon microliter plate (MS-3796F, Sumitomo Bakelite. JAPAN) and the plates were 

incubated overnight at room temperature. Following incubation, the reaction solution was 
removed with an aspirator. 

(3) Preparation of Standard mRNA 

A plasmid vector pKR2 containing human substance P receptor cDNA (provided by 
25 Dr. S. Nakanishi, Kyoto University, Japan) was linearized by Smal (Promega, WI) digestion. 

Substance P receptor-specific mRNA was prepared by the method described in Example 1 
Paragraph (3) using 77 RNA polymerase (Promega. WI). Concentrations of the resultant 
mRNA were determined by spectrophotometry and then diluted with DEPC-treated water 
to prepare standard mRNA solutions as described in Example 1 Paragraph (3). 
30 (4) Preparation of Cell Lvsates 

One line of human T lymphocyte-derived Jurkat cells which had been transfected 
with substance P receptor cDNA was examined for transcription of substance P receptor 
mRNA (hereafter called JSP. provided by Dr. D. G. Payan. University of California, San 
Francisco. CA). Another line of Jurkat cells which had been transfected with a vector that 
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didn't express substance P receptor (hereafter called JVec. a kind gift of Dr. D.G. Payan, 
ibid) was used as a negative control. Cells were collected by centrifugation, washed with 
PBS, and resuspended in lysis buffer at 1 x 10 7 cells/ml as described in Example 2. The 
ceiis were sheared by repeated passage through a 21 gauge needle, then incubated in a slow- 
5 shaking waterbath at 45° C for 1 hour. 

(5) Hybridization of the First Nucleotide Probe to mRNA 

In order to remove possible RNase contamination in the microtiter wells containing 
the bound first nucleotide probe, 250 \i\ of lysis buffer was added and incubated at 45° C for 
1 hour as described in Example 1. The buffer aspirated off, and 50 \lI of lysis buffer 
10 containing 0.5 M NaCl and various concentrations of standard human substance P receptor- 

specific mRNA was added into each well, and incubated at 50°C for 1 hour to allow 
hybridization as described in Example 1. 

(6) Hybridizati on of the Labeled Second Nucleotide Probe to a PolWAl Tail of mRNA. 
and Measuremen t of Chemical Activities of the Labeled Second Nucleotide Probe 

15 After hybridization of the first nucleotide probe to bound mRNA as described in 

Paragraph (5) above, the hybridization solution was removed by aspiration and individual 
wells were washed once with 250 ^1 of lysis buffer. Five \i\ of lysis buffer containing 0.5 M 
NaCl and 1 \i\ of the biotinylated second nucleotide probe solution (35 pmol/^1) was added 
into each well and incubated at room temperature for an additional 1 hour to allow the 

20 second hybridization as described in Example 1. 

Chemical activities of the labeled second nucleotide probe was also measured as 
described in Example 1. 

(7) Preparation of Calibration Line 

The calibration line of standard substance P receptor-specific mRNA is shown in Fig. 

25 9. 

(8) Evaluation of Reproducibility 

Concentrations of substance P receptor-specific mRNA in JSP and JVec cell lysates 
were measured twice. Table 7 shows concentrations of mRNA obtained from the calibration 
line as described in Paragraph (7). Cell numbers from each sample were counted by a 
30 phase hemacytometer (American Scientific Products. IL). As shown in Table 14, substance 

P receptor-specific mRNA concentrations were found in JSP cells that had been transfected 
with substance P receptor cDNA, but were not detected from JVec cells transfected with 
only a control vector. This shows that the invented method provides good quantification of 
a specific mRNA with high reliability 
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TABLE 14 

Concentrations of Substance P Receptor-Specific mRNA in Jurkat Cells Transfected with 
Substance P Receptor cDNA or Onlv a Control Vector. 

5 



No. of Assay 


mRNA Concentration 




JSP [pg/cell] 


JVec [fg/cell] 




(Transfected) 


(Control) 


1 


135 


<0.1 




62.4 


<0.I 



Oligonucleotide Primers and Probes for the Defection of Jun Oncogenes 
15 One aspect of the present invention relates to DNA primers and probes for the 

detection of jun oncogenes, such as by PCR. Expression of the jun gene can serve as a 

marker for cell growth. 

It is well known that certain oncogenes, such as jun or fos are most rapidly expressed 

when cells are stimulated with outside stimuli to proliferate. Therefore, detection or 
20 quantification of levels of expression of jun or fos oncogenes will be a good marker for 

cellular mitogenic activity. 

Jun oncogenes are one of oncogenes first reported by Maki, et al. (Proc. Natl. Acad. 

Sci. USA. 84:2848-2852, 1987), and current intensive molecular cloning identified at least 

three subtypes of jun oncogenes, such as jun-B, c-jun and jun-D. Based on nucleotide 
25 sequence comparison, three subtypes of jun oncogenes in both humans and mice are 

different (Ryder K. Nathans D, Proc. Natl. Acad. Sci. USA 85:8464-8467. 1988; Ryder K. 

et al., Proc. Natl. Acad. Sci.. USA 85:14*7-1491, 1988; Hattori K. et al, Proc. Natl. Acad. 

Sci. USA 85:9148-9152. 1988; Schueire J. et al.. Cell. 59:987-997, 1989). However, because 

it is still unclear how these three different jun oncogenes are involved in cell growth, one 
30 has to analyze all of the three genes during cell growth. 

Northern blot analysis is widely accepted in order to detect specific genes (Sambrook 

J. et al. Molecular cloning. A laboratory manual. 2nd eds. pp. 7.39-7.52), and some 

investigators have already adapted this method to detect jun oncogenes (Sherman, et al., 

Proc. Natl. Acad. Sci. USA. 87:5663-5666. 1990: Oursier MJ et aL Proc. Natl. Acad. Sci. 
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10 



USA, 88:6613-6617. 1991). In this method. mRNA is first purified from tissues/cell of 
researchers' interest, separated on agarose gel electrophoresis. After electrophoresis, 
mRNA is transferred onto membranes, and hybridized with radioactive probes to identify 
positive band(s) or. autoradiography. However, this meihuu is not so sensitive that it is 
difficult to identify specific signals if materials have only small quantity of eenes. 
Alternatively, because reverse PCR ("Molecular Cloning") is also used to detect a wide 
variety of genes from different tissues/cells, this technique is also applicable for the 
detection of jun oncogene. In this method, mRNA is first convert to cDNA by reverse 
transcriptase, then specific gene fragments are amplified by PCR by using one set of primers 
(sense and anti-sense primers). Amplified gene can be seen on agarose gel electrophoresis 
in the prepense of ethidium bromide. This method is much more sensitive than the 
conventional Northern blot analysis, but in order to detect three different subtypes of jun 
oncogenes, one has to prepare six primers. Also, because of different nucleotide sequences 
between humans and mice, one has to prepare human and mice primers separately. 
15 ° ne S oaI we accomplished with the present invention was the design of 

oligonucleotide sequences for both sense and anti-sense primers that could be readily used 
for detecting the minuscule amounts of the jun gene found in vivo. We also wanted these 
primers to be useful for detecting both human and murine jun sequences. Thus, we selected 
primers which fulfilled the below conditions: 
20 ( a ) A common nucleotide sequence among three subtypes of jun oncogenes (jun- 

B, c-jun and jun-D), with maximum 4 base mismatches among them; 

(b) A common nucleotide sequence among humans and mice, with maximum 4 
base mismatches between them; 

(c) At least 5 bases from 3' end are 100% identical without any mismatches 
25 among three subtypes ajid between humans and mice; 

(d) A nucleotide sequence length in both sense primer and anti-sense primer 
from 17 to 50; 

(e) A difference of T m between sense primer and anti-sense primer is within 2 
degree C, 

30 (0 L* ck of complementary structure with more than 4 bases long in either the 

sense or the anti-sense primer; 

(g) Lack of complementary structure with more than 4 bases long between the 
sense and the anti-sense primer; 
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(h) Both the sense and the anti-sense primer exist in the area of coding 
sequence; 

(i) The length of DNA to be amplified is greater than 200 bases: and 

(j) The nucleotide sequence homology of amplified genes among three subtypes 
5 of jun oncogenes is not greater than 809c. and after gene amplification by 

the PCR. one can characterize each subtype of jun ongenes by Southern 
Blot. 

DNA fracments (both sense and anti-sense primers) which satisfy the conditions as 
shown above were investigated and the several candidate oligonucleotides were synthesized 
10 The test results showed the sense primer S'-CCCTGAAGGAGGAGCCGCAGAC-S* and 
anti-sense primer 5-CGTGGGTCAAGACTCTGCTTGAGCTG-3* were desirable 
candidates. 

Table 15 shows the nucleotide sequence homology between a preferred sense primer 
and three subtypes of jun oncogenes in both humans and mice. 
15 TABLE 15 



20 



Sense primer 
5 1 -CCCTGAAGGAGGAGCCGCAGAC-3 1 (SEQ ID NO: 310) 



Mouse jun-B — T-T--A (SEQ ID NO: 311) 

c-jun A (SEQ ID NO: 3 12) 

Human jun-B — T-C A (SEQ ID NO: 313) 

25 c-jun T (SEQ ID NO: 314) 

indicates identical base to the sense priiner. 

In addition. Table 16 shows the nucleotide sequence homology between the sense 
primers and three subtypes of jun oncogenes in both humans and mice. 



30 
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TABLE 16 



Antisense primer 
3 1 -GTCGAGTTCGTCTTTCAGAACTGGGTGC-5 1 (SEQ ID NO: 315] 



Mouse jun-B C A- (SEQ ID NO: 316) 

c-jun T (SEQ ID NO: 317) 

jun-D T G-G-C (SEQ ID NO: 318) 

10 Human jun-B T C (SEQ ID NO: 319) 

c-jun T (SEQ ID NO: 320) 



-: indicates identical base to the sense primer. 

Also, Table 1 7 shows nucleotide sequence homology of amplified genes between each 
15 of subtypes of jun oncogenes. 

TABLE 17 

Subtypes Homology (%) 

20 jun-B & c-jun 72.0 

jun-B & jun-D 73.6 
c-jun & jun-D 75,3 



25 


Table 18 shows the 


homology between human and mouse subtypes. 






TABLE 18 




Subtypes 


Homology between human & mouse (%) 


30 


jun-B 


93.0 




c-jun 


91.9 



40 



The primer sequences can be modified in a variety of ways. For example, a primary 
amine residue, nucleotide sequence recognized by a restriction enzyme, or RNA promoter 
35 sequences for further modification after' PCR amplification, can be added to the 5' end of 

the primers. These modifications are represented herein by the designation "X". Thus, for 
example, a modified primer as shown in Table 16. can be represented as 
5 '-XCGTG G G TCA AG A CTTCTG CTTG A G CT G -3 ' ca n also be used, in which X represents 
the primary amine residue, nucleotide sequence recognized by the restriction enzvme. or 
RNA promoter sequences for further modification after PCR amplification. 

The attachment of a restriction site at the 5* end (X) is useful for cloninc of 
amplified genes. The attachment of RNA promoter sequences at 5* end (X) is useful in 
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RNA transcription and RNA transcription-based amplification. The attachment of primary 
amine at the 5* end is useful for coupling reaction to labeling compounds or solid supports, 
so that the amplified gene can be easily quantified. 

For the above-fiieriiioneu nucleotide sequence (X) which undergoes cutting by a 
5 restriction enzyme, instead of those nucleotide sequences with higher frequency such as Eco 

Rl or Bam HI. it is often desirable to use longer recognition sequences (e.g. 8 bases) which 
appear in lower frequency such as the nucleotide sequence, 5'-GCGGCCGC-3\ which is 
recognized by NotI or the nucleotide sequence, 5'-TTAATTAA-3\ which is recognized by 
Pad. The reason for this is the relatively common nucleotide sequences recognized by 
10 EcoRI or BamHI can result in false positive reaction during PCR amplification because of 

priming of primers to undesirable DNA sequences. 

The above-mentioned nucleotide sequence (X) containing an RNA promoter can be 
any of a variety of such sequences, such as the T7, SPfi or T3 RNA promoter sequences. 

The above-mentioned primary amine residues (X) can be added at the 5 1 end during 
15 oligonucleotide synthesis. 

The DNA fragments for the sense primer and anti-sense primer in this invention can 
be easily synthesized in DNA synthesizer which are available in the market. These 
synthesized oligonucleotides can be purified by high pressure liquid chromatography or gel 
electrophoresis. 

20 The test material to he analyzed is usually a total RNA or purified mRNA from cells 

or tissues. If desired, cells or tissues can be tested in their natural state without any 
pretreatment; however, in the case of drug testing, this can he accomplished by reacting the 
drug with cells or tissues in the test tube, or the drug can be administered (intravenous 
injection, subcutaneous injection, intramuscular injection, oral administration, intra- 

25 abdominal injection) to a laboratory animal, and after a set time period, cells or tissues can 

be removed for RNA/mRNA purificatipn. 

Methods of purification of total RNA or mRNA are well known in the art. 
Examples of such methods are the standard protocol described in Sambrook, et al., 
"Molecular Cloning, A Laboratory Manual, 2nd Ed., K Cold Spring Harbor Press, Cold Spring 

30 Harbor, NY ( 1989), and a commercially available kit such as FastTrack from Invitrogen (San 

Diego, CA). In either case, in order to avoid introducing RNase. hands should preferably 
be protected with vinyl gloves, and the instruments used for experiments not be touched 
with bare hands. Also, any glass containers to be used in the experiments should be heated 
prior to use at approximately 250 degrees C for at least 4 hours. Furthermore, 0.1% diethyl 
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pyrocarbonate (hereafter abbreviated DEPC) is added to the water which is to be used, and 
after incubation at 37 degrees C overnight, treated water is autoclaved. Other basic 
procedures are undertaken as outlined in the "Molecular Cloning," pp. 7.3-7.5, Vanadyl 
nbcnucleoSide complex can also be added to samples as described above. 
5 The method of synthesis of cDNA from mRNA template using reverse transcriptase 

is described in the "Molecular Cloning;* pp. 8.11-8.13. 

For PCR amplification, synthesized cDNA is mixed with sense primer, antisense 
primer, 4 types of deoxynucleotides (dATP, dCTP, dGTP and dTTP), Taq polymerase, 
inorganic salts, and other necessary materials, and PCR reaction is undertaken in the 

10 Thermal cycler (Perkin-Elmar Cetus). A typical PCR method is described in "Molecular 
Cloning/ pp. 14.2-14.33). 

In order to analyze the amplified gene, it is appropriate to use electrophoresis. 
After the amplified gene undergoes electrophoresis in agarose gel, the DNA is stained with 
ethidium bromide. The amplified DNA band will be then visible under fluorescent light as 

15 described in "Molecular Cloning." pp. E.3-E.4. After taking photographs of the DNA band, 

it is also possible to quantify the intensity of each bands by scanning photographs and 
analysis of scanned picture by commercially available system, such as Stratascan (Stratagene, 
La Jolla). 

Furthermore, after agarose gel electrophoresis, amplified genes can be transblotted 

20 onto membranes, and subtypes specific genes can be detected by hybridizing with labeled 

probes followed by exposing labeled signals, such as 32P or chemiluminescence to either 

Polaroid films or X-ray films (Southern blot). 

EXAMPLE 7 

Synthesis of Sense and Anti-Sense Primers and 
25 Amplification of Mouse Clones of Jun Oncogenes 

The sense primer SEQ ID NO;9 (s943-2) and the anti-sense oligonucleotide SEQ 

ID NO:10 (AS1 132-2) were synthesized by a synthesizer 380 B type (Applied Biosystems 

Co.). After treatment with ammonium hydroxide at 55°C overnight, synthesized 

oligonucleotides were dryed in a Speed-Vac (Savant Co.), and the concentration was 

30 adjusted at 1 microgram/ml with water, and was stored at -20 degrees C until use. 

1 microliter of one of the three types of mouse jun clones (jun-B. c-jun, or jun-D, 

obtained from ATCC), 1 microliter of sense primer. 1 microliter of anti-sense primer, 5 

microliters of 10X buffer for PCR (Promega), 1 microliter of 25 mM magnesium chloride, 

4 microliters of 10 mM dNTP mix. and 0.5 microliters of Taq polymerase (Promeea) were 
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missed, and watei was added up to a total volume of 50 microliters. After adding of two 
drops of mineral oLL PCR was undertaken using the thermal cycler, model 480. After 
reaction mixture was heated at 95 degree C for 10 min. PCR was carried out with the 
following cycles 30 times: annealing at 55 degrees C for 1.5 minutes, extension at 72 degrees 
5 C for 4 minutes, and denaturing at 95 degrees C for 1.5 minutes. 

After PCR, 10 microliters of the sample were mixed with 1 microliter of lOx loading 
buffer (025% bromophenol blue, 0.25% xylenecyanol FF, and 15% FicolL Type 400), and 
electrophoresis was carried out in 1.5% agarose gel containing 5 microgram/ml ethidium 
bromide. After electrophoresis, the amplified DNA bands were visualized by an ultraviolet 
10 light. 

As a result, it was found that in each of the cases where the clones jun-B, c-jun, and 
jun-D were used, a single band for the amplified DNa was observed at the position of 
approximately 270 bp. This indicates that the above-mentioned set of primers can recognize 
and amplify all three types of mouse jun oncogenes, such as jun-B. c-jun, jun-D. 

15 (1) Alternative Synthesis of PCR primers for amplification oFiun pene specific mRNA 

Sense (jun-s) and antisense Qun-as) oligodeoxyribonucleotides specific for jun 
oncogenes were synthesized in a DNA synthesizer as described above. The nucleotide 
sequences were as follows: 
jun-s 5 T -CCCTGAAGGAAGAGCCGCAGAC-.V (Seq rD NO: 11) 

20 jun-as S'-CGTGGGTCATGACTTTCTGCTTGAGCTGO* (Seq ID No: 12) 

(2) PCR amplification of cDNA 

To 1 ul of ss-cDNA as described in paragraph (IV), 0.1 ug (1 ul) each of jun-s and 
jun-as oligonucleotides, 5 ul of 10X PCR buffer (Promega). 1 ul of 25 mM MgCU, 4 ul of 
10 mM dNTP mixture (Promega), 0.5 ul Taq polymerase (Promega). and 36.5 \x\ of DEPC- 

25 water were mixed. Two drops of mineral oil were overlayered onto the reaction mixture to 

prevent evaporation. The reaction mixture was first heated at 95° C for 10 minutes, then 
PCR was carried out in a Thermal cycler (Model 480, Perkin-Elmer Cetus) with 30 cycles 
of annealing at 55° C for 1.5 minutes, extension at 72° C for 4 minutes and denaturing at 
95° C for 1.5 minutes as described above. 

30 Ten ul of the resultant PCR products were then mixed with 1 ul of 10X loading 

buffer, and applied to a 1.5% agarose gel containing 5 ug/ml ethidium bromide. 
Electrophoresis was carried out in a IX TBE buffer at 100 V/6.5 cm width for 
approximately 1 hour. After electrophoresis, fluorescent DNA bands were recorded onto 
Polaroid film. Lanes 1-8 of Figure 17 indicate duplicate experiments of a 1:1000 dilution, 
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full strength (equivalent to 1 ml of blood). 1:10. and 1:100 dilution respectively of a 

leukocyte lysate. As shown in Figure 17. the jun oncogene was amplified from the 1:10 and 

1:100 dilution displaying an estimated size of 370 base pairs. This size is similar to the 

372bp theoretical value. 

5 EXAMPLE S 

The effect of preireatment with EGF on expression of 
iun oncogenes in human mononuclear leukocytes 

(1) Pretreatment of human leukocytes with EGF 

40 ml of phosphate buffered saline (PBS) is added to 20 ml of heparinized human 
10 blood and mixed. 10 ml each of this sample is overlayered onto 3 ml of IsoLymph. then 

centrifuged for 30 minutes at 400 x g. After washing the pellet three times with PBS, the 
pellet is resuspended in 3 ml of PBS. 1 ml each of this sample is then placed in three tubes 
(No. 1 to No. 3). 1 ml of PBS Is placed in the fourth tube as control. The total of four 
tubes are incubated for 10 minutes at 37 degrees C then EGF (Epidermal Growth Factor) 
15 is placed in tube No. 1 at a final concentration of 30 ng/ml. After 15 minutes, EGF is 

placed in the same manner in tube No. 2, and incubated for another 5 minutes. After 5 
min, RNA is extracted from all 4 tubes simultaneously. Thus, the time period of the 
pretreatment of human leukocytes by EGF is 20 minutes for tube No. 1. 5 minutes for tube 
No. 2, 0 minutes for tube No. 3. 
20 (2) Extraction of RNA from cells 

The above four tubes are taken out and undergo centrifugation by microfuue for 10 
seconds. Then, the supernatant is discarded and the below described lysis buffer is added 
to the pellet, and after thorough mixing, this is incubated for 30 minutes at 45 degrees C. 

Contents of lysis buffer: 
25 10 mM EDTA pH 8.0 

0.59c SDS (bacteria removed by non-bacterial filter) 

0.2 M NaCl 

DEPC treated water 

RNA inhibitor 500 unites/ml 
30 Vanadyl Complex 10 mM 

Proteinase K 200 microgram/ml 
(3) Purification of mRNA 

After 5 M NaCl is added to the above cell lysates to obtain a final concentration of 
0.5 M, oligo (dT) cellulose (Stratagene Co.) is added, and undergoes a reaction for 30 
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minutes at room temperature. Then, after washing the cellulose in 10 ml binding buffer (20 
mM Tris-HCl, pH 7.6, 1 mM EDTA. 0.5 M NaCl) 5 times. 035 ml of DEPC treated water 
is added, and mRNA is eluted from the solid phase. Then. 53 microliters of 2 M sodium 
acetate and 2.5 times the volume of ethanol are added, and after cooling in a dry ice for 20 
5 minutes, this undergoes centrifugation at 15,000 rpm for 20 minutes. After washing the 

pellet in 75% ethanol one time, this is dried, and then dissolved in 10 microliters of DEPC 
treated water. 

(4) Synthesis of cDNA 

50 mM Tris-HCl (pH 83), 75 mM KCL 3 mM magnesium chloride, 10 mM DTT, 
10 0.5 mM dNTP (dATP, dCTP, dGTP. dTTP), 50 micrograms/ml oligo (dT) primer, and 

10.000 units/ml reverse transcriptase are added to 10 microliters of mRNA obtained from 
above to a total volume of 20 microliters, and this undergoes a reaction for one hour at 37 
degrees C. After the reaction, 20 microliters of a phenol:chloroform:isoamyl alcohol mixture 
is added, and cooled for 20 minutes in a dry ice to precipitate the cDNA. After 
15 centrifugation for 20 minutes at 10,000 rpm, the pellet is washed one time in 75% ethanol. 

Then, after drying, this is dissolved in 20 microliters of autoclaved water, and stored at -20 
degrees C. 

(5) PCR 

1 microliter each of sense primer and anti-sense primer (1 mg/ml) for jun gene 
20 amplification are added with 2 microliters of cDNA. After mixing this with 50 mM KC1. 10 

mM Tris-HCl (pH 8.4). 2.0 mM magnesium chloride, 100 micrograms/ml gelatin, and 0.2 
mM dNTP. 2.5 units of Taq polymerase are added (final volume is 50 ml). After reaction 
mixture was heated at 95 degree C for 10 minutes. PCR was carried out with the following 
cycles 30 times: annealing at 55 degree* C for 1.5 minutes, extension at 72 degrees C for 4 
25 minutes, and denaturing at 95 degrees-C for 1.5 minutes. 

(6) Agarose gel electrophoresis 

After completing PCR. 10 microliters of the reacted solution is taken and undergoes 
electrophoresis in the same method as described in Example 7. The results showed that in 
the leukocytes which did not undergo treatment with EGF (i.e., sample which underwent 0 
30 treatment time in tube No. 3) was found to have a minimal band for the amplified DNA at 

the position of about 270 bp size. However, it was found that the band of amplified DNA 
was increased after 5 min, then returned to the basal levels within 20 min. 

EXAMPLE 9 

Effect of pretreatment with PHA on jun gene expression in human leukocytes 
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ln the place of EGF. PHA (at a final concentration ot 10 micrograms/ml) was 
utilized and the pretreatment times were set at 0 minutes. 5 minutes. 15 minutes, and 30 
minutes. Other procedures were identical to those of Example S. As a result, with the 
pretreatment with PHA done at 15 minutes, it was found that the band for the amplified 
5 DNA at the position of approximately 270 bp was maximized, but (the time periods) 

afterwards, this band decreased. 

As with the other mRNAs that can be analyzed in accordance with the methods of 
the present invention, a variety of jun sequences could serve as sense or antisense primers 
for PCR methods or as probes for the detection of mRNA as described herein. A method 

10 for identification of such sequences that are either common to a variety of jun mRNAs or 

specific to a particular species is provided hereinbelow.' In the preferred embodiment of this 
method of identification, a computer program is used to identify the sequences. Through 
use of such a program, we have identified a large number of both common and specific 
primers and probes. Provided as Tables 19 through 24 are various sense sequences 

15 identified through the use of such a program that are useful as jun gene probes and primers. 

All of the sequences listed in these tables are useful within the context of the PCR 
methods of the present invention. The complementary antisense sequences are also useful 
in certain aspects of the invention. As will be known by those having ordinary skill in the 

20 art, for common probes that are similar, but not identical to target sequences, stringency 

conditions can be varied (e.g. by changes in temperature and salinity) so that such probes 
will hybridize or fail to hybridize with a particular target sequence. Thus, also included 
within the present iovention are sequences thar are capable of hybridizing with the same 
sequences as either the sense sequences listed or their anti-sense counterparts. Additional 

25 probes for jun genes include the following; 

Common jun gene probes 

5'-CCATGTCGATGGGGGACAGCGG-3 (SEQ ID NO: 13) 
5 , .CTGTTTAAGCTGCGCCACCTG-3 , (SEO ID NO: 14) 
^-GTCTGCGGCTCCTCCTTCAGGG^XSEQ ID NO: 15) 
30 S'-CGTGGGTCAAGACTTTCTGCTTGAGCTGO' (SEQ ID NO: 16) 

Specific probes 

B type: f'-CACTTGGTGGCCGCCAG-.V (SEQ ID NO: 17) 
C type: 5'-GAGCATGTTGGCCGTGG-.V (SEQ ID NO: 18) 
Human D type: S'-GATGCGCTCCTGCGTGT-.V (SEQ ID NO: 19) 
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Mouse D type: 5'-GCCTGTTCTGGCTTTTGAGGG-.V (SEQ ID NO: 20) 
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Locus 
GenSank 



Pos 



Sequences (5' -3 ' ) 
CCGCTGTCCCCCATCGACAT6G 



10 



15 



20 



25 



30 



35 



B: humjunca 


1189 


SEQ 


ID 


NO: 


21 


C: humjuna 


198: 


SEQ 


to 


NO: 


22 


D: hunjundr 




SEQ 


ID 


NO: 


23 


House 












B : Dius j unba 


1 079 


SEQ 


ID 


NO: 


24 


C: ffusjunc 




SEQ 


ID 


NO: 


25 


C: muscjun 


1646 


SEQ 


ID 


NO: 


26 


C: musjun 


1084 


SEO 


ID 


NO: 


27 


D: musjund 


927 


SEQ 


ID 


NO: 


28 


0: musjunda 


782 


SEQ 


ID 


NO: 


29 


0: musjundr 


793 


SEO 


ID 


NO: 


30 


Rat 












C: at j una pi 


1082 


SEQ 


ID 


NO: 


31 


C: ratrjg9 


2984 


SEO 


ID 


NO: 


32 


Chicken 












C: chkjun 


1470 


SEQ 


ID 


NO: 


33 


Quai I 












C: quljun 


1186 


SEQ 


ID 


NO: 


34 


Drosophila 












C: drojun 


1038 


SEQ 


ID 


NO: 


35 



Highest matched sequences in EMBL 
SDNAM2G Yeast NAM2 gene (SEQ 10 N0:36) 
PRK2TRFB Plasmid PK2 trfB ope (SEO ID NO:37) 
DMSYT D.meLanogaster synap (SEO ID NO: 38) 



---G A 

--C 

•--T G 

--TG A 

••C T 

-C T 

-C T 

G 

G 

--CT 

-CT 

*-C T-*T 

*-C T--T 

A-CG-TAAT T 

A*A GAAT 

GT GC-T- 

G-A UC 
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Table 20. Jjn-carmon antisense primer (AS1132-2). 



Lceus 



Sequences C 5 ' -3 ■ ) 
Pos . 



GTCCACCGCGTCGAOTTTGTC 



10 



15 



20 



25 



30 



Hunan 

B: humjunca 
C: humjuna 
D: humjundr 



House 

8: musjunba 
musjunc 
nuscjun 
musjun 
musjund 
mtis junda 
nusjundr 



C: 
C: 
C: 
D: 
D: 
0: 



Rat 

C: atjunap! 
C: ratrjg9 

Chicken 
C: chkjun 

Qua i I 

C: quljun 



1378 (SEQ ID MO: 39) 
2170 (SEO ID HO: 40) 
1132 (SEQ ID NO: 41 ) 



1268 

1535 

1273 

1533 

982 

1116 

971 



(SEQ ID NO: 
(SEQ ID HO: 
(SEQ ID NO: 
(SEQ ID NO: 
(SEQ ID NO: 
(SEQ ID NO: 
(SEQ ID NO: 



42) 
43) 
44) 
45) 
46) 
47) 
48) 



1271 (SEQ ID NO: 49) 
3173 (SEQ ID NO: 50) 



1659 (SEQ ID NO: 51) 



1375 (SEQ ID NO: 52) 



G G 

G--C-- 

-T G — C-- 

T C-- 

— T C-- 

T C-- 

GC-G 

A--T C- 

A — T C- 



35 



40 



Drosoohi la 
C: drojun 



1227 (SEQ ID NO: 53) 



Highest matched sequences in EMBL 
HSATFA Hunan raRNA. for ATF-a (SEQ ID NO: 54) 
ECDCM E. co Li dom gene CSEQ ID NO: 55) 
E CD CM A E . coli dom (SEQ ID NO: 56) 
OCIGKU Rabbit Ig germkine (SEQ ID NO: 57) 
OCIG05 Rabbit Ig k2 L chain (SEQ ID NO: 58) 
0CK1B4 Rabbit Ig kappa L (SEQ ID NO: 59) 
OCIGKCG Rabbit Ig kappa2 J-C (SE0 ID NO: 60) 
OC1GKC02 Rabbit Ig kappal J-C (SEO ID NO: 61) 



—A- 



GC-C--C — 



C- G---A-C 

AA ACCA 

AA ACCA 

G G-CGA 

G-CGA 

--G G-CGA 

G G-CGA 

G G-CGA 



45 
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Table 21. Jun-B specific probe (B-125B). 



CenBank 
B1258(5'-3') 



CTGGCGGCCACCAAGTG 



10 



15 



20 



hunan 

B: HUNJUNCA 
C: HUH J UN A 
0: HUHJUNOR 
Mouse 

8: MUSJUNBA 
C: KUSJUNC 
WJSCJUN 
HUSJUN 
D: MUSJUND 
MUSJUNOA 
MUSJUNDR 

Rat 

C: RATRJG9 
RATJUNAP 
Chicken 

C: CHKJUN 
Cuai L 

C: QULJUN 
Drosophi la 

C: DROJUN 



{ SEQ ID NO: 62) 

A-C--T---T (SEO ID NO: 63) 

A TCCAA {SEO ID NO: 64) 



A-T--C---T 

A-T--C---T 

A-T--C---T 

C CCCG 

T CCCA 

T GCCA 



(SEQ ID NO: 65) 

(SEQ ID NO: 66) 

(SEQ ID NO; 67) 

(SEO ID NO: 68) 

(SEQ ID NO: 69) 

(SEQ ID NO: 70) 

(SEQ ID NO: 71) 



A-C--7---T 

-C G- 



-T GCCAA (SEO ID NO: 72) 

(SEQ ID NO: 73) 

(SEQ ID NO: 7£») 



---ccc 



GGC AG TG- (SEO ID NO: 75) 



G T - -AT 



(SEO ID NO: 76) 



30 



35 



Homologous sequences in GenBank 

J04695 figure 2. Nueleo. G A G (SEO ID NO: 77) 

M2788^ figure 2. NucLeo. A GC (SEQ ID NO: 78) 

KUHCNPG2 green cone photo. T-AA ( SEO ID NO: 79) 

HUMCNPR2 green cone photo. T-AA (SEQ ID NO: 80) 

HUMPIGMF2 colour-blind pho. T-AA (SEO ID NO: 81) 
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Tabte 22. c-jun soecific probe (C2147). 



10 



15 



25 



30 



Subtype GenBank name 



C2147.DNA 
CCACGGCCAACATGCTC 



Human 

B- HUMJUNCA -G---A---C--CAGC (SEO 10 NO: 82) 

C- HUMJUNA (SEO ID HQ: 83) 

D; HUMJUNDR G-G-C----G (SEO ID NO: 84) 
House 

B- MUSJUNBA -G C CAGC (SEO ID NO: 85) 

C - MUSJUN (SEO ID NO: 87) 

MUSCJUN (SEO »> H0: Bfl > 

MUSJUNC (SEO ID NO: 89) 

D- HUSJUNO . — C----G-C-— G (SEO ID NO: 90) 

MUSJUNDR — . C — -G-C-—G (SEO ID NO: 91) 

MUSJUNDA -..-C-.--G-C — C (SEO ID NO: 92) 

^C: RAT JUNAP CSEQ ID NO: 93) 

RATRJC9 — -C (SEO ID NO: 94) 

Chicken 

C: CHUUN (SEO ID NO: 95) 
Quail 

C: QULJUK T (SEO ID NO: 96) 

D rosoph ila 

C; DROJUN CCTACA ACC (SEQ 10 NO: 97) 

Homologous sequences in GenBank 

MXBPALPA L . enzymogenes CG-- (SEQ ID NO: 98) 

MXBPALP L .en2ymogenes CG-- (SEQ ID NO: 99) 



35 
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Table 23. Hunan jun-D specific prooe (HUM0965) 



10 



15 



20 



25 



30 



Subtype 


GenSank name 


HUMD965 
ACACGCAGGAGCGCATC 








Hufndn 












8: 


HUM J UK CA 




( CP D 


1 V 


NU . 


C: 


HUMJUNA 


•GT-C G — 


( SEQ 


ID 


NO: 


0: 


HUHJUNOR 




(SEO 


ID 


NO: 


Mouse 












B: 


MUSJUNBA 




(SEQ 


ID 


NO: 


C; 


MUSJUNC 


-GT-T G-- 


(SEO 


ID 


NO: 




MUSCJUN 


-GT-T G--- 


(SEO 


ID 


NO: 




MUSJUN 


-GT-T G--- 


(SEO 


ID 


NO: 


D: 


MUSJUND 


A- -A 


(SEO 


ID 


NO: 




MUSJUNOA 


A--A 


(SEQ 


ID 


NO: 




MUSJUNDR 




(SEO 


ID 


NO: 


Rat 












C: 


RAT JUNAP 


-GT-T G--* 


(SEO 


ID 


NO: 




RATRJC9 


-GT-T- G — 


(SEQ 


ID 


NO: 



Chicken 












C: CHKJUN 


-GT 


- -A- A 


(SEQ 


ID 


NO 


Ouai I 












C: OULJUN 






(SEQ 


ID 


NO 


Drosophi La 












C: DROJUN 


G--A--T-- 




(SEQ 


ID 


NO 


Homologous sequences 


in GenBank 











CELPOU1 

ECOPR1AY 

SINOCK82 

TRPPROC 

ECOCYSJIHA 



C.elegens RNA 
E.col i pr i mo . 
Ockelbo 82 g. 
Treponema pa. 
E.COli NADPH- 



• A- AA 



• ACA- 

C 

G 



-TGC 
-TGC 
-G-G 



(SEO ID NO 

(SEQ ID NO 

( SEQ ID NO 

(SEQ ID NO 

(SEQ ID NO 



100) 
101 ) 
1C2) 

103) 
104) 
105) 
106) 
107) 
108) 
109) 

110) 
111) 

112) 

113) 

1U) 

115) 
116) 
117) 
118) 
119) 
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Table 24. House jun-D specific probe (HUSD1063). 



5 


Subtype 


GenSank name 


MUSD1063.0NA 
CCCTCAAAAGCCAGAACACCG 












Hunan 
















9: 


KUMJUNCA 


-G GGC-G G-G- 


(SEQ 


ID 


NO: 


i2u> 




C: 


HUMJUNA 


AA--GC-C GC 


(SEQ 


10 


NO: 


121) 




0: 


HUHJUNDR 


G--T G- 


(SEQ 


ID 


NO: 


122) 


10 


















Mouse 
















8: 


HUSJUNBA 


TA TCC TCTG C 


(SEQ 


ID 


NO: 


123) 




C: 


MUSJUNC 


AA--GC-T GC 


(SEQ 


ID 


NO: 


124) 






MUSJUN 


AA--GC-T GC 


(SEO 


ID 


NO: 


125) 


15 




MUSCJUN 


AA- - GC- T GC 


(SEO 


ID 


NO: 


126) 




D : 


HUSJUNDA 






t n 


NU, 








MUSJUNDR 




(SEO 


ID 


NO: 


128) 






MUSJUNO 




(SEO 


ID 


NO: 


129) 




Rat 














20 


C: 


RATRJG9 


G GCC — GC-C — n 


(SEQ 


ID 


NO: 


130) 






RATJUNAP 


G GCC — GC-C — TT 


(SEQ 


ID 


NO: 


131) 




Chicken 
















C: 


CHtCJUN 


A G-GG- A-A G 


(SEO 


ID 


NO: 


132) 




Quail 














25 


C: 


QULJUN 


C* -G- T TG--G-A 


(SEO 


ID 


NO: 


133) 




Orosophi la 














C: 


DROJUN 


A GGA TGTGGCGC 


(SEQ 


ID 


NO: 


134) 


30 


Homologous sequences in 


GenBank 










H27221 




TG T--T --- 


(SEQ 


ID 


NO: 


135) 




DROAMY 


D.erecta 


G G-A--T --- 


(SEQ 


ID 


NO: 


136) 




0ROAMYQ D.erecta 


G C-A--T 


(SEQ 


ID 


NO: 


137) 




HUMIGCMUDE Iimujnoglob. 


A--C GTT 


(SEO 


ID 


WO: 


138) 



35 
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TaD.e Z5. Human-rodent common jun-B specific probes. 







Name in 
GenBank 


B-5(K (5 ' *3' ) 
CACGACTACAAACTCCTGAAAC 


Seq. ID No: 


B ■ 739 ( 5 1 * 3 1 ) 
GGACAGTACTTTTACCCCCG 


3 Cv} . IV NO . 


3 1 tc 

(bo) 


nOube 
















jun-B 


5 


KUSJUNBA 




60C 




321 


251 








c-JUN 




MUSCJUW 


A- AT • - A- T AT A 


601 




322 




jun-D 




MUSJUND 


T G-C--C-T TTCC 


602 


ACG T-C-C GAA 


323 




Hunan 












^ 




jun-B 




HUM J UN C A 




603 




325 


251 








C- jun 


10 


HUHJUNA 


AGTA- -CC — GA 


604 


A- -C-GCGG AC 


326 




jun-D 




HUHJUNDR 


T G-C--C-T TTTG 


605 


TC T-C-C AA 


327 




GenBank 


i > ■ 


HUSMET ] 
RATMYHCD2 


T---CG GTA 


606 


C AT---A 


328 





*: Size of PGR products using B-504 and B-739 as primers. 

**: The highest matched sequences in GenBank (release 68.0) next to the jun genes. 
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Table 26. Hunan- roaent common c-jun specific prooes. 





Marre i n 


C-2101 C5'-3') 
GAGGAAAAAGTGAAAACCTTGAAAGC 


Seq. ID No: 


C-2219 (5'-3') 
GCCAACTCATGCTAACGCAG 


Seq. ID No: 


SUe 
(bo) 


House 
















jun-B 


5 


HUSJUMBA 




6U 


GT-GC GG-GTC 


329 












615 




330 


138 


C-JUN 




HUSCJUH 






jun-D 




HUSJUMD 




616 


T A-GCC--A-G A 


331 




Human 
















Jun-B 




KUHJUNCA 




617 


G-C-CC-AG GC 


332 






10 






618 




333 


138 


C-jun 


HUMJUNA 






jun-D 




HUHJUNDR 




619 


T G-GCO-A-G A 


334 




GenSank 




















HUMCASTA 




620 












HUSOCT22 








335 





13 



*: Size of PCR praaucts using C-2101 and C-2219 as primers. 

**: The highest matched sequences in GenBank (release 68.0) next to the jun genes. 
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Table 27. Human-rodent common jun-D specific probes. 



Mouse 
* jun-B 5 

C-JUN 
^ jurvD 
Human 
jun-B 
C-jun 10 
jun-D 
GenBank 



Name in 
GenBank 



MUSJUNBA 

MUSCJUN 

MUSJUND 

HUMJUNCA 
HUM J UNA 
HUMJUNDR 

HUHTRKR 
HUNBSSl 



0-916 <5<-3') Seq. ID No: 

GACCTCCCCAGCTTCGG 



CTG-G ACT - - 628 

T-A AG---AGAC 629 

630 

C---C CC C 631 

T-A---TT CG-- 632 

633 

TG 634 



0-1153 (V-3< ) 
AAAGTCCTCAGCCACGTCAAC 



Sec. ID No: 



-G---A-G-C---T G - 336 

A'G-A T--- 337 

338 

•G---A-G-C G- 339 

A "G-A 340 

341 

"A A--CTG 342 



Size of PGR products using 0-916 and 0-1153 as primers. 

The highest matched sequences in GenBank (release 68.0) next to the jun genes. 



Size 
(bo) 



258 



258 
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Identifvj n g Pn)he s PQR Primers 

Probes and PCR primers for use in the methods of the present invention can be 
identified in any way known to the arc. Preferably, however, such probes and prftrreiyare 
identified by a computer. We have developed a novel computer system for identifying the 
5 sequences to be used in such probes and primers. This system is an automated system 
which allows.ihe user to calculate and design extremely accurate oligonucleotide probes 
and PCR primers. 

The software 0 f the present invention runs under Microsoft Windows® on IBM® 
compatible personal computers (PCs). This invention allows a researcher to design 

10 oligonucleotide probes based on the GenBank database of DNA and mRNA sequences. 

The present invention further allows examination of probes for specificity or commonality 
with respect to a user-selected target gene sequences. Hybridization strength between a 
probe and a target subsequence of DNA or mRNA can he estimated through a 
hybridization strength model. Quantitatively, hybridization strength is given as the melting 

15 temperature (Tm). 

Two models for estimating hybridization strength models are supported by this 
invention: I ) the Mismatch Model and 2) the H-Site Model. In either case, the user can 
select the following calculations for each probe, results of which are then made available 
for display and analysis: 1) Sequence, Melting Temperature (Tm) and Hairpin 

20 characteristics (a hairpin is a nucleoside sequence that is homologous to itself and can 

"fold back" with one portion of the probe hybridizing to another portion of the same 
probe): 2) Hybridization to other species within the preparation mixture; and (3) Location 
and Tm for the stropcest hybridizations. The results of the invention's calculations are 
then displayed on a Mitsuhashi Probe Selection Diagram (MPSD) which is a graphic 

25 display of all potential hybridizations between the target mRNA and the probe sequences 

in the preparation. 

The Main OligoProbe Design Station dialog window controls all user-definable 
settings in the program. The user is offered a number of options at this window. The File 
option allows the user to print, print in color, save selected probes, and exit the program. 
30 The Preparation option allows the user to open and create preparation (PRP) files. The 
Models option allows the user to chose between the two hybridization models currently 
supported by the OligoProbe DesignStation: 1) the H-Site Model and 2) the Mismatch 
Model. 



If the user selects the H-Site Model option, the melting temperature tor each probe 

and the nucleation threshold parameters can be set. The nucleation threshold is the 

number of base pairs constituting a nucleation site (a subsequence with an exact match) 

If the user selects the Mismatch Model option, the probe length and mismatches (N) can 

5 be set. 

Mismatch Model 
* — i 

The Mismatch Model is used for designing DNA and mRNA probes utilizing 
sequence database information from sources such as GenBank. In this Model, 
hvbridization strength is related only to the number of base pair mismatches between a 
10 probe and its target. Generally, the more mismatches a user allows when setting 

parameters, the more probes will he identified. The Mismatch Model does not take into 
account the GC content of candidate probes so there is no calculation of the probe s 
binding strength. 

The basic technologies employed by the Mismatch model are hashing and continuous 
15 seed filtration. Hashing involves the application ot an algorithm to the records in a set of 

data to obtain a symmetric grouping of the records. When using an indexed set of data 
such as a database, hashing is the process of transforming a record key to an index value 
for storing and retrieving a record. The Mismatch Model is essentially a quick process for 
determining exact and inexact matching between DNA and mRNA sequences to support 
20 the Mitsuhashi Probe Selection Diagram (MPSD). 

The algorithm used by the Mismatch Model is based on the Waterman-Pevzner 
Algorithm (WPALG), which is a computer-based probe selection process. Essentially, this 
is a combination of n£w and improved pattern matching processes. See Hume and Sunday 
(1991. Ref. 4), Landau euil (1986-1990. Rets. r>. 7, K). Grossi and Luccio (1989, Ref. 3), 
25 and Ukkonen ( 1982. Ref. 14). 

There are three principal programs that make up the Mismatch Model in this 
implementation of the invention. The first is designated by the inventors as "k^diff." 
WPALG uses k diff to find all locations of matches of length greater than or equal to one 

* (1) (length is user-specified) with less than or equal to k number of mismatches (k is also 
30 user-specified) between the two sequences. If a candidate oligonucleotide probe fails to 

* match that well, it is considered unique, k diff uses hashing and continuous seed filtration, 
and looks for homologs by searching GenBank and other databases with similar file 
formats. The technique of continuous seed filtration allows for much more efficient 
searching than previously implemented techniques. 
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A seed is defined in this invention to be a subsequence having a length equal to the 
longest exact match in the worst case scenario. For example, suppose the user selects a 
probe length (I) of 18. with 2 or fewer mismatches (k). It a match exists with 2 
mismatches, then there must be a perfectly matching subsequence of length equal to 6. 

5 Once the seed length has been determined, the Mismatch Model looks at all substrings of 

that seed length (in this example, the seed length would be 6). finds the perfectly matched 
base pair subsequence of length equals 6, and then looks to see if this subsequence extends 
to a sequence of length equal to the user selected probe length (i.e.. 18 in this example). 
If so, a candidate probe has been found that meets the user's criteria. 

10 Where the seed size is large (i.e.. a long string of unique nucleotides), the program 

allocates a relatively large amount of memory for the hash table. This invention has an 
option that allows memory allocation for GenBank entries just once at the beginning of 
the program, instead of reallocating memory for each GenBank entry. This reduces input 
time for GenBank entries by as much as a factor of two (2). but this method requires the 

15 user to know the maximum GenBank entry size in advance. 

A probe is found to hybridize if it has k or fewer mismatches with a target sequence 
from the database or file searched. The hit extension time for all appropriate parameters 
of the Mismatch Model has been found by experimentation to be less than thirty-five (35) 
seconds, except in one case where the minimum probe length (I) was set to 24 and the 

20 maximum number of mismatches (k) was set to four (4). This situation would rarely be 

used in real gene localization experiments because the hybridization conditions are too 
weak. 

H-Site Model 

In this embodiment of the invention, the second hybridization strength model is 
25 termed the H-Site Model. One aspect of the H-Site Model uses a generalization of an 

experimental formula to analize nucleotide binding strength. The basic formula on which 
this aspect of the model is built is as follows: 

Tm = 81.5 - 16.6(log[Na]) - .63 ^(formamide) + .41 (^(G + C)) - 600 / N 
In this formula. log[Na] is the log of the sodium concentration. <£(G + C) is the fraction 
30 of matched base pairs which are G-C complementary, and N is the probe length. This 

formula relates the fact that melting temperature is a Junction of both probe length and 
percent GC content. This basic formula has been modified in this invention to account 
for the presence of mismatches. Each percent of mismatch reduces the melting 
temperature by an average of 1.25° (2°C for an AT mismatch, and 4°C for a GC 
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mismatch). This formula is, however, an approximation. The actual melting temperature 
might potentially differ trom this approximation, especially tor short probes or probes with 
a relatively large number ol mismatches. 

Hybridization strength in the H-Site Model is related to each of the following factors: 
5 1) "binding region": 2) type of mismatch (GC or AT substitution ); 3) length of the probe; 

4) GC content of the binding region; and 5) existence of a Enucleation site" (a subsequence 
with an exact match). The type of mismatch and GC content of the binding region from 
each sequence contributes to a candidate probe's binding strength. The binding strength 
from each probe is thereby determined enabling the user to select an optimal probe. 

10 The fundamental assumption of the H-Site Model is that binding strength is mostly 

determined by a paired subsequence of the probe and target, called the binding region. 
If the subsequence binding region contains more GC pairs than AT pairs, the binding 
strength will be higher due to the greater number of hydrogen bonds between G and C 
bases (three bonds) in comparison to A and T bases (two bonds). Thus. GC rich probes 

15 have a higher melting temperature and subsequently form stronger hybridizations. 

In the H-Site Model the program determines optimal probes, ideally without any 
mismatches to the target gene. With this model, however, a candidate probe can have 
more AT mismatches if the sequence is GC rich. The amount of allowable AT mismatches 
in a specific sequence is determined in the present invention program by looking primarily 

20 at subsequence regions of the probe and target that match without penalizing the probe 

for areas that mismatch. If the mismatches are located at either or both of the ends of 
the binding region, there is little effect on the overall stability of the base-pairing. 
Centrally located mismatches in the binding region are much more deleterious, as this will 
significantly lower the binding strength of the probe. 

25 The formula cited above for the melting temperature applies within the binding 

region. The length of the probe is used to calculate percentages, but all other parameters 
of the formula are applied to the binding region only The H-Site Model further assumes 
the existence ot a nucleation site. The length of this nucleation site may be set by the 
user. Typically, a value of 8 to 10 base pairs is used. To complete the H-Site Model the 

30 binding region is chosen so as to maximize the melting temperature Tm among all regions 

containing a nucleation site, assuming one exists (otherwise. Tm = 0). 

The H-Site Model is more complex than the Mismatch Model discussed above in that 
hybridization strength is modeled as a sum of multiple subsequence contributions, with 
matches generally providing positive binding energy and mismatches generally providinc 
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negative binding energy. The exact binding energies to he used depend only on the 
matched or mismatched pair. These coefficients may be specified by the user, although 
in the current version of this invention these coefficients are not explicitly user-selectable, 
but rather are selectee to best fit the hybridization strength formulas developed by itakura 
5 gt_a!(l!>84. Ref. 5). Bolton and McCarthy ( 1962. Ref. 2). Benner et^l ( 1973. Ref. 1). and 

Southern (1975. Ref. 13). 

A unique aspect of the H-Site Model is that hybridization strength is determined by 
the optimal binding region between the candidate probe and binding locus. This binding 
reeion is called the hybridization site, or h-site. and is selected so as to maximize overall 
10 hybridization strength, so that mismatches outside the binding region do not detract from 

the estimated hybridization strength. Several other unique features of the H-Site Model 
include the fact that it is more oriented toward RNA and especially cDNA sequences than 
DNA sequences, and the tact that the user has control over preparation and environmental 
variables. 

15 Tne ern phasis on RNA and cDNA sequences allows the user to concentrate on coding 

regions of genes, rather than necessitating sorting through all of a genomic sequence for 
the desired probe. The enhanced user control over environmental and preparation 
variables allows the user to more accurately simulate laboratory conditions that closely 
correspond with any experiments he or she is conducting. Further, this implementation 
20 of the invention does some preliminary preprocessing of the GenBank database to sort out 

and select the cDNA sequences. This is done by locating a keyword (in this case CDS) 
in each GenBank record, thereby eliminating any sequences containing introns. 

The Mitsuhashi Probe Selection Diagram (MPSD). Figure 28. is a key feature of this 
invention, as it is a unique way of visualizing the results of the probe designed by the 
25 Mismatch and H-Site Models. It is a graphic display of all of the hybridizations of 

candidate oligonucleotide probes and the target with all sequences in the preparation. 
Given a gene sequence database and a target mRNA sequence, the MPSD graphically 
displays all of the candidate probes and their hybridization strengths with all sequences 
from the database. In the present implementation, each melting temperature is displayed 
30 as a different color, from red (highest Tm) to blue (lowest Tm). The MPSD allows the 
user to see visuallv the number of false hybridizations at various temperatures for all 
candidate probes, and the sources of these false hybridizations (with a loci and sequence 
comparison). A locus may be a specific site or place, or. in the genetic sense, a locus is 
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anv of the homologous parts of a pair of chromosomes that may be occupied by aUelic 
genes. 

These probes may then be used to test for the presence of precursors of specific 
proteins in livinc tissues. The oligonucleotide probes designed with this invention may be 
5 used for medical diagnostic kits, DNA identification, and potentially continuous monitoring 

of metabolic processes in human beings. The present implementation of this computerized 
desicn tool runs under Microsoft® Windows™ v. 3.1 (made by Microsoft Corporation; 
Redmond. Washington) on IBM® compatible personal computers (PC's). 

The H-Site Model of this invention is unique in that it offers a multitude of 
10 information on selected probes and original and distinctive means of visualizing, analyzing 

and selectinc amone candidate probes designed with the invention Candidate probes are 
analyzed usine the H-Site Model for their binding specificity relative to some known set 
of mRNA or DNA sequences, collected in a database such as the GenBank database. The 
first step involves selection of candidate probes at some or all the positions along a given 
15 tarcet. Next, a melting temperature model is selected, and an accounting is made of how 

manv false hybridizations each candidate probe will produce and what the melting 
temperature of each will be. Lastly, the results are presented to the researcher along with 
a unique set of tools for visualizing, analyzing and selecting among the candidate probes. 
This invention is both much faster and much more accurate than the methods that 
20 are currently in use. It is unique because it is the only method that can find not only the 

most specific and unique sequence, but also the common sequences. Further, it allows the 
user to perform many types of analysis on the candidate probes, in addition to comparing 
those probes in various ways to the target sequences and to each other. 

Therefore, it is the object of this invention to provide a practical and user-friendly 
25 system that allows a researcher to design both specific and common oligonucleotide 

probes, and to do this in less time and with much more accuracy than currently done. 

This invention is employed in the form best seen in Figure 20. There, the 
combination of this invention consists ot an IBM® compatible personal computer (PC), 
running software specific to this invention, and having access to a distributed database with 
30 the file formats found in the GenBank database and other related databases. 

Unless defined otherwise, all technical and scientific terms used herein have the same 
meaninc as commonly understood by one of ordinary skill in the art to which this invention 
belongs. Although any methods and materials similar or equivalent to those described 
herein can be used in the practice or testing of the present invention, the preferred 
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methods and materials arc now describee. All publications mentioned hereunder are 
incorporated herein by reference. 

The preferred computer hardware capable of operating this invention involves of a 
system with at least the following specifications (Figure 20): 1 ) an IBM® compatible PC. 

5 generally designated 1A. IB. and 1C. with an 80486 coprocessor, running at 33 Mhz or 

faster: 2) 8 or more MB of RAM. 1A: 3) a hard disk IB with at least 200 MB ot storage 
space, but preferably I GB; 4) a VGA color monitor (1C) with graphics capabilities of a 
size sufficient to display the invention s output in readable format, preferably with a 
resolution of 1024 x 768: and 5) a 580 MB CD ROM drive 5 ( IB of Figure 20 generally 

10 refers to the internal storage systems included in this PC. clockwise from upper right, two 

floppy drives, and a hard disk). Because the software of this invention preferably has a 
Microsoft® Windows™ interface, the user will also need a mouse 2. or some other type of 
pointing device. 

The preferred embodiment of this invention would also include a laser printer 3 
15 and/or a color plotter 4. The invention may also require a modem (which can be internal 

or external) if the user does not have access to the CD ROM versions of the GenBank 
database 8 (containing a variable number of gene sequences 6). If a modem is used, 
information and instructions are transmitted via telephone lines to and from the GenBank 
database 8. If a CD ROM drive 5 is used, the GenBank database (or specific portions of 
20 it) is stored on a number of CDs. 

The computer system should preferably have at least the Microsoft® DOS version 5.0 
operating system running Microsoft® Windows™ version 3.1. All of the programs in the 
preferred embodiment ot the invention were written in the Borland® C++ (Borland 
International. Inc.: Scotts Valley. CA) computer language. It should be noted that 
25 subsequently developed computers, storage systems, and languages may be adapted to 

utilize this invention and vice versa. 

This inventive computer program is designed to enable the user to access DNA. 
mRNA and cDNA sequences stored either in the GenBank or in databases with similar 
file formats. GenBank is a distributed flat file database made up of records, each record 
30 containing a variable number of fields in ASCII file format. The stored database itself is 

distributed, and there is no one database management system (DBMS) common to even 
a majority of its users. One general format, called the line type format, is used both for 
the distributed datahase and for all of GenBanks internal record keeping. All data and 
system files and indexes for GenBank are kept in text files in this line type format. 
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The primary GenBank database is currently distributed in ;i multitude of files or 
divisions, each of which represents the genome of a particular species (or at least as much 
of it as is currently known and sequenced and publicly available). The GenBank provides 
a collection of nucleotide sequences as well as relevant bibliographic and biological 
5 annotation. Release 72.0 (6/92) of the GenBank CD distribution contains over 7 1.000 loci 

with a total of over ninety-two (92) million nucleotides. GenBank is distributed by 
IntelliGenetics, of Mountain View. CA. in cooperation with the National Center for 
Biotechnology Information. National Library of Medecinge. in Bethesda, MD. 

10 1. Overall Description of th e OlicoProbe DesignStation 

a General Theory 

The intent of this invention is to provide one or more fast processes for performing 
exact and inexact matching between DNA sequences to support the Mitsuhashi Probe 
Selection Diagram (MPSD) discussed below, and other analysis with interactive graphical 
15 analysis tools. Hybridization strength between a candidate oligonucleotide probe and a 

target subsequence of DNA. mRNA or cDNA can be estimated through a hybridization 
strencth model. Quantitatively, hybridization strength is given as the melting temperature 
(Tm). Currentlv. two hybridization strength models are supported by the invention: 1) 
the Mismatch Model and 2) the H-Site Model. 
20 b. Inputs 

i Main Oli'joProhe DesiiznStation Dialog Window 

The Main OlicoProbe DesignStation dialog window, Figure 21. controls all user- 
definable settings. This window has a menu bar offering five options: 1) File 10; 2) 
Preparation SO; 3) Models 70; 4) Experiment 90; and 5) Help 50. The File 10 option 
25 allows the user to print, print in color,. save selected probes, and exit the program. The 

Preparation SO option allows the user to open and create preparation (PRP) files. 

The Models 70 option allows the user to chose between the two hybridization models 
currentlv supported by the OlicoProbe DesignStation: 1) the H-Site Model 71 and 2) the 
Mismatch Model 75. If the user selects the H-Site Model 71 option, the left hand menu 
30 of Figure 21C is displayed and the user sets the following model parameters: 1) the 

melting temperature Tm 72 for which probes are being designed (i.e.. the melting 
temperature that corresponds to a particular experiment or condition the user desires to 
simulate); and 2) the nucleation threshold 7.V which is the number of base pairs 
constituting a nucleation site. If the user selects the Mismatch Model 75 option, the right 
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hand menu of Figure 21C is displayed and the user sets the following model parameters: 
1) probe length 76. which is the number of base pairs in probes to be considered: and 2) 
mismatch N 77. which is the maximum number ot mismatches constituting a hybridization. 
Computation ot the user s request takes longer with the H-Siie Model if the threshold 73 
settinc is decreased, but longer with the Mismatch Model if the number of mismatches K 
77 is increased. 

In addition, for hoth Model options the user chooses the target species 11 DNA or 
mRNA for which probes are being designed and the preparation 12. a file of all sequences 
with which hybridizations are to be calculated. A sample of a target species file is shown 
in Figure 47 (humbjunx.cds). while a sample of a preparation file is shown in Figure 48 
(junmiwseq). Each of these inputs is represented by a file name and extension in standard 
DOS format. In the target species and preparation fields, the file format follows the 
GenBank format with each ot the fields having a default file extension. Pressing the "OK" 
button 91 (Figure 21C) will initiate processing, while pressing the "Cancel" button 93 will 

15 stop the processing. 

The Experiment 90 option and the Help 50 option are expansion options not yet 
available in the current implementation of the invention, 
c. Processing 

Ficure 22 is a flow chart of the overall OligoProbe DesignStation Program, illustrating 
its sequence and structure. Generally, the main or "control" program of the OligoProbe 
DesignStation performs overall maintenance and control functions. This program, as 
illustrated in Figure 22. accomplishes the general housekeeping functions 51. such as 
defining global variables. The user-friendly interlace 53. carries out the user-input 
procedures 55. the file 57 or database 59 access procedures, calling of the model program 
25 62 or 63 selected by the user, and the user-selected report 65 or display 67. 69. 71 and 73 

features. Each of these features is discussed in more detail in later sections, with the 
exception of the input procedures, which involves capturing the user's set-up and control 
inputs. 

d. Outputs 

30 l The Mitsuhashi Probe Selection Diagram Window 

The Mitsuhashi Probe Selection Diagram (MPSD). Figure 23. is a key feature of the 
invention as it is a unique way of visualizing the results of the program's calculations. It 
.is a graphic display of all ot the hybridizations of probes with the target oligonucleotides 
in the preparation. Specifically, given a nucleotide sequence database and a target mRNA. 
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the MPSD graphically displays ail of the candidate probes iind their hybridization strengths 
with all sequences from the nucleotide database. The MPSD allows the user to visualize 
the number of false hybridizations at various temperatures for all candidate probes, and 
the sources of these talse hybridizations (with a loci and sequence comparison). 
5 For each melting temperature selected, a graph showing the number of hybridizations 

for each probe is displayed. In the preferred embodiment, the graphs are color coded. 
In this implementation of the invention, the color red 123 identifies the highest melting 
temperature and the color blue 124 identifies the lowest melting temperature. Each 
mismatch results in a reduction in the Tm value. The melting temperature is also a 
10 function of probe length and percent GC content. Within the window, the cursor 125 

shape is changed from a vertical line bisecting the screen to a small rectangle when the 
user selects a particular probe. The current probe is defined to be that probe under the 
cursor position (whether it be a line or a rectangle) in the MPSD window. More detailed 
information about the current probe is given in the Probelnfo and Matchlnfo windows. 
15 discussed below Clicking the mouse button 2 once at the cursor 125 selects the current 

probe. Clicking the mouse button 2 a second time deselects the current probe. Moving 
the cursor across the screen causes the display to change and reflect the candidate probe 
under the current cursor position. 

The x-axis 1 10 of the MPSD. Figure 23. shows the candidate probes' starting positions 
20 along the given mRNA sequence. The user may "slide" the display to the left or right in 

order to display other probe starting positions. The y-axis 1 15 of the MPSD displays the 
probe specificity, which is calculated by the program. 

The menu options 116, 117. 118. 119. and 120 available to the user while in the 
MPSD. Figure 23. and are displayed along a menu bar at the top of the screen. The user 
25 can click the mouse 2 on the preferred .option to briefly display the option choices, or can 

click and hold the mouse button on the option to allow an option to be selected. The user 
may also type a combination of keystrokes in order to display an option in accordance with 
well-known computer desk top interface operations. This combination usually involves 
holding down the ALT key while pressing the key representing the first letter of the 
30 desired option (i.e. F. P, M. Eor H). 

The File option 116 allows the user to specify input files and databases. The 
Preparation option 117 allows the user to create a preparation file summarizing the 
sequence database. The Models option 1 IS allows the user to specify the hybridization 
model (i.e.. H-Site or Mismatch) and its parameters The Experiment option 1 \9 and the 
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Help option 120 are not available in the current implementation of this invention. These 
options are part of the original Main OligoProbe DesignStation dialog window. Figure 21. 

Areas on the graphical display of the MPSD. Figure 11. where the hybridizations tor 
the optimal probes are displayed are lowest and most similar, such as shown at 121. 
indicate that the particular sequence displayed is common to all sequences. Areas on the 
graphical display of the MPSD where the hybridizations for the optimal probes are 
displaved are highest and most dissimilar, such as shown at 122. indicate that the particular 
sequence displayed is extremely specific to that particular gene fragment. The high points 
on the MPSD show many loci in the database, to which the candidate probe will hybridize 
(i.e., manv false hybridizations). The low points show few hybridizations, at least relative 
to the riven database. Specifically, the sequence shown at 121 would reflect a probe 
common to all of the gene fragments tested, such that this probe could be used to detect 
each of these genes. The sequence shown at 122 would reflect a probe specific to the 
particular genJ fragment. «>ch that this probe could be used to detect this particular gene 
15 and no others. 

ii The Prohelnto and M-'itehlnfo Window 

The combined Probelnfo and Matchlnfo Window. Figure 24. displays detailed 
information about the current candidate probe. The upper portion of the window is the 
Probelnfo window, and the lower portion is the Matchlnfo window. The Probelnfo 
window portion displays the following types of information: the target locus (i.e.. the 
mRNA. cDNA. or DNA from which the user is looking for probes) is displayed at 131. 
while the preparation used for hybridizations is displayed at 132. In the example shown 
in Figure 24. the target locus 131 is the file named HUMBJUNX.CDS. which is shown as 
being located on drive F in the subdirectory MILAN. The preparation 132 is shown as 
being the file designated JUNMIX.PRP. which is also shown as being located on drive F 
in the subdirectory MILAN. The JUNMIX.PRP preparation in this example is a mixture 
of human and mouse jun loci. 

The current and optimal probe s starting position is shown at 135. The current 
candidate oligonucleotide probe is defined at 136. and is listed at 137 as having a length 
of 21 bases. The melting temperature for the probe 13(> as hybridized with the targets is 
shown in column 140. The melting temperature for the optimal prohe is given as ft 1.7 
degrees C at 138. The Probelnfo Window Figure 24 also displays hairpin characteristics 
of the probe at 139. In the example shown, the Probelnfo Window shows that there are 
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four (4) base pairs involved in the worst hairpin, and that the worst hairpin has a length 
of one (1) (see Figure 24. at 139) 

The Maichlnto Window portion displays a list of hybridizations between the current 
probe and species within the preparation file, including hybridization loci and hybridization 
5 temperatures. The hybridizations are listed in descending order by melting temperature. 

The display shows the locus with which the hybridization occurs, the position within the 
locus, and the hybridization sequence. 

In the Matchlnfo window portion, the candidate probe 13(> is shown at 150 as 
hybridizing completely with a high binding strength. This is because the target DNA is 
10 itself represented in the database in this case, so the candidate probe is seen at 150 to 

hybridize with itself (a perfect hybridization). The locus of each hybridization from the 
preparation 132 are displayed in column 141. while the starting position of each 
hybridization is given in column 142. The calculated hybridizations are shown at 145. 
iii. The ProhesEdit Window 

15 The ProbesEdit Window. Figure 25, is a text editing window provided for convenient 

editing and annotation of OligoProbe DesignStation text file output. It is also used to 
accumulate probes selected from the MPSD, Figure 23. by mouse button 2 clicks. 
Standard text editing capabilities are available within the ProbesEdit Window. The user 
may accumulate selected probes in this window (see 155 for an example) and then save 
20 them to a file (which will bear the name of the preparation sequence with the file 

extension of "prh" 15(>. or may be another file name selected by the user). A sample of 
this file is shown in Figure 26A. 

iv. Miscellaneous Output 

The present embodiment of this invention also creates two output files, currently 
25 named "test.out" and "test 1. out", depending upon which model the user has selected. The 

first file, "test.out". is created with both the Mismatch Model and the H-Site Model This 
file is a textual representation of the Mitsuhashi Probe Selection Diagram (MPSD). It 
breaks the probe sequence down by position, length, delta Tm. screensN. and the actual 
probe sequence (i.e.. nucleotides). An example of this file created by the Mismatch Model 
30 is shown in Figure 40, and example created by the H-Site Model is shown in Figure 44a. 

The second file, "test 1. out", is created only by the H-Site Model. This file is a textual 
representation of the Probelnfo and Matchlnfo window that captures all hybridizations, 
alone with their locus, starting position, melting temperature, and possible other 
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hybridizations. A partial example of this file is shown in Figure 44b ( 10 paces out of a 
total of 100 pages created by the H-Site Model). 

2. Description of the Mism atch Model Program 

5 a Overview 

In this invention, one of the hybridization strength models is termed the Mismatch 
Model (see Figure 21 for selection of this model). The basic operation of this model 
involves the techniques of hashing and continuous seed filtration, as defined earlier but 
described in more detail below. The essence of the Mismatch Model is a fast process for 

10 doing exact and inexact matching between nucleotide sequences to support the Mitsuhashi 

Probe Selection Diagram (MPSD). There are a number of modules in the present 
implementation of the Mismatch Model contained in this invention, the most significant 
of which are shown in the flow chart in Figure 27 and in more detail in Figures 28 through 
38. The main kdiff module shown in the flow chart in Figure 28 is a structured program 

15 that provides overall control of the Mismatch Model, calling various submodules that 

perform different functions. 

b. Inputs 

The user-selected input variables for this model are minimum probe length 76 (which 
is generally from 18 to 30) and maximum number of mismatches 77 (which generally is 
20 from 1 to 5). These inputs are entered by the user in the Main OligoProbe DesignStation 

Dialog Window. Figure 21C. 

c. Processing 

i. kdiff Program 

Some terms of art need to be defined before the processing performed by this module 
25 can be explained. A hash table basicajly is an array or table ot data. A linked list is a 

classical data structure which is a chain, of linked entries and involves pointers to other 
entry structures. Entries in a linked list do not have to be stored sequentially in memory, 
as is the case with elements contained in an array. Usually there is a pointer to the list 
associated with the list, which is often initially set to point to the start of the list. A 
30 pointer to a list is useful for sequencing through the entries in the list. A null pointer (i.e.. 

a pointer with a value of zero) is used to mark the end of the list. 

As the flow charts in Figures 27 and 28 illustrate, the general process steps and 
implemented functions of this model can be outlined as follows: 
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Siep 1: First, create a hash lahie and linked list trom the query (Figure r. nashing 
module 222). 

Step 2: Next, while there are stiti GenBank entries available tor searching (Figure 

27. assembly mouuie 250): 
5 Step 2a: Read the current GenBank entry (record) sequence of user- 

specified length (Figure 27, seqioad module 232), or read the current 
sequence (record) from the file selected by the user (Figure 27. readl 
module 234). 

Step 2b: For the current sequence for each position of the sequence 
10 from the first position (or nucleotide) to the last position (or nucleotide) 

(incrementing the position number once each iteration of the loop) 
(Fieure 27. q_colour module 242). 

Siep 2c: set the variable dna_hash equal to the hash of the current 
position of the current sequence (Figure 27. q_colour module 242). 
15 Step 2d: While not at the end of the linked list for dna_hash (Figure 

27. q_colour module 242). 

Step 2e: set the query pos equal to the current position of dna_hash in 
the linked list (Figure 27. q_colour module 242) and 
Step 2f: Extend the hit with the coordinates (query_pos. dna_pos) 
20 (Figure 27. hit_ext module 244), 

Step 2g: If there exists a kjnismatch in the current extended hit 
(Ficure 27. colour module 246). then 

Step 2h: print the current hit (Figure 27. q_colour module 242), and 
repeat from Step 2. 

25 As this illustrates, there are three (3) basic looping or iteration processes with 

functions being performed based on variables such as whether the GenBank section end 
has been reached (the first "WHILE" loop. Step 2). whether the end of the current DNA 
entry has been reached (the "FOR" loop. Step 2b). and whether the end of the dna_hash 
linked list has been reached (the second "WHILE" loop. Step 2d). A "hit" will only be 

30 printed if there are k_mismatches in the current extended hit. 

Figures 28 through 38 illustrate the functions of each of the modules of the present 
embodiment of this invention, all of which were generalized and summarized in the 
description above. Figure 28. which outlines the main "k_diff module, shows that this 
module is primarily a program organization and direction module, in addition to 
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performing routine "housekeeping" functions, such as defining the variables and hash tables 
251. checking .1 the user-selected gene sequence tile is open 252. extracting needed 
identification information from the GenBank 253. and ensuring valid user input 254. This 
module also performs a one-time allocation ot memory for the gene sequences, and 
allocates memory for hit information, hashing, hybridization and frequency length profiles 
and output displays. 255 & 256. The -k__difT module also initializes or "zeros out" the 
hashing table, the linked hashing list and the various other variahles 257 in preparation for 
the hashing function. In addition, this module forms the hash tables 258 and extracts a 
sequence and finds the sequence length 259. 

One of the most important functions performed by the "k_diff module is to define 
the seed (or kernel or kjuple) size. This is done by setting the variable kjuple equal to 
(min_probe_length-max-_mismatch_#)/(max_mismatch + # + 1) Figure 28 at 265. Next, 
if the remainder of the aforementioned process is not equal to zero 266. then the value 
of the variable kjuple is incremented by one 267. The resulting value is the size of the 
15 seed. The module then reads the query 268 and copies the LOCUS name 269 for 

identification purposes (a definition of the term locus is given earlier in the specification). 

The "k_diff" module Figure 28 also calls the "assembly" module 260. writes the results 
to a file 261a. plots the results 26 lb (discussed below), calculates the hairpin characteristics 
262 (i.e.. the number of base pairs and the length of the worst.hairpin) and the melting 
temperature (Tm) for each candidate probe 263. and saxes the results to a file 264. 

The screen graphs are plotted 261b by converting the result values to pixels, filing a 
pixel array and performing a binary search into the pixel array. Next, given the number 
of pixels per probe position and which function is ot interest to the user (i.e.. the three 
mismatch match numbers), the program interpolates the values at the value of 
25 (pixelsPerPositionN-l) and computes .the array of pixel values for drawing the graph. 

These values are then plotted on the MPSD. 

The "hashing" module. Figure 29. performs hashing ot the -query. In other words, it 
creates the hash table and linked list of query positions with the same hash. The variable 
has jable[i] equals the position of the first occurrence of hash i in the query. If i does not 
30 appear in the query. hash_table[i] is set to zero. 

The "tran" module. Figure 30. is called by the "hashing" module 271. and performs the 
hashing of the sequence of kjuple (kernel or seed) size. If the kjuple exists (i.e.. its 
length is greater than zero), the variable uns is set equal to uns ALF+p 291. The variable 
p represents the digit returned by the "let jlig" module Figure 3 1 that represents the 
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nucleotide being examined. ALF is a constant that is set by the program in this 
implementation to equal four. The query pointer is then incremented, white the size of 
k tuple (the seed) is decremented 292. This process is repeated until the sequence of 
k tuple has been entirely ha>hed. Then ihe t: tran module returns the variable 

5 current_hash 293 to the "hashing" module Figure 29. 

The "let dig" module. Figure 31, is called by the "tran" module 291, and transforms 
the nucleotides represented as the characters "A". T'. "U". "G" and M C in the GenBank 
and the user s query into numeric digits for easier processing by the program. This 
module transforms "a" and "A" into "0" 301. V. "T\ V and W U" into "1" 302. V and "G" 

10 into "2" 303, and "c" and "C" into "3 M 305. If the character to be transformed does not 

match anv one of those listed above, the module returns M -l" 305. The "hashing" module. 
Figure 29. then calls the "update" module 272. Figure 32. which updates the hash with a 
sliding window (i.e.. it forms a new hash after shifting the old hash by T). The remainder 
of old_hash divided by power_l is calculated 311 (a modulus operation), the remainder is 

15 multiplied by ALF 312 (i.e.. four), and then the digit representing the nucleotide is added 

to the result 313. The "update" module then returns the result 314 to the "hashing" module 
Figure 29. 

If the current hash has already occurred in the query, the program searches for the 
end of the linked list for the current hash 273 and marks the end of the linked list for the 

20 current hash 274. If the current hash has not already occurred in the query, the program 

puts the hash into the hash table 275. The resulting hash table and linked list are then 
returned to the "kjJiff module. Figure 28 at 258. 

The "assembly" module. Figure 33. extracts sequences from the GenBank and 
performs hit locating and extending functions. This module is called by the "k_diff" module 

25 Figure 28 at 260 if the user has chosen to use the database to locate matches. The output 

from the "assembly" module (Figure 33,) tells the user that the section of the database 
searched contains E number of entries 32 1 of S summary length 322 with H number of hits 

323. Further, the program tells the user that the number of considered l-tuples equals T 

324. The entry head line is also printed 326. 

30 The "seqload" module. Figure 34, is called by the "k_diff module Figure 28 at 259 

once the query hash table and linked list have been formed by the "hashing 1 module Figure 
29. The "seqload" module Figure 34 checks to see if the end of the GenBank file has been 
reached 327. and, if not. searches until a record is tound with LOCUS in the head-line 328. 



WO 93/15221 



PCT/US93/00977 



-94- 



Next, the LOCUS name is extracted 329 for identification purposes, and the program 
searches for the ORIGIN field in the record 331). 

The program then extracts the current sequence 33 1 from the GenBank and performs 
two passes on each sequence. The first is to determine the sequence length 332 and 
5 allocate memory for each sequence 333. and the second pass is to read the sequence into 

the allocated memory 334. Since the sequences being extracted can contain either DNA 
nucleotides or protein nucleotides, the "seqload" module can recognize the characters "A", 
T\ "U '. "G". and "C. The bases "A". "T". "G" and "C are used in DNA sequences, while 
the bases "A". "U". "G" and "C" are used in RNA and mRNA sequences. The extracted 
10 sequence is then positioned according to the type of nucleotides contained in the sequence 

335. and the process is repeated. Once the end of the sequence has been reached, the 
"seqload" module returns the sequence length 33<> to the "k_diff" module Figure 28. 

If the user has chosen to use one or more files to locate matches, rather than the 
database, the "read 1" module. Figure 35. rather than the "seqload" module Figure 34. is 
15 called by the "k_diff" module Figure 28. The "read I" module. Figure 35. reads the 

sequence from the user specified query file 341 and allocates memory 342. This module 
also determines the query length 343. extracts sequence identification information 344. 
determines the sequence length 345. transforms each nucleotide into a digit 346 by calling 
the "let dig" module Figure 31. creates the query hash table 347 by calling the "digjet: 
20 module Figure 36. and closes the file 348 once everything has been read in. 

First, the "readl" module Figure 35 allocates space for the query 342. To do this, the 
"ckalloc" module. Figure 35 at 342. is called. This module allocates space and checks 
whether this allocation is successful (i.e.. is there enough memory or has the program run 
out of memory). After allocating space, the "readl" module Figure 35 opens the user- 
25 specified file 34') (the "ckopen" module. Figure 35 at 349. is called to ensure that the query 

Tile can be successfully opened 349). determines the query length 343. locates a record with 
LOCUS in the head-line and extracts the LOCUS name 344 for identification purposes, 
locates the ORIGIN Held in the record and then reads the query sequence from the file 
341. Next, the sequence length is determined 345. memory is allocated for the sequence 
30 342. and the sequence is read into the query file 350. If the string has previously been 

found, processing is returned to 344. If not. then each character in the query file is read 
into memory 350. 

The characters are transformed into digits 34ft using the "let_dig" module. Figure 31. 
until a valid dinit has been found, and then the hash table containing the query is set up 
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347 using tnc module "digjet". Figure 30, which transforms the digits into nucleotides 
represented by the characters "A" 371. T' 371. "G" 373. M C 374, and "X" 375 as a defauk. 
If the end ol the tile has not been reached, processing is returned to 344. If it has. the file 
is dosed 348 and the query is then returned to the "readl" module Figure 35 at 347. 

The "q_colour" module. Figure 27 (Figure 33 at 325). is called by the "assembly" 
module Figure 33 after the current sequence has been extracted from the GenBank. The 
"q_colour" module Figure 37 performs the heart of the Mismatch Model process in that 
it performs the comparison between the query and the database or file sequences If the 
module finds that there exists a long (i.e., greater than the minjiitjength) extended hit, 
it returns a "1" to the "assembly" module Figure 34. Otherwise, the "q_colour" module. 
Figure 37, returns a "0". 

In the "q_colour" moduie. Figure 37, all DNA positions are analyzed in the following 
manner. First, the entire DNA sequence is analyzed 391 to see whether each position is 
equal to zero 392 (i.e., whether it is empty or the sequence is finished). If it is not equal 
to zero 393, the "q_colour" module Figure 40 calls the "tran" module. Figure 30 described 
above, which performs the hashing of k_tuples. The "Iran" module Figure 30 calls other 
modules which transform the nucleotides represented by characters into digits for easier 
processing by the program and then updates the hash with a sliding window. If the 
position is equal to zero, the current_hash position is set to new_has after one shift of 
old_hash 390 by calling the "update" module Figure 32. 

If the nucleotide at the current_hash position is equal to zero, processing is returned 
to 391. If not. the query position is set equal to (nucleotide at current hash position - 1). 
Next, the "qjrolour" module Figure 37 looks for the current_hash in the hash table 394. 
If the current k tuple does not match the query 395. then the next k_tuple is considered 
395. and processing is returned to 391. If the current k_tuple does match the query, then 
the program checks the hit's (i.e.. the match's) vicinity 396 hy calling the "hit_e\*t H module. 
Figure 38 to determine if the hit is weak. The inventors have found that if the code for 
the module "hit ext" is included within the module "q colour", rather than being a separate 
module utilizing the parameter transfer machinery. 25*7 of CPU time can be saved. 

The "hit_ext" module Figure 38 determines the current query position in the hit's 
vicinity 42 1. determines the current DNA position in the hit's vicinity 422. and creates the 
list of mismatch positions (i.e.. the mismatch_Iocation_ahead 423. the 
mismatch locaiion_behind 423 and the kernel match location). If the hit is weak 424. the 
"hit ext" module Figure 3S returns "0" to the "q^colour" module Figure 37. If the hit has 
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a chance ro contain 425. the module returns "1" to the -q_.-.-mur module Figure 37. A hit 
has a chance to contain, and is therefore no: considered weak, if the 
mismatch locationjihead - the mismatch Jocation_behind is greater than the 
min hitjength. If not. it is a short hit and is too weak. 

If the "hit ext" module Figure 38 tells the "q_colour" module Figure 37 that the hit 
was not a weak one, then the M q_colour B module determines whether the current hit is long 
enouch 398 hv calling the "colour" module Figure 39. The "colour" module Figure 39 
performs queryjrolour modification by the hit data, starting at pos_querv and described 
by mismatch Jocation_ahead and mismatch Jocation_behind. After the variables to be 
used in this module are defined, variable isw_print (which is the switch indicating the hit 
leneth) is initialized to zero 430. The curjength is then set equal to the length of the 
extending hit 431 ( mismatch Jocation_behind[i] + mismatch Jocitionjthead[j]- 1). Next, 
if cur leneth is greater than or equal to the min_hitjength 432 (i.e.. the minimum 
considered probe size), the hit is considered long and isw_print is set equal to two 433. 
15 The value of is\v__print is then returned 434 to the "q_cnIour" module Figure 37. 

If the leneth of the extending hit is longer than the min_hit_length, the hit is 
considered lone 399. Otherwise, the hit is considered short. If the hit is short, nothing 
more is done to the current hit and the module begins again. If. on the other hand, the 
hit is considered long 399, the "q_colour M module Figure 37 prints the current extended hit 
20 400. The current extended hit can be printed in ASCII, printed in a binary file, or printed 

to a memory file. The "q_colour" module Figure 37 then repeats until the end of the 
linked list is reached. 

d. Outputs 

The output of the k_diff program may be either a binary file containing the number 
25 of extended hits and the k_mismatch hit locations (see Figure 40). or the output may be 

kept in memory without writing it to a file. See Section I(d)(iv) for more detail. 
3. Description of the H-Site Model Program 
a. Overview 

In this invention, the second hybridization strength model is termed the H-Site Model 
30 (see Ficure 21 for user selection of this model). The formula used in the H-Site Model 

is an expression of the fact that melting temperature Tm is a function of both probe length 
and percent of GC content. This basic formula has been modified in this invention to 
account for the presence of mismatches. Each percent of mismatch reduces the melting 
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temperature Tm by an average ot J .25 degrees (2 degrees C for ;in AT mismatch, and 4 
degrees C for a GC mismatch). 

In addition, this implementation of the invention does some preliminary preprocessing 
of the GenBank database to sort out and select the cDNA sequences. This is done by 
locating a keyword (in this case CDS) in each GenBank record. 

There are a number of modules in the present embodiment of the H-Site Model 
contained in this invention. Each step of the processing involved in the H-Site Model is 
more fully explained below, and is accompanied by detailed flow charts, 
b. Inputs 

There are two basic user-selected inputs lor the H-Site Model (see Figure 21C): 1) 
the melting temperature Tm 22 for which probes are being designed (i.e.. the melting 
temperature that corresponds to a particular experiment or condition the user desires to 
simulate); and 2) the nudeation threshold 23. which is the number of base pairs 
constituting a nudeation site. The user is also required to select the 1) target species 11 
15 gene sequence(s) (DNA. mRNA or cDNA) for which probes are being designed; 2) the 

preparation 12 of all sequences with which hybridizations are to be calculated; and 3) the 
probe output file 13. The preparation file is the most important, as discussed below, 
c. Organization of the H-Site Model Program 

The current implementation of the H-Site Model program of this invention is 
20 distributed between five files containing numerous modules. The main file is designated 

by the inventors as B ds.cpp M in its uncompiled version. This file provides overall control 
to the entire OligoProbe DesignStation invention. It is divided into six sections. Section 
0 defines and manipulates global variables. Section i controls general variable definition 
and initialization (including the arrays and memory blocks). It also reads and writes 
25 buffers for user input selections, and constructs multi buffers. 

Section 2 sets up and initializes various "snippet" variables (see section below for a 
complete definition of the term snippet), converts base pair characters to a representation 
that is 96 base pairs long and to ASCII base pair strings, and performs other sequence file 
manipulation such as comparing snippets. This section also reads the sequence format file, 
reads base pairs, checks tor and extracts sequence identification information (such as 
ORIGIN and LOCUS) and filters out sequences beginning with numbers. 

Section 3 involves preparation file manipulation. This section performs the 
preprocessing on the PRP file discussed above. It also merges and sorts the snippet files. 
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creates a PRP file and sorts it. and outputs the sorted snippets. Next, this section streams 
through the PRP file. 

Section 4 contains the essential code tor H-Site Model processing (see Figures 41 
through 43 for details, discussed below). Streams are set up. and then RIBI comparisons 
5 are performed tor hybridizations (see tile "ribi.cpp" for definitions of RIBI search 

techniques). Next, probes are generated, binding strength is converted to melting 
temperature, and hybridizations are calculated and stored (including hybridization 
strength). Lastly, other H-Site calculations are performed. 

Section 5 is concerned with formatting and presenting diagnostic and user file 
10 (test.out. test Lout, and test2.out files) output. This section also handles the graphing 

functions (the MPSD diagram in particular). In addition, this section calculates the hairpin 
characteristics tor the H-Site Model candidate probes. 

The second H-Site Model file, designated as "ds.h" defines data variables and 
structures. Section 1 of this file concerns generic data structures (including memory blocks 
15 and arrays, and file inputs ;ind outputs). Section 2 defines the variables and structures 

used with sequences, probes and hybridizations. Section 3 defines variables and structures 
concerned with protocols (i.e.. function prototypes, graphing, etc.). 

The third H-Site Model file, designated as Tuncdoc.txf . contains very detailed 
documentation for this implementation of the H-Site Model program. Numerous variables 
20 and structures are also defined. The flow of the program is clearly shown in this file. 

The fourth H-Site Model file, designated as "ribi.h" handles the sequence 
comparisons. The fifth and last H-Site Model file, designated as 'nbi.cpp", performs 
internal B-Tree indexing. Definitions of Red-black Internal Binary Index (RIBI) searching 
are found in this file. Definitions are also included for the concepts keyed set. index. 
25 binary tree, internal binary index, paths, and red-black trees. Implementation notes are 

also included in this file. 

d. Processing 

Implementation of the H-Site Model in this invention is done in three stages. First, 
the invention creates the preparation (PRP) file, which contains all relevant information 
30 from the sequence database. This is the preprocessing stage discussed above. Next, the 

target is prepared by the program. Lastly, the invention calculates the MPSD data using 
the PRP file and target sequence to find probes. 
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j Creation of The Preprocessed Preparation File 

Ficure 41. Step I: The program first opens the sequence database for reading into 
memory 461. 462. Step 2: Next, as sequence base pairs are read in 462. "snippets" are 
saved to disk 463. along with loci information A snippet is a fixed-length subsequence of 
5 a preparation sequence. The purpose of snippets is to allow the user to examine a small 

portion of a preparation sequence together with its surrounding base pairs. Snippets in 
the implementation of this invention are 96 base pairs long (except for snippets near the 
end or beainning of a sequence, which may have fewer base pairs). The "origin" of the 
snippet is in position 40. For snippets taken near the beginning of a sequence, some of 
10 the initial 40 bases are undefined. For snippets near the end of a sequence, some of the 

final 55 bases are undefined. Snippets are arranged in the preparation file (PRP) in sorted 
order (lexicographical order beginning at position 40). In this invention, the term 
"lexicographical order" means a preselected order, such as alphabetical, numeric or 
alphanumeric. In order to conserve space, snippets are only taken at every 4th position 
15 of the preparation sequence. 

Step 3: The snippets are merge sorted 464 to be able to search quickly for sequences 
which pass the "screen", discussed below. Step A: The merged file is prepended with 
identifiers tor the sources ot the snippets 465. This is done to identify the loci from which 
hybridizations arise. 
20 ii- Tartiet Preparation 

Ficure42. Step 1: The target sequence file is opened 471 and read into memory 472. 
For each position in the target mRNA. the probe defined at that starting position is the 
shortest subsequence starting at that position whose hybridization strength is greater than 
the user specified melting temperature Tm. Typically, the probes are of length 18 to 50. 
25 Step 2: Four lists of "screens" are formed 473. 474. 475. each shifted by one base pair 475 

to correspond to the fact that snippets are only taken at every four base pairs. A screen 
is a subsequence of the target mRNA of length equal to the screening threshold specified 
by the user. The screens are then indexed 476 and sorted in memory 477. 

iii Calculation of the MPSD Data 

30 Fiuure 43. Step 3: This step is the heart of the process. Step 3a: The program 

streams through the following five items in sync, examining them in sequential order: the 
snippet file and the four lists of screens 481-4K4 Step 3b: Each snippet is compared to 
a screen 4K5. Step 3c If the snippet does not match, whichever stream is behind is 
advanced 486 and Step 3b is repealed. If the snippet does match. Step 4 is performed. 
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Step 4: It a snippet and a marching .screen were tound in Srep 3b 487. the 
hybridi2arion strength of the binding between the sequence containing the snippet and all 
of the probes containing the screen is calculated (see Step 5 ). Double counting is avoided 
by doing this. only for the first marched screen containing the probe. Each pair of bases 
5 is examined and assigned a numerical binding strength. An AT pair would be assicned a 

lower binding strength than a GC pair because AT pairs have a lower melting temperature 
Tm. The process is explained more fully below at Srep 5b, 

Step 5: The hybridization strengths between sequence and all the probes containing 
it are calculated using a dynamic programming process. The process is as follows: Step 

10 5a: Begin at the position of the first probe containing the given screen bur not containing 

any orher screens which srarr at an earlier position and also match the sequence. This is 
done to avoid double counting. Two running rorals are maintained; a) boundSrrength, 
which represents rhe hybridization strength contribution which would result if the sequence 
and probe were to match exactly for all base pairs to the right of the current position, and 

15 b) unboundStrength. which represents the strength of the maximally binding region. Step 

5b: At each new base pair, the variable boundStrength is incremented by 71 if the 
sequence and probe match and the matched base pair is GC 489. incremented bv 30 if the 
matched base pair is AT 490 (i.e.. this number is about 42.259?: of the first number 71), 
and decremented by 74.5 if there is not a match 488 (i.e.. this number is about 5% larger 

20 than the first number 71). Step 5c: If the current boundStrength exceeds the current 

unboundStrength 491 (which was originally initialized to zero), a new binding region has 
been found, and unboundStrength is set equal to boundStrength 492. Step 5d: If the 
current boundStrength is negative. boundStrength is reset to zero 493. Step 5e: If the 
current position is at the end of a probe, the results (the hybridization strengths) are 

25 tallied for that probe. Step 5f: If the current position is at the end of the last probe 

containing the screen, the process stops. 

Step 6: A tally is kept of the number and melting temperature of the matches for 
each candidate probe, and the location of the best 20 candidates, using a priority queue 
(reverse order by hybridization strength number) 494. Step 7: A numerical "score" is kept 

30 for each preparation sequence by tallying the quantity exp (which can be expressed as Ee* 

Tm ) for each match 495. where Tm is the melting temperature for the "perfect" match, the 
probe itself. In other words, the probe hybridizes "perfectly" to its tarcet. 

Step 8: Hairpins are calculated by first calculating the complementary probe. In 
other words, the order of the bases in the candidate probe are reversed (CTATAG to 
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GATATC). ;inci complementary base pairs are substituted (A tor T. T for A. G tor C and 
C for G, chancing GATATC to CTATAG in the above example). Next, the variable 
representing the maximum hairpin length for a candidate probe is initialized to zero, as 
is the variable representing a hairpins distance For each offset, the original candidate 
probe and the complementary probe just created are then aligned with each other and 
compared. The longest match is then found. If any two matches have the same lencth, 
the one with the longest hairpin distance (i.e.. the number of base pairs separatinc the 
match) is then saved. 

Step 9: The preparation sequences are then sorted 496 and displayed in rank order, 
from best to worst 497. Step 10: The resulting MPSD, which includes all candidate 
probes, is then displayed on the screen. Step 1 I: The best 20 matches are also printed 
or displayed in rank order, as the user requests 497. 
e. Outputs 

The outputs of the H-Site Model are fully described in Section l(d)(iv), above, and 
illustrated in Figures 23 through 25. Samples of the two output files created by the H-Site 
Model are shown in Figures 44A and 44B. 

4 - Descriptio n of the Mitsubishi Probe Selection Diagram Processing 

Once the Mitsuhashi Probe Selection Diagram (MPSD) data has been calculated by 
the H-Site Model program (see stage three and Figure 43. discussed above), it is necessary 
to convert this data to pixel format and plot a graph. An overview of this process is shown 
in Figure 45. First, the program calculates the output (x.y) ranges 500. Next, these are 
converted to a logarithmic scale 501. The values are then interpolated 502. and a bitmap 
is created 503. Lastly, the bitmap is displayed on the screen 504 in MPSD format 
(discussed above in section l(e)(i)). A sample MPSD is shown in Figure 23. 
25 5 * Description of the Matchlnfo Window Processing 

The Probelnfo and Matchlnfo windows are discussed in great detail in Section 
l(e)(ii), and a sample of these windows is shown in Figure 24. An overview of the 
processing involved in creating the Matchlnfo portion of the window is given in the flow 
chart in Figure 46. First, as the user moves the MPSD cursor 570 (seen as a vertical line 
bisecting the MPSD window), the program updates the position of the candidate probe 
shown under that cursor position 521. Next, based upon the candidate probe's position, 
the program updates the sequence 522 and hairpin information 523 for that probe. This 
updated information is then displayed in an updated match list 524. shown in the 
Matchlnfo window. 
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EfTect of the Invention 

The invented method offers rapid characterization ot mRNA in approximately 5 
hours without the need for mRNA purification. Because conventional methods of mRNA 
quantification, for example Northern blots, require over 70 hours to complete, the method 
5 of the present invention significantly reduces the length of time needed for mRNA 

analysis. Furthermore, because various iypes of first nucleotide probes can be immobilized 
in each well of microliter plate and identified using a single kind of labeled second 
nucleotide probe, the present invention advantageously can provide the simultaneous 
quantification of multiple varieties of mRNA. 

10 As shown in the examples above, the methods of the present invention provide high 

reliability and reproducibility. Therefore, these methods can be utilized for quantification 
and analysis of many mRNAs in addition to those exemplified here. Thus, this invention 
can be used as a rapid method of analyzing the pathogenesis of various diseases, as well 
as assisting in their diagnosis. Many such methods will suggest themselves to those with 

15 ordinary skill in the art upon a review of the present specification. Accordingly, the scope 

of the present invention is to be determined with reference to the appended claims, and 
is not to be limited to those methods specifically exemplified herein. 
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SEOUENCE LISTING 



(1) GENERAL INFORMATION: 

5 

(i) APPLICANT: Akiiaya. Tatsuo 

Miisuhashi. Masato 
Cooper, AJIan 

10 (ii) TITLE OF INVENTION: METHOD AND REAGENT FOR MEASURING 

MESSENGER RNA 

(in) NUMBER OF SEQUENCES: 711 

15 (iv) CORRESPONDENCE ADDRESS: 

(A) ADDRESSEE: KNOBBE, MARTENS. OLSON. AND BEAR 

(B) STREET: 620 NEWPORT CENTER DRIVE SIXTEENTH FLOOR 

(C) CITY: NEWPORT BEACH 

(D) STATE: CA 

20 (E) COUNTRY: USA 

(F) ZIP: 92660 

(v) COMPUTER READABLE FORM 

(A) MEDIUM TYPE: Floppy disk 

25 (B) COMPUTER: IBM PC compatible 

(C) OPERATING SYSTEM: PC-DOS/MS-DOS 

(D) SOFTWARE: Patent In Release #1.0. Version #1.25 

(vi) CURRENT APPLICATION DATA: 
30 (A) APPLICATION NUMBER- 

(B) FILING DATE: 

(C) CLASSIFICATION: 

(viii) ATTORNEY./AGENT INFORMATION: 
35 (A) NAME: Ahman. Daniel E 

(B) REGISTRATION NUMBER: 34.1 15 

(C) REFERENCE/DOCKET NUMBER: HITACH1.006CP2 

(ix) TELECOMMUNICATION INFORMATION: 
40 (A) TELEPHONE: 714-760-0404 

(B) TELEFAX: 714-760-9502 



45 



(2) INFORMATION FOR SEQ ID NO:l: 



(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 22 bases 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: simile 
50 (D) TOPOLOGY: linear 
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(ii) MOLECULE TYPE: cDNA to mRNA 

(iii) HYPOTHETICAL: NO 

5 (iv) ANTi-SENSE: NO 

(vi) ORIGINAL SOURCE: 

(A) ORGANISM. Gs mRNA NUCLEOTIDE PROBE 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: I: 
TTCATCCTCC CACAGAGCCT TG 
15 (2) INFORMATION FOR SEQ II) NO:2: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 22 bases 

(B) TYPE: nucleic acid 

20 (C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: cDNA to mRNA 
25 (iii) HYPOTHETICAL: NO 

(iv) ANTI-SENSE: NO 

(vi) ORIGINAL SOURCE: 
30 (A) ORGANISM: Gi-2 mRNA NUCLEOTIDE PROBE 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO:2: 

35 ATGGTCAGCC CAGAGCCTCC GG 

(2) INFORMATION FOR SEQ ID NO:3: 

(i) SEQUENCE CHARACTERISTICS: 
40 (A) LENGTH: 22 bases 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: sinele 

(D) TOPOLOGY: linear 

45 (ii) MOLECULE TYPE: cDNA to mRNA 

(iii) HYPOTHETICAL: NO 



50 



(iv) ANTI-SENSE: NO 



WO 93/15221 



PCT/US93/00977 



-105- 

(vi) ORIGINAL SOURCE: 

(A) ORGANISM: G PROTEIN SENSE PCR PRIMER 



5 (xi) SEQUENCE DESCRIPTION: SEQ ID NO:3: 

AG CACCATTG TGAAGCAGAT GA 22 
(2) INFORMATION FOR SEQ ID NO:4: 

10 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 22 bases 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: sincle 
15 (D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: cDNA to mRNA 

(iii) HYPOTHETICAL: NO 

20 

(iv) ANTI-SENSE: YES 

(vi) ORIGINAL SOURCE: 

(A) ORGANISM: G PROTEIN ANT1SENSE PCR PRIMER 

25 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO:4: 
CTCTGGCCTC CCACATCAAA CA tt 

30 

(2) INFORMATION FOR SEQ II) NO:5: 

(i) SEQUENCE CHARACTERISTICS: 
(A) LENGTH: 22 bases 

35 (B) TYPE: nucleic acid 

(C) STRANDEDNESS: sincle 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: cDNA to mRNA 

40 

(iii) HYPOTHETICAL: NO 

(iv) ANTI-SENSE: NO 

45 (vi) ORIGINAL SOURCE: 

(A) ORGANISM: SUBSTANCE P RECEPTOR mRNA PROBE 



50 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO:5 
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GGACTTATGA GAAAGCGTAC CA 22 
(2) INFORMATION FOR SEQ ID NO:6: 

5 (i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 22 bases 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: simile 

(D) TOPOLOGY: linear 

10 

(ii) MOLECULE TYPE: cDNA to mRNA 

(iii) HYPOTHETICAL: NO 

15 (iv) ANTI-SENSE: NO 

(vi) ORIGINAL SOURCE: 

(A) ORGANISM: HUMAN B2 RECEPTOR mRNA PROBE 



20 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO:fi: 
ATGCTGGCCG TGACGCACA GCA 23 
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VVE CLAIM : 

1. A highly sensitive, quantitative and rapid method for delecting and 
quantifying mRNA in a sample without the need to purify mRNA from cells, said method 
comprising: 

5 (a ) identifying a polynucleotide sequence that is unique to said mRNA; 

(b) immobilizing a first polynucleotide to an insoluble support, said first 
polynucleotide having a first sequence complementary to said sequence unique to 
said mRNA; 

(c) incubating said sample with said insoluble support under conditions 
10 wherein said unique sequence will hybridize with said first polynucleotide, thereby 

immobilizing mRNA present in said sample to said insoluble support: 

(d) washing non-immobilized components of said sample from said 
insoluble support; 

(e) labeling mRNA on said support in a manner that label is 
15 incorporated onto said support related to the amount of mRNA on said support; 

and 

(f) measuring the amount of label immobilized on said support. 

2. The method of Claim 1. wherein step (e) comprises: 

i. incubating a labeled probe with said insoluble support under 
20 conditions wherein said labeled probe will hybridize with mRNA on said support, 

said labeled probe bearing a label and being complementary to a part of said 
mRNA other than a part of said mRNA complementary to said first probe, thereby 
immobilizing to said insoluble support the label on said labeled probe which has 
hybridized with said mRNA; and 
25 ii. washing any non-immobilized labeled probe from said insoluble 

support. 

3. The method of Claim 2, wherein said labeled probe comprises poly-d(T). 

4. The method of Claim 1, wherein step (e) comprises labeling 
oligonucleotides on said support with a nucleic acid stain. 

30 5. The method of Claim 4. wherein said nucleic acid slain is selected from the 

group consisting of ethidium bromide, yoyo-1 and toto-1. 

6. The method of Claim 2. wherein said label comprises a radionuclide. 

7. The method of Claim 2, wherein said label comprises biotin. 

8. The method of Claim 7, wherein step (f) comprises; 
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incubation with enzyme-conjugated streptavidin: and 
measurement of the amount of enzyme bound to said insoluble support. 
9. The method of Claim 1. wherein said mRNA codes for a p receptor, and 
wherein said first polynucleotide comprises the sequence of SEQ ID NO:6. 
5 10. The method of Claim 1, wherein said mRNA codes for the a subunit of Gs 

protein, and wherein said first polynucleotide comprises a sequence complementary to a 
sequence selected from the group consisting of SEQ ID NO: 1, 145-146. 209-212, 268-287. 
705 and 70S. 

11. The method of Claim I, wherein said mRNA codes for Gi-2 protein, and 
10 wherein said first polynucleotide comprises a sequence complementary to a sequence 

selected from the group consisting of SEQ ID NO:703. 806. 141 143, 187-197 and 238-253. 

12. The method of Claim 1, wherein said mRNA codes for substance P 
receptor, and wherein said first polynucleotide comprises the sequence of SEQ ID NO:5. 

13. The method of Claim 1, wherein said mRNA codes for Gi-3 protein, and 
15 wherein said first polynucleotide comprises a sequence complementary to a sequence 

selected from the group consisting of SEQ ID NO:7()4. 707. 144. 198-208 and 254-267. 

14. The method of Claim 1, wherein said mRNA codes for Gi-1 protein, and 
wherein said first polynucleotide comprises a sequence complementary to a sequence 
selected from the group consisting of SEQ ID NO:702, 805, 141, and 224-237. 

20 15. The method of Claim 1, wherein said mRNA codes for G protein, and 

wherein said first polynucleotide comprises a sequence complementary to a sequence 
selected from the group consisting of SEQ ID NO:7()0. 701. 139. 140. 148-159. 162-172, 
176-186 and 3-4. 

16. The method of Claim 1, wherein said first sequence hybridizes to said 
25 unique sequence of said mRNA and not to other polynucleotides present in said sample. 

17. The method of Claim 1, wherein said sample comprises a cell lysate. 

18. The method of Claim 2. wherein said label can be measured by light 
emitted therefrom, and wherein step (f) comprises measuring the amount of light emitted 
by said label. 

30 19. The method of Claim 18. wherein said label comprises biotin and wherein 

step (f) additionally comprises adding streptavidin bound to alkaline phosphatase to said 
insoluble support. 

20. The method of Claim 18. wherein the step of measuring the amount of light 
comprises: 
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recording the amount of light on film; and 
measuring the exposure of this film using a densitometer. 
21. The method of Claim 18, additionally comprising adding ATTOPHOS and 
aring flourescense emitted using a fluorimeter. 
5 22. The method of Claim 21, wherein the step of measuring the amount of light 

comprises use of a fluorimeter. 

23. The method of Claim 1, wherein step (a) comprises use of a computer 
program. 

24. The method of Claim 23, wherein step (a) comprises use of an H-site 

10 model. 

25. The method of Claim 24, wherein identification of the polynucleotide 
sequence using said H-site model comprises the steps of: 

specifying a minimum melting temperature for the first polynucleotide and 
said polynucleotide sequence; 
15 specifying a nucleation threshold that places a minimum value on the 

number of base pairs at any nucleation site; 

determining the melting temperatures (T m ) of the first polynucleotide and 
said polynucleotide sequence at every possible hybridization point; and 
selecting the nucleotide sequence having the highest T m value. 
20 26 - The method of Claim 25 wherein said melting temperature is determined by 

the formula: T m = 81.5 - 16.6(log[Na]) - 0.63%(formamide) + 0.41(%(G + C))- 600/N, 
wherein Log[Na] is the log function of the sodium concentration, 0.063 %(formamide) is 
the concentration of formamide, %(G + C) is the percentage of matched GC base pairs, 
and N is the probe length. 

25 27 - The method of Claim 2, wherein the labeled probe is identified using an 

H-site model comprising the steps of: 

specifying a minimum melting temperature for the labeled probe and the 
part of said mRNA complementary to a part of said mRNA other than a part of 
said mRNA complementary to said first probe; 
30 specifying a nucleation threshold that places a minimum value on the 

number of base pairs at any nucleation site; 

determining the melting temperatures (T m ) of the first nucleotide probe 
and said part of said mRNA complementary to a part of said mRNA other than 
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a part of said mRNA complementary to said first probe at every possible 
hybridization point: and 

selecting the nucleotide probe of the proper length having the lowest T m 

value. 

5 28. The method of Claim 27 wherein said melting temperature is determined by 

the formula: T m « 81.5- 16.6(Iog[Na]) - 0.63(%(formamide)) + 0.41(%(G + C))-600/N, 
wherein log[Na] is the log function of the sodium concentration, 0.063%(formamide) is the 
concentration of formamide, %(G + C) is the percentage of matched GC base pairs, and 
N is the probe length. 

10 29. A method of inhibiting RNase activity in a sample containing RNA in a 

buffer system containing SDS and EDTA. comprising adding vanadyl ribonucleoside 

complex to said sample. 

30. The method of Claim 29, wherein the sample contains at least a five-fold 

molar excess of EDTA over RNA in the sample. 
15 31. The method of Claim 29, wherein the vanadyl ribonucleoside complex is 

added in combination with proteinase K. 

32. The method of either Claim 29 or 31. wherein the vanadyl ribonucleoside 
complex is added in combination with RNasin. 

33. A polynucleotide-immobilized support that is -useful in detecting and 
20 quantifying the amount of p receptor mRNA in a sample, comprising: 

an insoluble support; 

a first polynucleotide bound to said support, said first polynucleotide 
comprising a # sequence that is complementary to a sequence unique to p receptor 
mRNA. 

25 34. A polynucleotide-immobilized support as in Claim 33, wherein said 

sequence unique to P receptor mRNA comprises the sequence of SEO ID NO:6. 

35. A polynucleotide-immobilized support that is useful in detecting and 
quantifying the amount of jun oncogene mRNA in a sample, comprising: 

an insoluble support: 

30 a first polynucleotide bound to said support, said first polynucleotide 

comprising a sequence that is complementary to a sequence unique to jun 
oncogene mRNA. 

36. A polynucleotide-immobilized support as in Claim 35, wherein said 
sequence unique to jun oncogene mRNA comprises a sequence complementary to a 
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sequence selected from the group consisting of SEQ ID NO: 11-35. 39-53, 62*138. 310-342. 
600-606, 614-620, and 628-634. 

37. A polynucleotide-immobilized support that is useful in detecting and 
quantifying the amount of G protein mRNA in a sample, comprising: 

5 an insoluble support; 

a first polynucleotide bound to said support, said first polynucleotide 
comprising a sequence that is complementary to a sequence unique to G protein 
mRNA. 

38. A polynucleotide-immobilized support as in Claim 37. wherein said 
10 sequence unique to G protein mRNA comprises a sequence selected from the group 

consisting of SEQ ID NO:l, 2, 700-709, 805, 806. 139-159, 162-172 and 176-309. 

39. A polynucleotide-immobilized support that is useful in detecting and 
quantifying the amount of Substance P receptor mRNA in a sample, comprising: 

an insoluble support; 

15 a firsl polynucleotide bound to said support, said first polynucleotide 

comprising a sequence that is complementary to a sequence unique to Substance 
P receptor mRNA. 

40. A polynucleotide-immobilized support as in Claim 39, wherein said 
sequence unique to Substance P receptor mRNA comprises the sequence of SEQ ID 

20 NO:5. 
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FIG.2 6A 



PROBE: C:\HITACHI\JUNMTX.PRP 
HYBRIDIZATION: C : \HITACHI\HUMBJUNX . CDS 
Length - 374 Hairpin =35 
Locus Pos Tm 

humbjunx 374 61.47 

musbjunx 365 61. 47 

ftumdiunx 41 34 .82 t g-g — agt 

humbjunx 182 31. 12 a gtgg — gc 

huindjunx 602 31. 12 c c-ggg-gc 

huindjunx 602 31. 12 c — -——c-ggg-gc 

PROBE : C : \HITACHI\ JUNMIX . PRP 
HYBRIDIZATION : C : \HITACHI\HUMBJUNX . CDS 
Length » 377 Hairpin = 2 14 
Locus Pos Tn 

humbjunx 377 61.55 

musbjunx 3 68 61.55 — 

huindjunx 383 28. 12 tg-cg-c — g 

musdjunx 38 3 28.12 tg-ca-c — g 

musdjunx 38 3 28.12 tg-ca-c — g 

PROBE: C:\HITACHI\JUNMIX.PRP 
HYBRIDIZATION: C:\HITACHI\HUMBJUNX.CDS 
Length = 389 Hairpin =33 
Locus Pos Tm 

humbjunx 3 89 61.7 — 

muscjunx 314 56.65 -c 

musbjunx 380 50.85 1 — g 

hurncjunx 314 4 9. 35 -t g- 

huindjunx 395 33.85 tt-gc — ag 

musdjunx 395 33.85 tt-gc — aa 

hurncjunx 32 6 32.35 g-ttcgcc— — — — tg 

huindjunx 404 32.35 — ttcgcc — t- 

muscjunx 32 6 32.35 gcttcgcc — — tg 

musdjunx 253 30.85 gacg-gct-ct 

humbjunx 953 30.65 g — t — c-cagct- 

musdjunx 83 27.3 cc-gcggt-gt— — g 

PROBE: C:\HITACHI\JUNMIX.PRP 
HYBRIDIZATION: C : \HITACHI\HUMBJUNX . CDS 
Length =397 Hairpin =41 
Locus Pos Tm 

humbjunx 397 61.55 — — 

muscjunx 322 53.44 — g 

hurncjunx 322 4 5.33 g g 

musbjunx 388 41.38 1 — g 1 

humdjunx 214 3 6.83 ccccctgc 

humdjunx 99 3 6.16 eg — — gc-c 

musdjunx 261 3 4.55 -ct gatct 

humdjunx 4 00 33.27 c ag g 

musdjunx 400 33.27 c ag- a 

hurncjunx 33 4 32.28 — - — — tgcg — c- 

humdjunx 412 32.28 t-a-g-c- 

muscjunx 334 32.28 — tgcg — c- 

humbjunx 658 30. 17 cc-cc -gt 

humdjunx 241 28.95 -c — cacc-c 

humdjunx 34 2 28.95 c-cca-ca ag 

musbjunx 606 28.95 ct-a-ac 

musdjunx 229 28.95 -c-ctgcg-c 

musdjunx 91 26. 67 -gt gcc-ccg 
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PROBE: C:\HITACKI\JUNMIX.PRP 
HYBRIDIZATION : C : \ HITACHI \ HUMBJUNX . CDS 
Length «= 417 Hairpin = 2 15 
Locus Pos Tm 

humbjunx 417 60. 08 

musbjunx 4 08 55. 52 — c 

humdjunx 4 20 37. 3 c — g — — g t-a- 

musbjunx 61 29.0 g — -gg ca-cctgt- 

muscjunx 672 2 6.27 gc-gc a-g — aga — 

PROBE : C : \HITACHI\JUNMIX . PRP 
HYBRIDIZATION : C : \HITACHI\HUMBJUNX . CDS 
Length « 461 Hairpin =49 
Locus Pos Tn 

humbjunx 461 61. 63 — — 

musbjunx 452 61.63 — - 

musbjunx 452 61. 63 — — 

PROBE: C:\HITACHI\JUNMIX.PRP 
HYBRIDIZATION : C:\HITACHI\HUMBJUNX.CDS 
Length « 467 Hairpin « 2 13 



- Locus Pos Tn 

humbjunx 467 61.7 

musbjunx 458 51.6 c-g- 

humdjunx 32 29.35 tgagcgg— ~ gcgg- 

humdjunx 32 29.35 tgagcgg ~ gcgg- 



PROBE: C:\HITACHI\JUNMIX.PRP 
HYBRIDIZATION: C:\HITACHI\HUMBJUNX.CDS 
Length « 477 Hairpin -24 
Locus Pos Tm 

humbjunx 477 61.37 

humdjunx 489 34.93 c-c eg 

humdjunx 4 89 3 4 .93 c-c eg 

PROBE : C:\HITACHI\JUNMIX.PRP 
HYBRIDIZATION: C : \HITACHI \HUMBJUNX . CDS 
Length =4 87 Hairpin =33 
Locus Pos Tm 

humbjunx 487 61.14 

musdjunx 74 51.0 ct 



humdjunx 499 45. 64 1 g 

humdjunx 527 30.72 cc-c-c 

musdjunx 97 30.72 ttc-c -g 

musdjunx 580 30.72 -cc t-g 

musdjunx 637 30.72 cc-cc — g 

musdjunx 637 30.72 cc-cc — g 



PROBE: C:\HITACHI\JUNMIX.PRP 
HYBRIDIZATION : C : \HITACHI\HUMBJUNX . CDS 
Length «= 498 Hairpin = 32 
Locus Pos Tm 

humbjunx 498 61.26 

humbjunx 498 61.26 

PROBE : C : \HITACHI\ JUNMIX . PRP 
HYBRIDIZATION: C:\HITACHI\HUMBJUNX.CDS 
Length « 504 Hairpin =32 



Locus Pos Tm 

humbjunx 504 61.47 

musbjunx 495 40. 35 c — a - 1- 

humdjunx 609 35.29 eg egggg- 

humdjunx 609 3 5.29 eg egggg- 
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AAATGGAACAG CCCTTCTACC 
AATGG AACAG CCCTTCT A CCA 
ATGGAACAGCCCTTCTACCAC 
TCGAACAGCCCTTCTACCACG 
GGAACAGCCCTTCTACCACGA 
G AACAG CCCTTCTACCACG AC 
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ACAGCCCTTCTACCACGACGA 

CAGCCCTTCTACCACGACGAC 

AGCCCTTCTACCACGACGACT 

GCCCTTCTACCACGACGACTC 

CCCTTCTACCACGACGACTCA 

CCTTCTACCACGACGACTCAT 

CTTCTACCACGACGACTCATA 

TTCTACCACGACGACTCATAC 

TCTACCACGACGACTCATACA 

CTACCACGACGACTCATACAC 

TACCACGACGACTCATACACA 

ACCACGACGACTCATACACAG 

CCACGACGACTCATACACAGC 

CACGACGACTCATACACAGCT 

ACG ACGACTCATACACAG CTA 

CGACGACTCATACACAGCTAC 

G ACG ACTCATACACAG CTACG 

ACGACTCATACACAGCTACGG 

CG ACTCATACACAG CTA CGGG 

G ACTCATACACAG CTA CGGG A 

ACTCATACACAGCTACGGGAT 

CTCATACACAGCTACGGGATA 

TCATACACAGCTACG GG ATA C 

CATACACAG CTACGGG ATACG 

ATACACAGCTACGGGATACGG 

TACACAGCTACGGGATACGGC 

ACACAGCTACGGGATACGGCC 

CACAGCTACGGGATACGGCCG 

ACAGCTACGGGATACGGCCGG 

CAGCTACGGGATACGGCCGGG 

AG CTACGGG ATACGGCCGGGC 

G CTACG GGATACGGCCGGGCC 

CTACGGGATACGGCCGGGCCC 

TACGGGATACGGCCGGGCCCC 

ACGGGATACGGCCGGGCCCCT 

CGGGATACGGCCGGGCCCCTG 
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OligcPrcfce. DssignStation 

Probes : C : \HITACHI\HUKBJUNX . CDS 

Preparation : C : \HITACHI\ JUNMIX . PRP 

Locus pos Tm Locus pos Tm Locus pos 



atgtgcactaaaatggaacagcccttctac 
1 30 1 1 1 2 
humbjunx 1 60.7 6 
musbjunx 1 50.03 
muscjunx 1 30.07 
musdjunx 721 27.84 



tgtgcactaaaatggaacagcccttctac 

2 29 1112222234 
humbjunx 65533 60.68 
musbjunx 65533 49.58 
muscjunx 1 29.97 
musdjunx 721 27.66 



gtgcactaaaatggaacagcccttctac 

3 28 1112222234 
humbjunx 65533 60.60 
musbjunx 65533 49.10 
muscjunx 1 29.86 
musdjunx 721 27.47 



tgcactaaaatggaacagcccrttctacc 

4 28 1111222234 
humbjunx 65533 60.60 
musbjunx 65533 46.57 
muscjunx 1 29.86 
musdjunx 729 27.47 



gcactaaaatggaacagcccttctacc 
5 27 1 1 1 
humbjunx 5 60.51 
musbjunx 5 45.96 
muscjunx 1 29.75 
musdjunx 729 27.26 



WO 93/15221 



PCI7US93/00977 



FIG. V5 



Be9in MPSD 







Ccl cu 


late 


output 


<x,y> 


ron 


ges 



500 



Convert to 
I ogcir i thwi c 
scat e 



I n t erpo I ate 



kj<\\ ues 



Create 



b i troop 



■501 



-502 



-503 





504 




D i sp I ay 
b i twip 








End MPSD 

V- 



WO 93/15221 



55/59 



PCI7US93/00977 



FIG. H6> 



Begin Match I nfo 







Fol 1 ow 
cursor i n 
MPSD window 






Update 
pos i t i on of 
candidate probes 






Update sequence & 
hairpin Information 
for candidate 
probes 







Update 
Match 



End Match I nfo 



WO 93/15221 



5 6/59 



PCT/US93/00977 



FIG • 47 



LOCUS HUMBJUNX 1044 bp DNA 19-DEC-1991 

BASE COUNT 195 A 3 68 C 340 G 141 T 

1 ATGTGCACTA AAATGCAACA GCCCTTCTAC CACGACGACT CATACACAGC TACGGGATAC 
61 GGCCGGGCCC CTGGTGGCCT CTCTCTACAC GACTACAAAC TCCTGAAACC GAGCCTGGCG 
121 GTCAACCTGG CCGACCCCTA CCGGAGTCTC AAAGCGCCTG GGGCTCGCGG ACCCGGCCCA 
181 GAGGGCGGCG GTGGCCGCAG CTACTTTTCT GGTCAGGGCT CGGACACCGG CGCGTCTCTC 
241 AAGCTCGCCT CTTCGGAGCT GGAACGCCTG ATTGTCCCCA ACAGCAACGG CGTGATCACG 
301 ACGACGCCTA CACCCCCGGG ACAGTACTTT TACCCCCGCG GGGGTGGCAG CGGTGGAGGT 
3 61 GCAGGGGGCG CAGGGGGCGG CGTCACCGAG GAGCAGGAGG GCTTCGCCGA CGGCTTTGTC 
421 AAAGCCCTGG ACGATCTGCA CAAGATGAAC CACGTGACAC CCCCCAACGT GTCCCTGGGC 
481 GCTACCGGGG GGCCCCCGGC TGGGCCCGGG GGCGTCTACG CCGGCCCGGA GCCACCTCCC 
541 GTTTACACCA ACCTCAGCAG CTACTCCCCA GCCTCTGCGT CCTCGGGAGG CGCCGCGGCT 
601 GCCGTCGGGA CCGGGAGCTC GTACCCGACG ACCACCATCA GCTACCTCCC ACACGCGCCG 
661 CCCTTCGCCG GTGGCCACCC GGCGCAGCTG GGCTTGGGCC GCGGCGCCTC CACCTTCAAG 
721 GAGGAACCGC AGACCGTGCC GGAGGCGCGC AGCCGGGACG CCACGCCGCC GGTGTCCCCC 
781 ATCAACATGG AAGACCAAGA GCGCATCAAA GTGGAGCGCA AGCGGCTGCG GAACCGGCTG 
841 GCGGCCACCA AGTGCCGGAA GCGGAAGCTG GAGCGCATCG CGCGCCTGGA GGACAAGGTG 
901 AAGACGCTCA AGCCCGAGAA CGCGGGGCTG TCGAGTACCG CCGGCCTCCT CCGGGAGCAG 
961 GTGGCCCAGC TCAAACAGAA GGTCATGACC CACGTCAGCA ACGGCTGTCA GCTGCTGCTT 
1021 GGGGTCAAGG GACACCCCTT CTGA 

// 
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LOCUS HUMBJUNX 1044 bp DNA 19-DEC-1991 

3ASE COUNT 195 A 368 C 340 G 141 T 

3RIGIN 

1 ATGTGCACTA AAATGGAACA GCCCTTCTAC CACGACGACT CATACACAGC TACGGGATAC 
61 GGCCGGGCCC CTGGTGGCCT CTCTCTACAC GACTACAAAC TCCTGAAACC GAGCCTGGCG 
121 GTCAACCTGG CCGACCCCTA CCGGAGTCTC AAAGCGCCTG GGGCTCGCGG ACCCGGCCCA 
181 GAGGGCGGCG GTGGCGGCAG CTACTTTTCT GGTCAGGGCT CGGACACCGG CGCGTCTCTC 
241 AAGCTCGCCT CTTCGGAGCT GGAACGCCTG ATTGTCCCCA ACAGCAACGG CGTGATCACG 
301 ACGACGCCTA CACCCCCGGG ACAGTACTTT TACCCCCGCG GGGGTGGCAG CGGTGGAGGT 

3 61 GCAGGGGGCG CAGGGGGCGG CGTCACCGAG GAGCAGGAGG GCTTCGCCGA CGGCTTTGTC 
421 AAAGCCCTGG ACGATCTGCA CAAGATGAAC CACGTGACAC CCCCCAACGT GTCCCTGGGC 
481 GCTACCGGGG GGCCCCCGGC TGGGCCCGGG GGCGTCTACG CCGGCCCGGA GCCACCTCCC 
541 GTTTACACCA ACCTCAGCAG CTACTCCCCA GCCTCTGCGT CCTCGGGAGG CGCCGGGGCT 
601 GCCGTCGGGA CCGGGAGCTC GTACCCGACG ACCACCATCA GCTACCTCCC ACACGCGCCG 
661 CCCTTCGCCG CTGGCCACCC CGCGCAGCTG GGCTTGGGCC GCGGCGCCTC CACCTTCAAG 
721 GAGGAACCGC AGACCGTGCC GGAGGCGCGC AGCCGGGACG CCACGCCGCC GGTGTCCCCC 
781 ATCAACATGG AAGACCAAGA GCGCATCAAA GTGGAGCGCA AGCGGCTGCG GAACCGGCTG 
841 GCGGCCACCA AGTGCCGGAA GCGGAAGCTG GAGCGCATCG CGCGCCTGGA GGACAAGGTG 
901 AAGACGCTCA AGGCCGAGAA CGCGGGGCTG TCGAGTACCG CCGGCCTCCT CCGGGAGCAG 
961 GTGGCCCAGC TCAAACAGAA GGTCATGACC CACGTCAGCA ACGGCTGTCA GCTGCTGCTT 

1021 GGGGTCAAGG GACACGCCTT CTGA 

/ 

.OCUS HUMCJUNX 996 bp DNA 19-DEC-1991 

ASE COUNT 226 A 342 C 299 G 129 T 

RIG IN 

1 ATGACTGCAA AGATGGAAAC GACCTTCTAT GACGATGCCC TCAACGCCTC GTTCCTCCCG 
61 TCCGAGAGCG GACCTTATGG CTACAGTAAC CCCAAGATCC TGAAACAGAG CATGACCCTG 
121 AACCTGGCCG ACCCAGTGGG GAGCCTGAAG CCGCACCTCC GCGCCAAGAA CTCGCACCTC 
181 CTCACCTCGC CCGACGTGGG GCTGCTCAAG CTGGCGTCGC CCGAGCTGGA GCGCCTGATA 
241 ATCCAGTCCA GCAACGGGCA CATCACCACC ACGCCGACCC CCACCCAGTT CCTGTGCCCC 
301 AAGAACGTGA CAGATGAGCA GGAGGGGTTC GCCGAGGGCT TCGTCCGCGC CCTGGCCGAA 
361 CTGCACAGCC AGAACACGCT GCCCAGCGTC ACGTCGGCGG CGC^GCCGGT CAACGGGGCA 
421 GGCATGGTGG CTCCCGCGGT AGCCTCGGTG GCAGGGGGCA GCGGCAGCGG CGGCTTCAGC 

4 81 GCCAGCCTGC ACAGCGAGCC GCCGCTCTAC GCAAACCTCA GCAACTTCAA CCCAGGCGCG 
541 CTGAGCAGCG GCGGCGGGGC GCCCTCCTAC GGCGCGGCCG GCCTGGCCTT TCCCGCGCAA 
601 CCCCAGCAGC AGCAGCAGCC GCCGCACCAC CTGCCCCAGC AGATGCCCGT GCAGCACCCG 
661 CGGCTGCAGG CCCTGAAGGA GGAGCCTCAG ACAGTGCCCG AGATGCCCGG CGAGACACCG 
721 CCCCTGTCCC CCATCGACAT GGAGTCCCAG GAGCGGATCA AGGCGGAGAG GAAGCGCATG 
781 AGGAACCGCA TCGCTGCCTC CAAGTGCCGA AAAAGGAAGC TGGAGAGAAT CCCCCGGCTG 
841 GAGGAAAAAG TGAAAACCTT GAAAGCTCAG AACTCGGAGC TGGCGTCCAC GGCCAACATG 
901 CTCAGGGAAC AGGTGGCACA GCTTAAACAG AAAG TCATGA ACCACGTTAA CAGTGGGTGC 
961 CAACTCATGC TAACGCAGCA GTTGCAAACA TTTTGA 



3CUS 

■FINITION 
SESSION 

:ywords 

JURCE 
ORGANISM 



:ference 

AUTHORS 
JOURNAL 
STANDARD 

'.FERENCE 



HUMDJUNX 1044 bp ss-mRNA PRI 24-MAY-1991 

Hunan junD mRNA 

X56681 

jun-D gene; oncogene. 
Homo sapiens RNA. 
Homo sapiens 

Eukaryota; Animalia; Metazoa; Chordata; Vertebrata? Mammalia; 
Theria: Eutheria; Primates? Haplorhini; Catarrhinir Hominidae. 

1 (bases 1 to 1891) 
Shaul , Y. 

Unpublished (199 0) 
full automatic 

2 (sites) 
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FIG • 48 (2) 



AUTHORS Berger,!. and Shaul,Y. 
TITLE Structure and function of human jun-D 
JOURNAL Unpublished (19 90) 
STANDARD full staff_review 
COMMENT From EMBL 26 entry HSJUKDR; dated 18 -MAR- 1991. 

FEATURES Location/Qualifiers 
mRNA 1..1891 

/gene=" junD" 
/ ev i dence-EXPERIMENTAL 
CDS 175.. 1218 

/product^" junD protein" 
/gene= n junD" 
/codon_start=l 
polyA site 1891. .1891 

lASE COUNT" 162 A 405 C 360 G 117 T 

fRIGIN 

1 ATGGAAACAC CCTTCTACGG CGATGAGGCG CTGAGCGGCC TGGCCGGCGG CGCCAGTGGC 
61 AGCGGCGGCA CGTTCGCGTC CCCGGGCCGC TTGTTCCCCG GGGCGCCCCC GACGGCCGCG 
121 GCCGGCAGCA TGATGAAGAA GGACGCGCTG ACGCTGAGCC TGAGTGAGCA GGTGGCGGCA 
181 GCGCTCAAGC CTGCGCCCGC GCCCGCCTCC TACCCCCCTG CCGCCGACGG CGCCCCCAGC 
241 GCGCCACCCC CCCACGGCCT GCTCGCCTCT CCCGACCTGG GGCTGCTGAA GCTGGCCTCC 
301 CCCGAGCTCG AGCGCCTCAT CATCCAGTCC AACGGGCTGG TCACCACCAC GCCGACGAGC 
3 61 TCACAGTTCC TCTACCCCAA GGTGGCGGCC AGCGAGGAGC AGGAGTTCGC CGAGGGCTTC 
421 GTCAAGGCCC TGGAGGATTT ACACAAGCAG AACCAGCTCG GCGCGGGCCG GGCCGCTGCC 
481 GCCGCCGCCG CCGCCGCCGG GGGGCCCTCG GGCACGGCCA CGGGCTCCGC GCCCCCCGGC 
541 GAGCTGGCCC CGGCGGCGGC CGCGCCCGAA GCGCCTGTCT ACGCGAACCT GAGCAGCTAC 
601 GCGGGCGGCG CCGGGGGCGC GGGGGGCGCC GCGACGGTCG CCTTCGCTGC CGAACCTGTG 
661 CCCTTCCCGC CGCCGCCACC CCCAGGCGCG TTGGGGCCGC CGCGCCTGGC TGCGCTCAAC 
721 GACGAGCCAC AGACGGTGCC CGACGTGCCG AGCTTCGGCG AGAGCCCGCC GTTGTCGCCC 
781 ATCGACATGG ACACGCAGGA GCGCATCAAG GCGGAGCGCA AGCGGCTGCG CAACCGCATC 
341 GCCGCCTCCA AGTGCCGCAA GCGCAAGCTG GAGCGCATCT CGCGCCTGGA AGAGAAAGTG 
901 AAGACCCTCA AGAGTCAGAA CACGGAGCTG GCGTCCACGG CGAGCCTGCT GCGCGAGCAG 
961 GTGGCGCAGC TCAAGCAGAA AGTCCTCAGC CACGTCAACA GCGOCTGCCA GCTGCTGCCC 
1021 CAGCACCAGG TCCCGGCGTA CTGA 



3CUS MUSBJUNX 103 5 bp DNA 19-DEC-1991 

\SE COUNT 210 A 333 C 333 G 159 T 

IIGIN 

1 ATGTGCACGA AAATCCAACA GCCTTTCTAT CACGACGACT CTTACGCAGC GGCGGGATAC 
61 GGTCGGAGCC CTGGCAGCCT GTCTCTACAC GACTACAAAC TCCTGAAACC CACCTTGCCG 
121 CTCAACCTGG CGGATCCCTA TCGGGGTCTC AAGGGTCCTG GGGCGCGGGG TCCAGGCCCG 
181 GAGGGCAGTG GGGCAGGCAG CTACTTTTCG GGTCAGGGAT CAGACACAGG CGCATCTCTG 
241 AAGCTAGCCT CCACGGAACT GGAGCGCTTG ATCGTCCCCA ACAGCAACGG CGTGATCACG 
3 01 ACGACGCCCA CGCCTCCGGG ACAGTACTTT TACCCCCGTG GGGGTGGCAG CGGTGGAGGT 

3 61 ACAGGGGGCG GCGTCACCGA GGAGCAGGAG GGCTTTGCGG ACGGTTTTGT CAAAGCCCTG 
421 GACGACCTGC ACAAGATGAA CCACGTGACG CCCCCCAACG TGTCCCTGGG CGCCAGCGGG 

4 81 GGTCCCCAGG CCGGCCCAGG GGGCGTCTAT GCTGGTCCGG AGCCGCCTCC CGTCTACACC 
541 AACCTCAGCA GTTACTCTCC AGCCTCTGCA CCCTCTGGAG GCTCCGGGAC CGCCGTCGGG 
601 ACTGGGAGCT CATACCCGAC GGCCACCATC AGCTACCTCC CACATGCACC ACCCTTTGCG 
661 GGCGGCCACC CGGCACAGCT GGGTTTGAGT CGTGGCGCTT CCGCCTTTAA AGAGGAACCG 
721 CAGACCGTAC CGGAGGCACG CAGCCGCGAC GCCACGCCGC CTGTGTCCCC CATCAACATG 
781 GAAGACCAGG AGCGCATCAA AGTGGAGCGA AAGCGGCTGC GGAACAGGCT GGCGGCCACC 
841 AAGTGCCGGA AGCGGAAGCT GGAGCGCATC GCGCGCCTGG AGGACAAGGT GAAGACACTC 
901 AAGGCTGAGA ACGCGGGGCT GTCGAGTGCT GCCGGTCTCC TAAGGGAGCA AGTGGCGCAG 
961 CTCAAGCAGA AGGTCATGAC CCATGTCAGC AACGGCTGCC AGTTGCTGCT AGGGGTCAAG 

1021 GGACACGCCT TCTGA 
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7IG. 48 (3) 



w~ ml* 19-DEC-1991 
ccas HUSCJUNX 1005 bp M| ? 

ASE COUNT 223 A 3" 

RICIK » irvrefiAAAC GACCTTCTAC GACGATGCCC TCAACCCCTC GTTCCTCCAG 

1 ATGACTCCAA AGATGGXXAC GA CCTAAGATCC TAAAACAGAG CATGACCTTG 
61 TCCCAGAGCG STCCCTACCG CC CCACCTCC GCGCCAAGAA CTCGGACCTT 

121 AACCTGGCCG ACCCGGTGGG CACTCXGAMj cc CGGAGCTGGA GCGCCTGATC 

181 CTCACGTCGC CCGACGTCGG ACACCGACCC CCACCCAGTT CTTGTGCCCC 

241 ATCCAGTCCA GCAATGGGCA fAJ^g^ GCCGAGGGCT TCCTGCGCGC CCTGGCTGAA 
301 AAGAACGTGA CCGACGAGCA GSAGGGCTTC ^CXCCCST! CAGCGGGGCG 

361 CTGCATAGCC ^AACACGCT GCAGGCGCTG GCGGCCGTGG TCGCTACAGC 

421 GGCATGGTGG CTCCCCCGGT J^ggST* GCCAACCTCA GCAACTTCAA CCCGGGTGCG 
481 GCCAGCCTGC ACAGTGAGCC TCCGGTCTAC CGCTGGCCTT TCCCTCCCAG 

541 CTGAGCAGCG GCGGTGGGGC GCCCTCCTAT "ggg^r jQcccokAC* GATCCCCGTG 
601 CCGCAGCAGC XGCAGCAGCC Jf""^ SgotSgA CCGTGCCGGA CATGCCCGGA 
661 CAG»CCCGC GGCTGCAAGC CJJgJ£J£ £gTCTCAGG AGCGGATCAA GGCAGACAGC 
721 GAGACGCCGC CCCTGTCCCC ™*2£=£ AAGTGCCGGA AAAGGAAGCT GGAGCCGATC 
781 ""SSES 32555 a^SgSSS ACXCCGAGCT CGCATCCACG 

!£ gSacSgI ££S£KS SraSSa ctpaagcaga axgicatgaa ccaccttaac 
\a SSSSS IStcatgct aacgcagcag ttgcaaacgt tttga 



h „ nMA 19-DEC-1991 
TO HUSDJUKX 1026 bp DMA 

,E COUNT 172 A 382 C 

GIN rexeCAGGCG CTGAGCGCCC TGGCTGCGGG TGCGTCGAGC 

1 ATGGAAACGC CCTTCTATCG ggggg SgTGGCTTCG CGCCCCCGGG CCGCGCTTTC 
61 GTCGCTGGTG ScCATGCTG AAGAAAGACG C6CTCACGCT CAGCCTGGCG 

121 CCCGGGGCGC MCCSMSAG SaACCACCC TCGCCCACXG CACCTTCTGC GCTGCGCCCC 
181 GAGCAGGGAG CGGCGGGATT GAAACCACCG Jgf GGCTGCTCAA ACTCGCSTCG 

241 GACGGCGCCC CCfiACGGGCT S5K§£ AACGGCCTGG TGACCACTAC CCCGACCAGT 
301 CCGGAGCTGG ^GCTGAT Cggjgg JSSaGGAGC AGGAGTTCGC CGAAGGCTTC 
361 ACGCAGTTCC TCTACCCGAA GCTCGCAGCC GTGCGGCCAC CGCGGCCACC 

421 GTCAAGGCGC TGGAGGACCT CCACAACCAA CCACCCCCGC GGCCACGGAG 

481 TCAGCCGCTC CCCCGCCTCC CCGGGTCGCG CCGGGCCCCC TGGGGGCGCG 

541 ACCCCCGTCT ACCCCAACCT GAGCACTTTC G ^ CGCGCTGGGG 

601 GCCACCGTGG CTTTCGCCGC JGAGCCACTG A6GACGAGCC GCAGACCGTG 

661 CCGCCGCCAC CTCCGCATCC ^CGCGCCTG ^ CCATCQACAT GGACACGCAA 

721 CCGGACGTGC CGAGCTTCGG CGA»«gJ CGCAACCGCA TCOCCGCCTC CAAATGCCGC 
781 GAACGCATCA AGGCGGAGCG CAAQAGGCTG ^ TCXAGACCCT CAAAAGCCAG 

841 AAGCGCAAGC TGGAGCGTAT ^GCG^ ctgcGCGAGC AGGTGGCGCA GCTCAAACAG 
|°6i gcc™ SgSgSgI CAGCTGCTGC CCCAGCACCA GGTCCCGGCG 

1021 TACTGA 
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